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Preface 

This  topic  was  researched  at  the  request  of  the  Operations 
Evaluation  Group,  Assistant  Chief  of  Staff,  Studies  and  Analysis,  USAF. 
The  performance  of  a  laser  guided  stand-off  missile  in  the  close-atr- 
support  mission  is  dependent  on  a  complex  interdependent  set  of  factors 
including  laser  design,  target  characteristics,  receiver  design,  and 
atmospheric  conditions.  Although  the  relationship  among  all  these 
factors  is  developed,  this  work  addresses  in  detail  only  the  problem 
of  atmospheric  transmission  in  determining  maximum  missile  lock-on 
range.  An  applied  approach  is  taken  in  which  factors  afiecting 
atmospheric  transmission  are  discussed,  these  factors  are  then  related 
to  atmospheric  observables,  and  these  observables  are  then  related  to 
an  atmospheric  transmission.  Art  attempt  is  made  to  utilize  existing 
atmospheric  models  and  widely  accepted  concepts  to  apply  to  this 
Important  operational  problem.  Where  model  theories  are  not  supported 
by  measurements,  sensitivity  analyses  are  undertaken.  It  has  become 
obvious  to  the  author  in  the  course  of  this  investigation  that  the 
model  predictions  when  compared  to  actual  measurements  are  still 
Inexact  due  to  the  inaccuracies  in  the  measurement  and  reporting  of 
meteorological  conditions,  the  difficulty  cf  measuring  aerosol 
composition,  size,  and  distribution,  and  the  many  approximations  which 
must  be  made  to  reasonably  calculate  the  atmospheric  transmission. 

I  would  like  to  express  my  sincere  appreciation  to  the  many 
people  without  'whose  assistance  such  an  indepth  investigation  could 
never  have  been  accomplished  in  the  aiiotted  time  period.  In  particular, 


,”T^' 


1,1  J  »,U  JAfU  W 


r,7Frzr‘ 


GEP/PH/7**-l» 


I  would  like  to  thank  Colonel  Ed  Battle  who  sponsored  this  work  and 
who  provided  cc-icinued  support  throughout  the  entire  effort.  A  spe¬ 
cial  debt  of  gratitude  is  owed  to  Major  Paul  Try  of  Headquarters, 

Air  Weather  Service,  for  his  hours  of  assistance  and  consultation 
throughout  this  project  and  in  particular  for  his  providing  the  initial 
Mie  Theory  computer  program  and  expert  advice  as  to  changes  that  could 
be  made  in  the  program.  Many  interesting  discussions  on  various 
approaches  to  the  problem  were  held  with  Captain  Neil  McQuage,  Staff 
Meteorologist,  Air  Force  Avionics  Laboratory.  Mr.  Paul  Markley  and 
Mr.  Phil  Miller  of  Aeronautical  Systems  Division  and  Mr.  Frank  McCann 
of  the  Avionics  Laboratory  provided  valuable  literature  sources. 

f  also  thank  Doctors  Louis  Elterman,  Robert  Fean,  Robert 
McClatchey,  Erik  Shettle,  and  Frederic.  Volz  of  Air  Force  Cambridge 
Research  Laboratory  for  their  expert  guidance  and  the  valuable  infor¬ 
mation  they  provided.  Without  their  detailed  and  candid  advice  about 
shortcomings  of  present  models  and  without  results  of  the  most  recent 
of  their  investigations  (some  of  it  unpublished),  this  work  would  have 
been  considerably  limited.  In  addition,  appreciation  goes  to  the 
following  at  Air  Force  Environmental  Technical  Applications  Center  who 
provided  climatological  data  and  advice  for  this  project:  Majors  Tom 
Stanton  and  David  Lydon,  Captain  Ken  Nozaki,  and  Lt  Gerald  Riley. 

Major  Carl  Case,  my  thesis  advisor  at  Air  Force  Institute  of  Technology, 
provided  valuable  insight  and  direction  that  enabled  me  to  clarify  the 
problem.  My  wife,  Bonnie,  and  son,  Charles,  have  given  their  unselfish 


support  and  understanding  throughout 


the  days  necessary  io  compiete 


this  work. 


Charles  H.  Coolidge,  Jr. 
Captain,  USAF 
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Abstract 

Laser  designator  characteristics,  target  composition  and  size, 
laser  designator-target-receiver  geometry,  receiver  characteristics, 
and  atmospheric  conditions  limit  the  maximum  range  at  which  a  laser 
target  designator  system  may  lock-on  and  track  a  target.  This  report 
shows  the  relationship  of  these  factors  and  develops  a  model  for  the 
prediction  of  maximum  lock-on  range  for  a  1 .06y  laser  target  designator 
system  where  atmospheric  attenuation  is  the  limiting  system  factor. 

A  transcendental  lock-on  range  equation  is  developed  which  is 
a  function  only  of  atmospheric  conditions  and  the  designator-target- 
receiver  geometry.  Several  qualifying  assumptions  are  made  in  the 
derivation  of  this  equation:  (1)  the  typical  military  target  is  a 
diffuse  Lambertian  surface  with  ten  percent  reflectivity;  (2)  the 
laser  energy  is  reflected  only  from  the  target  surface;  and  (3)  the 
reflected  radiation  is  entirely  within  the  field  of  view  of  the 
receiver. 

This  report  assumes  the  validity  of  the  Beer-Lambert  Law  for 
describing  atmospheric  transmission  of  laser  radiation.  The  deter¬ 
mination  of  atmospheric  attenuation  is  reduced  to  a  calculation  of  the 
attenuation  due  to  aerosol  scattering  and  absorption,  molecular 
scattering,  molecular  absorption,  and  scattering  by  fogs,  clouds,  and 
rain.  Theoretical  calculations  of  the  attenuation  coefficients 
demonstrate  that  fogs,  clouds,  and  rain  are  the  largest  atmospheric 
attenuators.  For  dry  atmospheric  conditions,  the  chief  attenuating 
mechanism  of  1.06y  laser  radiation  is  the  aerosol.  Molecular 
absorption  proves  to  be  of  lesser  importance  while  molecular  (Rayleigh) 
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scattering  is  negligible. 

Theoretical  calculations  of  molecular  absorption  coefficients 
are  made  using  atmospheric  absorption  line  intensities  compiled  by 
Air  Force  Cambridge  Research  Laboratories.  Molecular  absorption 
coefficients,  which  are  a  function  of  temperature,  pressure  and 
relative  humidity,  are  related  to  mean  seasonal  atmospheric  profiles. 
Calculations  of  the  aerosol  attenuation  coefficient  are  made  using 
Mie  Theory;  these  calculations  are  limited  to  single  .scattering  theory. 
The  calculations  are  based  on  a  continental  aerosol  model  using  re~ 
fractive  indices  typical  of  those  reported  for  Central  Europe,  although 
sensitivity  analyses  are  conduced  using  other  models.  Aerosol 
attenuation  is  related  to  ground  level  horizontal  visibility.  A 
boundary  (mixed)  layer  concept  Is  used  to  describe  the  vertical 
distribution  of  aerosols. 

Mie  Theory  calculations  are  also  made  to  determine  the 
attenuation  coefficient  for  clouds  where  clouds  are  assumed  to  be 
homogeneous.  Attenuation  coefficients  for  rain  are  calculated  using 
the  geometrical  approximation  for  Mie  Theory  where  rain  showers  are 
assumed  to  be  homogeneous  and  monodlsperse.  Rain  attenuation 
coefficients  are  related  to  rain  rates  reported  by  the  meteorologist. 
Attenuation  coefficients  for  fog  are  related  to  visibility  and  are 
those  used  in  a  previous  Monte  Carlo  study. 

Equations  and  Tables  are  provided  for  ease  in  determining 
attenuation  coefficients  as  a  function  of  meteorological  observables. 
These  coefficients  may  then  be  used  to  solve  the  maximum  lock-on  range 
equation  for  various  tactical  scenarios.  The  model  is  specifically 
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designed  for  use  in  variable  laser-target-receiver  geometry  and  slant 
path  calculations. 


xii 


gep/phM-^ 


ATMOSPHERIC  TRANSMISSION  OF  1 .06  MICRON  LASER  RADIATION: 
APPLICATION  TO  STAND-OFF  MISSILE  PERFORMANCE 

!.  Introduction 


Background 

Laser  target  designator  systems  locate  targets  using  the  same 
basic  principles  used  by  conventional  radar;  that  is,  the  system  de¬ 
tects  radiation  reflected  from  an  object  (target).  The  purpose  of  the 
laser  target  designator  system,  however,  is  to  designate  where  the 
target  is  located  relative  to  the  receiver  rather  than  to  measure  the 
range.  The  laser  transmitter,  also  called  designator  or  illuminator, 
illuminates  the  target  with  pulses  of  radiant  energy  in  a  narrow  beam 
and  the  receiver  detector  discerns  (in  its  relatively  wide  field  of 
view)  where  the  laser  energy  is  being  reflected,  thereby  locating  the 
target  (Ref  2*»:13-6). 

The  receiver  detector  for  the  laser  target  designator  system 
detects  the  narrow-band  reflected  laser  radiation.  Once  the  radiation 
is  detected  by  the  receiver  sensor,  then  tracking  is  accomplished  by 
one  of  two  methods.  The  receiver  may  track  the  "one  brightest  spot  in 
the  field  of  view"  or  may  track  the  "centroid  of  entire  illumination 
in  the  field  of  view"  (Ref  52:23).  This  detection  and  tracking  scheme 
of  the  laser  designator  receiver  system  is  then  unlike  the  broad-band 
detection  systems  which  detect  target-background  contrast.  These 
broad-band  detection  systems  generally  use  a  high-resoluticn  TV  imaging 


system  which  provides  a  target  Image  on  the  TV'  screen.  The  system 


operator  then  must  place  a  tracking  gate  over  the  display  to  begin 


I 


ssf 


v*  , 1  „,  j'.'.a, j  nnmn^ngnpiw  wumuwj 


GEP/PH/7V-A 


receiver  tracking. 

The  laser  target  designator  system,  on  the  other  hand,  requires 
that  the  target  be  initially  acquired  by  seme  other  means.  This  Im¬ 
plies  that  the  target  must  have  been  recognized  by  the  system  operator 
using  some  means  (usually  visual)  other  than  the  designator  system. 

The  laser  designator  may  then  be  directed  at  the  target  (object).  Once 
the  laser  Is  illuminating  the  target,  then  receiver  lock-on  and  tracking 
Is  possible. 

Several  factors  determine  the  maximum  range  at  which  the 
receiver  may  acquire  and  track  a  target  once  the  target  has  been 
Illuminated.  These  factors  include:  detector  sensitivity  and  design, 
target  geometry  and  size,  laser  transmitter  (illuminator)  power  and 
beam  geometry  and  positions  of  the  receiver  and  illuminator  relative 
to  the  target.  The  last  factor  is  complicated  by  the  fact  that  the 
target  must  be  initially  acquired  by  means  other  than  the  designator 
system.  Because  the  laser  target  designator  requires  acquisition  of 
the  target  by  other  means,  many  tactical  system  uses  have  evolved. 
Possible  uses  include:  (1)  the  receiver  and  transmitter  are  co-located 
on  the  same  aircraft;  (2)  the  receiver  and  transmitter  are  located  on 
separate  aircraft;  and  (3)  the  transmitter  is  ground-based,  but  the 
receiver  is  an  aircraft.  Thus,  the  laser  radiation  may  traverse  two 
completely  different  paths.  This  report  addresses  this  general  two 
path  problem.  The  co-location  of  source  and  receiver  then  is  a 
simplification  of  this  more  general  case. 

An  additional  complicating  factor  is  atmospheric  attenuation 
of  the  laser  radiation.  The  laser  radiation  is  scattered  and  absorbed 
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by  atmospheric  aerosols  and  molecules.  Other  atmospheric  phenomena 
such  as  clouds,  fog,  and  precipitation  also  scatter  and  absorb  laser 
radiation.  In  fact,  this  atmospheric  attenuation  is  a  serious  system 
limiting  factor  in  the  effective  operational  use  of  a  laser  target 
designator  system. 

The  determinat'on  of  atmospheric  effects  on  the  laser  radiation 
may  be  reduced  to  a  determination  of  the  effects  of  the  separate 
attenuating  mechanisms.  These  mechanisms — molecular  scattering, 
molecular  absorption  and  aerosol  scattering  and  absorption — are  then 
related  to  meteorological  observables.  Using  these  relationships,  one 
can  determine  the  total  atmospheric  attenuation  of  radiation  as  a 
function  of  radiation  path  and  the  meteorological  observables.  Similar 
relationships  may  be  used  to  determine  the  attenuation  of  laser 
radiation  due  to  the  other  attenuating  mechanisms:  clouds,  fog,  and 
precipitation. 

Purpose  of  the  Report 

The  purpose  of  this  report  is  to  investigate  the  atmospheric 
attenuation  of  laser  radiation  and  to  determine  methods  for  prediction 
of  laser  target  designator  system  lock-cn  ranges.  This  involves 
Investigating  the  various  attenuation  factors  and  then  associating 
these  factors  with  measurable  atmospheric  quantities.  This  report 
includes  the  investigation  of  several  different  atmospheric  attenuation 
models  to  determine  the  assumptions  of  each.  The  modeis  also  are 
compared  with  each  other  and  with  results  of  flight  test  data. 
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Scope  of  Report 

Since  atmospheric  attenuation  of  radiation  is  sensitive  to 
the  wavelength  of  the  laser  radiation,  then  it  is  necessary  to  limit 
the  frequencies  for  which  the  atmospheric  attenuation  is  determined. 

This  report,  therefore,  investigates  the  attenuation  of  iaser  radiation 
for  one  wavelength,  1.06p.  In  addition,  this  investigation  is  limited 
to  atmospheric  attenuation  models  which  predict  laser  target  designator 
system  performance  in  a  continental  geographical  area  such  as  Germany, 
it  is  also  limited  to  vertical  models  of  the  lowest  20,000  feet  of  the 
troposphere. 

The  report  compares  various  models  for  the  prediction  of 
atmospheric  attenuation.  An  attempt  is  made  to  utilize  existing  models 
whenever  possible.  Where  existing  models  are  unsatisfactory  for  this 
particular  tactical  application,  then  they  are  modified  to  conform  with 
recent  measurements.  The  report  discusses  in  detail  the  attenuating 
mechanisms  and  the  relationship  of  these  mechanisms  to  the  attenuation 
models. 

it  does  not,  however,  examine  the  details  of  the  various  sensor 
systems.  For  an  assumed  sensor  design,  this  report  determines  the 
maximum  system  lock-on  range  where  atmospheric  attenuation  is  the 
limiting  system  factor.  In  addition,  two  atmospheric  related  factors 
not  included  in  this  investigation  are:  (l)  beamspreading  due  to 
turbulence;  and  (2)  change  in  the  refractive  index  of  air  with  change 
in  altitude. 

The  first  of  the  two  factors  is  not  investigated  because  several 
references  indicated  that  small  pockets  of  air  rising  and  falling  in 
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the  path  of  the  beam  would  cause  local  defocusing  of  the  beam  but  that 
the  energy  reflected  by  the  target  will  not  fluctuate  significantly 
(Ref  6  and  Ref  52).  For  the  tactical  problem  being  considered,  then, 
turbulence  is  a  negligible  factor;  for  other  systems,  especially 
communication  systems,  turbulence  is  not  necessarily  negligible 
(Ref  56:B-9). 

The  change  in  the  refractive  index  of  air  with  altitude  is  not 
discussed  due  to  the  negligible  effect  that  this  refraction  has  on 
laser  radiation  at  these  altitudes.  This  factor  only  becomes  signifi¬ 
cant  for  long  slant  paths  when  the  angle  between  the  receiver  and  the 
vertical  is  larger  than  80°  (Ref  14:^2). 


Assumptions 

When  modeling  something  as  complex  as  the  atmosphere,  certain 
assumptions  must  be  made.  These  assumptions  are  too  numerous  and  too 
lengthy  to  enumerate  at  this  point.  The  assumptions  will  be  pointed 
out  in  the  discussion  of  each  of  the  separate  models.  Two  important 
assumption*,  however,  are  worth  mention.  One  is  that  the  laser 
radiation  is  assumed  to  be  monochromatic.  The  other  is  that  the  beam 
geometry  is  neglected  and  the  attenuation  of  radiation  in  the  atmos¬ 
phere  is  then  a  function  only  of  the  attenuating  mechanisms. 

Development  of  the  Report 

In  Chapter  II,  the  variables  that  determine  maximum  lock-on 
range  in  an  atmosphere-free  environment  are  discussed.  An  atmosphere- 
free  laser  lock-on  range  equation  is  veloped.  Chapter  111  discusses 
atmospheric  effects  on  laser  radiation  and  modifies  the  laser  lock-on 
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range  equation  to  include  atmospheric  effects,  in  Chapter  IV,  the 
atmospheric  attenuation  of  laser  radiation  due  to  molecular  scattering 
and  absorption  is  discussed  in  detail.  Chapter  V  describes  atmos¬ 
pheric  aerosol  in  general.  This  general  knowledge  is  then  applied  In 
Chapters  VI  and  VII  in  the  development  of  the  model  for  prediction  of 
scattering  and  absorption  due  to  atmospheric  aerosol.  Chapter  VI 
discusses  the  ground  level  aerosol  attenuation  while  Chapter  VII 
investigates  the  vertical  aerosol  distribution. 

Chapter  VI I i  presents  methods  for  predicting  attenuation  in 
clouds,  rain,  and  fog.  The  use  of  the  model  is  described  in  Appendix 
C.  The  conclusions  and  recommendations  are  presented  in  Chapter  IX. 
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I l .  Development  of  the  Atmosphere-Free  Laser 
Lock~on  Range  Equation 

An  understanding  of  the  many  variables  that  determine  the 
maximum  range  at  which  a  laser  guided  receiver  can  acquire  and  track 
the  target  {lock-on  range)  is  essential  prior  to  an  investigation  of 
atmospheric  attenuation.  Excluding  the  atmosphere,  factors  which 
affect  the  maximum  lock-on  range  and  which  must  be  dL.  j  ,<±d  are: 

(1)  laser  designator  characteristics,  (2)  target  size  and  composition, 
and  (3)  receiver  characteristics.  In  tracing  the  path  of  laser  energy 
from  the  laser  to  receiver,  one  must  successively  investigate  the  laser, 
the  target,  and  receiver — all  of  which  may  limit  detection  range. 


Energy  From  Laser.  First,  one  must  know  peak  laser  power, 
pulse  duration,  transmitter  optics,  beam  profile  and  beam  width  in 
order  to  determine  the  power  density  within  the  beam.  If  there  is  no 
Intervening  medium,  then  the  radiant  power  or  radiant  flux  of  the  beam 
is  conserved.  The  radiant  intensity,  the  power  per  unit  solid  angle, 
will  change  due  to  the  change  in  beam  geometry,  and  it  may  not  be  equal 
across  the  beam  to  begin  with  depending  on  whether  beam  profile  is 
gaussian  or  constant  average  distribution.  Nevertheless,  total  power 
Is  conserved.  The  radiant  power  arriving  at  the  target  then  is  the 
peak  laser  power  transmitted  through  the  designator  optics  system 


where  is  the  peak  power  transmitted, 

P^  is  the  peak  power  of  the  laser, 
and  Tj  is  the  transmittance  of  the  designator  optics. 
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Actually,  the  beam  loses  some  of  Its  energy  due  to  scattering  and  ab- 
sorption  by  the  atmosphere.  For  purposes  of  our  present  discussion, 
however,  all  of  the  power  arrives  at  the  target  area. 

Energy  Loss  at  Target,  in  this  simplification,  the  target 
becomes  the  first  loss  of  energy.  Target  size,  geometry,  orientation, 
and  reflective  characteristics  determine  the  total  radiant  power 
reflected  into  the  direction  of  the  receiver.  If  the  beam  size  at  the 
target  site  is  larger  than  the  target,  then  the  problem  becomes  more 
complicated.  In  this  event,  the  beam  will  intersect  the  background 
and  reflect  from  these  surfaces  as  well  as  from  the  target.  This  is  a 
serious  laser  designator  problem  for  now  the  receiver  may  lock-cn  the 
background  in  lieu  of  the  target  depending  on  the  various  compositions 
of  the  target  and  background.  This  problem  is  most  clearly  illustrated 
with  the  following  example.  Suppose  that  the  target  were  a  long 
narrow  boat  on  the  open  water.  The  boat  might  be  narrower  than  the 
beam,  and  therefore  part  of  the  beam  might  "spillover"  into  the  back¬ 
ground  which  in  this  case  is  water.  The  water  being  highly  reflective 
or  at  least  more  reflective  than  the  vessel  might  cause  the  receiver 
to  be  guided  to  the  background,  thereby  missing  the  target.  Although 
this  problem  of  beam  geometry  versus  target  size  is  a  serious  system 
limitation,  i t  w?  1 1  be  assumed  for  the  purpose  of  this  work  that  the 
target  is  larger  than  the  laser  beam  intersecting  it.  The  beam,  then, 
Is  reflected  only  from  the  target. 

The  direction  and  amount  of  this  reflected  radiation  is 
extremely  important  if  the  detector  is  to  receive  as  much  of  it  a s 
possible.  Most  military  targets  do  not  have  highly  reflective  surfaces 
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but  absorb  and  diffuse  the  radiant  laser  energy.  Moreover,  their 
surfaces  are  usually  highly  Irregular,  often  with  edges  and  often  at 
many  different  angles  to  the  Incoming  laser  beam.  Thus,  In  the 
operational  use,  the  amount  and  direction  of  reflected  radiation  will 
vary  widely  depending  on  the  Illuminated  target  surfaces,  edges,  and 
composition  and  the  relative  positions  of  the  designator  and  receiver. 
The  following  illustrations  may  place  this  difficult  problem  In  better 
perspective.  Consider  the  relative  positions  of  the  designator  and 
receiver  with  respect  to  the  thatched  target  as  shown  In  Figure  1. 


Fig.  1.  Receiver  Oriented  at  more  than  90°  from  Perpendicular  to 

Target  Reflecting  Surface 


The  receiver  detects  no  laser  radiation  because  the  detector  Is 
oriented  at  more  than  90  degrees  from  the  perpendicular  to  the  surface 
If  this  same  target  were  oriented  as  In  Figure  2,  then  the  receiver 
might  see  more  of  the  projected  spot  size  than  the  illuminator. 


!f  the  targets  were  not  the  flat  surfaces  as  shown  in  Figures  1  and  2, 
but  were  highly  irregular  surfaces,  then  only  a  small  amount  of  the 
laser  radiation  might  be. reflected  into  the  viewing  angle  of  the 
receiver.  The  problem  is  further  corr.pl  icated  by  relative  positions  of 
the  designator  and  receiver  or  by  their  co-location. 

The  evaluation  of  the  bi-directional  reflectance,  the  ratio  of 
radiant  power  (or  radiant  flux  as  it  is  often  called)  reflected  into  a 
particular  solid  angle  to  the  radiant  power  incident  from  a  specific 
solid  angle,  is  extremely  complex.  This  reflection  is  often  a 
specular  reflection  from  a  hard  surface  and  the  reflected  radiation  is 
usually  orders  of  magnitude  less  in  all  other  directions  (Ref  5).  The 
Target  Signature  Analysis  Center  (TSAC) ,  University  of  Michigan,  has 
investigated  the  bi-directional  reflectance  of  many  materials  of 
tactical  interest.  TSAC's  mathematical  model  for  bi-directional 
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reflectances  consists  of  weighted  sums  of  diffuse  reflection  from  a 
hard  surface,  specular  reflection  from  a  hard  surface  and  multiple 
diffuse  volume  reflection  from  a  thin  volume  of  scattering  dielec¬ 
trics  on  the  surface  (Ref  52:15).  The  three  parts  of  the  model  are 
weighted  according  to  experimental  results.  In  order  to  address  the 
reflectance  problem  at  length,  then  a  maximum  lock-on  range  would  have 
to  be  calculated  for  each  target  of  tactical  importance. 

Although  the  bi-directional  reflectance  problem  is  extremely 
complicated,  certain  assumptions  may  be  made.  TSAC  has  published  the 
1.06  micron  laser  reflective  characteristics  of  targets  of  military 
significance.  Variations  in  reflective  properties  were  found  for  new 
and  old  olive  drab  paints,  surface  conditions  (wet,  dry  or  dusty),  and 
target  composition.  Reflectances  of  2.9%  to  12.1%  from  clean  surfaces 
of  a  military  jeep  were  reported  by  Mardis  (Ref  52:19).  Dusty  sur¬ 
faces  nad  reflectances  as  high  as  60%.  These  targets,  in  addition, 
exhibited  less  specular  reflectance  than  many  other  non-military 
targets  and  background  materials.  For  the  purpose  of  this  report,  a 
reasonable  assumption  for  the  target  reflective  characteristics  is  that 
the  target  is  a  perfectly  diffusing  reflector  of  incident  energy,  a 
Lambertian  emitter  or  reflector,  with  a  reflectance  of  10%.  This 
assumption  compares  favorably  with  those  made  by  others  (Ref  24:13-8). 

Briefly,  a  Lambert  source  emits  radiation  completely  randomly 

and  this  radiation  propagates  in  all  directions,  from  the  source  with 

equal  radiance,  N.  Then  for  the  Lambert  surface,  the  bi-directional 

reflectance,  pt  ,  summed  over  all  solid  angles  yields  the  total 
u 

reflectance,  p^,  given  by 


II 


«  /pbdn  *=  pb/dn 


(2) 


where  dft  is  a  differential  solid  angle.  The  reflectance  into  the 
solid  angle  subtended  by  the  receiver  then  is 


pb 


Pt 

—  cos  e  ft 
IT  t 


(3) 


where  ft  is  the  receiver  solid  angle  and  0^  is  the  angle  between  a  nor¬ 
mal  to  the  reflecting  surface  and  a  line  connecting  the  surface  to  the 
receiver  as  shown  in  Figure  3*  The  solid  angle  subtended  by  the 
receiver  at  a  distance  Rr  is 


w 


where  Ar  is  the  area  of  the  receiver  optics.  The  area  of  the  receiver 

optics,  Ar  =  Trdr?  ,‘where  df  is  the  diameter  of  the  receiver,  is  multi- 
!| 

plied  by  cos  $r.  The  angle  <{>r  is  the  angle  between  the  line  connecting 
the  receiver  and  reflecting  target  surface  and  the  normal  to  the 
receiver  surface  as  shown  in  Figure  3.  The  factor,  cos  tf>r,  then  is 
Included  for  those  occasions  when  the  receiver  is  not  looking  directly 
at  the  target  surface. 

For  the  tactical  application  considered  in  this  report,  the 
receiver,  in  general,  views  the  Lambert  reflective  surface  (target) 
as  a  point  source  (Ref  k2) .  The  receiver,  then  being  a  long  distance 
from  the  target,  effectively  sees  the  entire  surface  and  this  surface 
subtends  a  very  small  angle  in  the  receiver's  relatively  wide  field  of 
view. 
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Combining  Eqs  (l),  (3)  and  (4),  then  the  power  or  radiant  flux 
incident  on  the  receiver  optics,  P'r,  is 


P' 


r 


P  T.p  cos0  A  cosA 
p  d  t  t  r 


r 


(5) 


This  radiant  flux  is  further  reduced  by  the  receiver  optics.  Hence, 
the  radiant  flux  received  by  the  sensor,  P  ,  is  given  by 


p  =  Pp1>tc°s6tArc°s*rTr  (6) 

r  7rRf2 

where  Tr  is  the  transmission  of  the  receiver  optics.  This  expression 
gives  the  power  detected  by  the  receiver  in  an  atmosphere-free 
environment  from  a  laser  of  peak  power  Pp  reflected  from  a  Lambertian 
surface  of  reflectivity  p^..  As  it  has  been  developed,  this  expression 
is  valid  for  the  condition  that  the  laser  is  the  sole  emitter  of  this 
radiation  in  the  viewing  angle  of  the  receiver. 

Unfortunately,  there  are  other  sources  of  1.06y  radiation 
within  the  receiver's  field  of  view.  This  causes  the  receiver  to 
detect  a  background  signal  of  1.06y  radiation.  The  strength  of  this 
background  signal  and  the  design  characteristics  of  the  receiver 
detector  determine  the  minimum  detectable  power  at  the  receiver.  This 
is  usually  described  by  a  signal-to-noise  ratio  (Ref  2k:]3S)  given  by 

b2p«2r*g2  ,7) 
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where  S/N  Is  the  peak  signal -to-rr.'Dt-mean-square  noise  ratio  of  the 
receiver  detector  for  one  pulse, 

B  is  the  bandwidth  of  the  receiver  (Hz), 

3  Is  the  responsivity  of  the  photodetector  (amp/watt), 
e  Is  the  charge  on  an  electron  (1.6  x  10“19  coul), 
r  is  the  noise  factor  of  the  receiver, 

G  is  the  internal  gain  of  the  photodetector, 
is  the  detector  dark  current  (amp), 

K  is  Boltzmann’s  constant  (1.38  x  10"23  J/°K) , 

Pp  is  the  peak  power  received, 

’Pjj  is  the  power  received  from  background  radiation  (watt), 
is  the  detector  load  resistance, 
and  T  Is  the  Absolute  Temperature  of  the  detector  (°K) . 

The  radiant  flux  from  the  background  radiation,  P^,  is  found  to  vary 
depending  on  the  viewing  angle  of  the  detector.  The  reflected  laser 
power  received  is  essentially  monochromatic,  but  the  background  power 
received  is  much  broader  spectrally,  on  the  order  of  100  to  200  Angtroms 
with  current  detector  filters  in  use.  Thus,  we  must  investigate  the 
power  of  background  sources  over  the  broader  spectral  region.  We  shall 
discuss  the  broader  range  using  100  Angstroms  (A)  as  an  example. 

One  background  source  of  1.05SU  to  1.069y  radiation  is  the  sun. 
The  sun  is  obviously  a  background  source  problem  only  during  daylight 
hours,  but  the  problem  is  worthy  of  discussion.  Radiation  from  the 
sun  incident  on  the  earth's  atmosphere  is  scattered  and  absorbed  by  the 
atmospheric  aerosols  and  moiecuies.  The  solar  irradiance  at  the  mean 
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earth-sun  distance  incident  on  the  earth's  atmospheric  layer  is 
approximately  6.4*1  watts/m2  in  the  spectral  region  1  .059u  to  1,Q69p 
(Ref  61:16-6,  1 6-8)  (Ref  58:7)*  The  flux  density  of  the  solar  radi¬ 
ation  at  any  point  in  the  atmosphere  is  a  complicated  function  of  the 
path  traveled  by  that  radiation.  That  is,  the  atmosphere's  radiance 
at  the  point  of  observrtion  (detection)  is  a  function  of  the  sun's 
zenith  angle  and  the  receiver's  altitude  and  azimuthal  and  zenith 
angles  (which  together  determine  the  length  of  the  radiance  path),  the 
number  and  kinds  of  scatterers  in  the  path,  and  the  number  and  kinds  of 
absorbers  in  the  path.  The  geometry  of  this  problem  is  given  in  Figure 
4. 

There  is  no  significant  atmospheric  absorption  of  the  sun's 
broadband  (1  .059u  to  1  .069p)  radiation  (P.ef:48).  In  addition,  the 
emission  of  radiation  by  the  atmosphere,  and  in  particular  by  clouds, 
is  an  important  factor  only  for  detectors  that  utilize  the  electro¬ 
magnetic  spectrum  beyond  approximately  4  microns  (Ref  3:1314).  Thus, 
the  problem  is  reduced  to  determining  how  much  of  the  sun's  radiation 
is  incident  on  the  detector  or  is  scattered  into  the  receiver's  field 
of  view. 

There  are  three  sources  of  background  solar  radiation  that  may 
be  detected:  (1)  strongly  forward  scattered  solar  radiation  when  the 
sun  Is  within  the  viewing  angle  of  the  detector,  (2)  solar  radiation 
multiply  scattered  into  the  receiver's  viewing  (look)  angle,  and  (3) 
reflected  radiation  from  the  target  or  target  background  within  the 
receiver  field  of  view.  The  first  of  the  three  methods  is  the  largest 
background  radiation  source  and  a  most  important  laser  receiver  problem. 
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Fig.  k.  Sun-Target-Receiver  Geometry  where  6g  is  Sun's 
Zenith  Angle,  9  is  Receiver’s  Zenith  Angle, 

<j>  is  Receiver's  Azimuthal  Angle 
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In  general,  our  system  receiver  will  be  viewing  a  target  located  on 
the  ground.  If  the  receiver  is  at  a  high  zenith  angle  and  the  sun  is 
likewise  at  a  high  zenith  angle,  then  the  receiver  views  only  solar 
radiation  that  has  been  scattered  into  its  field  of  view  or  reflected 
from  the  target  area.  Although  the  radiation  from  the  sun  is  most 
intense  when  the  sun  is  near  zenith  due  to  the  shorter  path  through 
the  atmosphere,  the  amount  received  is  greatly  reduced  by  the  receiver's 
viewing  angle  and  direction  of  viewing.  A  significant  problem  exists, 
however,  in  the  event  that  the  sun  and  receiver  zenith  angles  approach 
a  maximum  and  the  azimuthal  angle  of  the  receiver  is  such  that  the 
receiver  opposes  (faces)  the  sun.  Even  though  the  amount  of  solar 
radiation  is  a  minimum  due  to  its  longer  path  through  the  atmosphere, 
this  solar  radiance  with  a  strong  forward  scatter  component  is  incident 
approximately  normally  on  the  detector  surface. 

According  to  Reference  61,  the  irradiance  of  the  sun  at  a 
zenith  angle  of  78.5°  is  b  watts/m2  at  sea  level.  Many  factors  affect 
the  magnitude  of  this  irradiance  such  as  cloud  cover,  temperature, 
precipi table  water  vapor  in  the  air,  pressure,  number  of  particles  in 
the  air,  but  this  reported  value  is  representative.  Such  irradiance, 
which  is  nearly  normally  incident  on  the  receiver,  results  in  a  serious 
system  limitation  especially  if  the  receiver  detector  is  otherwise 
capable  of  detecting  on  the  order  of  10"6  to  10“9  watts/m2.  This 


problem  is  complicated  by  the  fact  that  receivers  have  a  wide  field 
of  view  and  therefore  subtend  a  significant  portion  of  this  radiation 
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Reduction  of  the  field  of  view  is  a  method  of  alleviating  this  problem, 
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but  this  method  exacts  a  large  cost  in  target  acquisition  probability. 

Elimination  of  this  system  deficiency  may  be  accomplished  by 
restricting  the  tactical  use  of  the  weapon  system  such  that  during 
periods  of  low  sun  angles  the  receiver  approaches  the  target  in  such 
a  manner  as  to  subtend  the  smallest  receiver  view  angle  intersection 
with  the  forward  scatter  component  of  the  solar  irradiance.  This 
method  means  the  receiver  would  be  restricted  to  use  only  in  the  quad¬ 
rants  nearest  the  sun  (receiver's  back  to  the  sun)  and  to  high  angle 
deliveries— a  condition  unacceptable  In  a  hostile  tactical  environment. 
Another  method  of  eliminating  the  solar  background  radiance  is  designing 
the  receiver  logic  such  that  it  discriminates  against  continuous  returns 
as  opposed  to  pulsed  receptions.  The  design  of  system  logic  to  eliminate 
sensitivity  to  extended  returns  would  also  eliminate  the  problem  of 
detection  of  solar  reflected  radiation  from  the  target  and  target 
background  since  this  radiation  is  detected  as  a  steady  state  back¬ 
ground  component. 

Other  possible  sources  of  1 ,059u  to  1 ,069y  radiation  are  the 
moon,  stars,  and  hot  surfaces  within  the  field  of  view  of  the  detector. 
The  spectral  irradiance  of  the  night  sky  in  this  region  is  on  the 
order  of  10-12  less  than  the  day  sky  equivalent,  although  this  figure 
Is  affected  by  the  phases  of  the  moon  (Ref  2h:6-S).  The  laser 
designator  system  in  the  absence  of  scattered  solar  radiation,  then 
may  detect  nanowatts  of  laser  designator  energy  even  without  internal 
system  logic.  On  the  other  hand,  any  object  within  the  field  of  view 
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It  were  extremely  hot.  Spectral  curves  for  blackbodies  appear  in 
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References  2k  and  33  which  indicate  that  objects  with  temperatures  in 
excess  of  500°K  emit  significant  radiation  in  the  1.059P  to  1.069u 
region.  This  too  may  be  eliminated  by  the  design  logic;  but  in  the 
event  no  logic  utilized,  this  might  be  a  system  limitation. 

For  the  duration  of  this  report,  the  problem  of  background 
signal  will  be  assumed  to  have  been  eliminated  through  the  design  logic 
system  which  discriminates  the  extended  return  from  the  pulsed  return. 
If  this  assumption  were  not  possible  then  the  more  difficult  problem 
of  determining  the  background  signal  for  various  receiver  and  sun 
zenith  angles  must  be  addressed  at  length.  This  problem  is  more  com¬ 
monly  known  as  the  radiative  transfer  problem.  This  is  a  much  more 
complex  problem  and  cne  that  has  only  recently  been  addressed  for  the 
visible  region  (Ref  19). 

In  the  absence  of  a  background  signal  reaching  the  receiver, 
the  Johnson  noise  term  (2FKTB)  and  the  dark  current  term  (2eBljReG2) 
are  the  only  remaining  terms  in  the  denominator  of  Eq  (7).  The  design 
specifications  of  the  receiver  then  uniquely  determine  the  minimum 


power,  Pmr,  that  can  be  detected  by  the  receiver  for  a  given  signal- 

to-root-mean-square  noise  ratio.  This  value  of  P _  is  the  minimum 

detectable  laser  signal.  When  this  power  is  equal  to  the  power  at  the 
receiver,  Pr,  computed  in  Eq  (6),  then  this  determines  the  maximum 
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which  is  the  atmosphere-free  lock-on  range  equation. 

Unfortunately  the  signal-to-noise  ratio  is  not  constant  but 
varies  with  temperature  for  the  non-cooled  detector  heads  currently 
In  use.  A  realistic  signal-to-noise  ratio  (SNR)  may  be  computed, 
however,  based  on  typical  mission  profiles  and  temperatures. 

It  is  not  the  intent  of  this  report  to  explore  the  character¬ 
istics  of  various  receivers  and  lasers.  A  typical  laser  peak  power 
and  designator  optics  transmission  value  are  assumed.  The  target  is 
assumed  to  be  a  Lambert  emitter  with  reflectivity  of  10%.  Further,  if 
typically  achievable  values  are  assumed  for  the  receiver  transmittance, 
Tr,  and  are*  of  the  receiver  optics,  Ar,  then  the  right  hand  side  of 
Eq  (8)  is  reduced  to  a  constant  if  one  assumes  that  the  receiver  views 
the  target  normally  (cos  0t  =  l) ,  an  optimum  condition  .  Using  the 
de.-ign  parameters  from  Reference  2k,  then  the  right  hand  side  of  Eq 
(8)  is  reduced  .o  6.7  x  109m2  (Ref  2^4 : •  It  is  important  to 
recognize  that  this  value  is  for  one  particular  set  of  design  parameters; 
another  set  of  parameters  will  yield  a  different  value,  but  the  method 
of  solution  is  the  same  and  Eq  (8)  is  a  general  result  (for  the 
atmosphere-free  lock-on  range). 

in  summary,  we  have  discussed  the  laser,  the  target  reflective 
characteristics,  and  the  receiver  characteristics  and  their  relationships 
to  the  maximum  lock-on  range.  In  the  development  of  the  range  equation, 
we  have  assumed:  (1)  there  is  no  loss  of  laser  power  to  the  atmosphere, 
(2)  the  beam  at  the  target  is  smaller  than  the  target  and  is  reflected 
only  from  it  and  not  the  background,  (3)  the  target  is  a  perfectly 
diffuse  Lambertian  ref lector  wi th  reflectivity  of  10%,  {k)  the  receiver 
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IS  so  designed  so  as  to  discriminate  against  continuous  background 

returns,  and  (5)  that  the  target  is  entirely  within  the  field  of  view 
of  the  receiver. 
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Mi.  Atmospheric  Attenuation  of  Laser  Radiation 

As  the  maximum  lock-on  range  equation  was  developed  in  Chapter 
if,  two  important  laser  energy  losses  were  ignored — both  of  them  the 
result  of  atmospheric  attenuation.  As  the  radiant  power  travels 
through  the  atmosphere,  some  of  it  is  scattered  out  of  the  path  to  the 
target  or  receiver  and  some  of  the  radiation  is  absorbed.  The  scat- 
terers  and  absorbers  are  the  atmospheric  gaseous  molecules,  aerosols, 
and  water  droplets.  Unless  the  designator  and  receiver  are  co-lccated, 
then  the  two  paths  through  the  atmosphere  are  different;  and  therefore, 
the  energy  losses  due  to  the  atmospheric  attenuation  are  not  the  same. 
Moreover,  the  nature  of  the  radiation  transfer  in  the  two  paths  is  not 
the  same  due  to  the  assumption  that  the  target  is  a  Lambertian  surface. 
The  two  paths  will  be  addressed  separately.  The  two  losses  wll 1  be 
described  mathematically  and  then  related  to  the  maximum  lock-on  range 
equation  from  Chapter  U. 

In  describing  the  atmospheric  energy  loss,  this  report  uses 
the  following  terminology.  Atmospheric  attenuation  is  the  process  by 
which  energy  is  scattered  and  absorbed  by  the  atmospheric  constituents. 
Attenuation  coefficients  are  the  constants  used  to  describe  the 
attenuation  due  to  aerosols,  molecules,  rain  and  others.  Extinction 
coefficient  refers  to  the  sum  total  of  the  separate  attenuation 
coefficients.  That  is,  extinction  refers  to  the  sum  total  of  all 
scattering  and  absorption  (Ref  62:5),  although  this  atmospheric  process 
is  called  attenuation  in  this  report,  in  much  of  the  literature 
concerning  this  subject,  attenuation  and  extinction  are  used  inter¬ 
changeably. 
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Atmospheric  Effects  on  Laser  Radiation  (Laser  to  Target) 

Assuming  that  there  is  some  constant,  at,  by  which  the 
atmospheric  absorption  and  scattering  may  be  described  for  given  atmos¬ 
pheric  conditions,  then  the  intensity  of  radiation  at  any  point  in  space 
is  given  by  the  fundamental  law  of  extinction  of  Beer-Lambert  (or 
Bcurguer,  in  European  literature),  it  states  that  the  extinction 
process  is  linearly  independent  with  respect  to  the  intensity  of 
radiation  and  amount  of  matter,  provided  that  the  physical  state  (i.e., 
temperature,  pressure,  composition)  is  held  constant  (Ref  33:23).  This 
lav/  postulates  that  matter  in  the  optical  path  is  in  the  same  physical 
state  and  therefore  extinction  is  dependent  only  on  the  amount  of  mat¬ 
ter.  Then  the  monochromatic  spectral  transmission  of  radiation  in  the 
atmosphere  is  (Ref  ^9 : I C) 


-/L  cr  (£}d£ 
T  =*  e  °  t 


(9) 


where  ot  =  total  atmospheric  extinction  coefficient  per  unit  path 
length, 

L  =  total  path  length, 
and  dZ  -  incremental  path  length. 

Van  de  Hulst  in  Reference  62  states  that  for  diffuse  light,  Lambert’s 
law  is  valid  for  single  scattering  only  for  an  optical  depth  of  0.1 
or  less  where  one  optical  depth  is  defined  as  the  exponent  required 
to  reduce  the  transmission  to  I/e  of  its  original  value.  For  optical 
depths  of  between  0.1  and  0.3,  a  correction  term  is  necessary  to 
account  for  double  scattering  arid  for  optical  depths  in  excess  of  0.3 


Ik 


GEP/PH/74-4 


then  the  full  complexities  of  multiple  scattering  become  a  factor 
(Ref  62:6).  Van  de  Hulst's  work,  however,  refers  only  to  diffuse 
sources  of  light  and  not  to  the  narrow  collomated  laser  light.  Zuev, 
Kalanov,  and  Savel 'eo  in  Reference  73  have  investigated  the  applica¬ 
bility  of  the  Beer-Lambert  law  to  laser  propagation.  They  experi¬ 
mentally  determined  the  validity  of  the  exponential  law  for  laser 
radiation  up  to  optical  depths  of  12  (Ref  73:1^0).  For  optical  depths 
greater  than  9,  however,  multiple  scattering  may  be  significant  for 
targets  that  are  larger  than  9  times  the  diameter  of  the  laser  beam. 
Optical  depths  larger  than  12  were  not  experimentally  verified  by  Zuev 
because  of  light  source  and  receiver  limitations,  but  this  reference 
predicts  valid  ty  of  the  law  to  optical  depths  of  30. 

In  the  problem  of  the  laser  radiation  incident  on  the  target, 
the  loss  of  energy  can  be  determined  then  by  use  of  the  Beer-Lambert 
exponential  law.  The  power  incident  on  the  surface  of  the  target  is 
then 

-/Rd  crtU)dJl 

P..  =  P  T.e  °  (10) 

t  p  d 

where  at(£)  is  the  atmospheric  extinction  coefficient  per  unit  path 
length  as  a  function  of  position,  d£  is  the  incremental  path  length, 
and  Rj  is  the  distance  from  designator  to  target.  Thus,  the  radiant 
power  on  the  target  is  a  function  of  the  specific  atmospheric 
extinction  coefficient,  crt. 

Atmospheric  Effects  on  Reflected  Laser  Radiation  (Target  to 
Receiver).  On  the  surface  it  would  appear  that  there  need  be  no 
Justification  of  the  use  of  the  exponential  law  for  the  reflected  laser 
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radfatlon.  Since  our  model  has  assumed  a  diffuse  Lambert  reflector, 
however,  justification  is  essential.  Conversely,  if  our  model  had 
assumed  specular  reflection  by  the  target,  then  the  properties  of  the 
laser  radiation  would  have  been  preserved.  Such  an  assumption  is  not, 
however,  in  agreement  with  the  reflective  characteristics  of  most  mili¬ 
tary  targets  as  we  have  previously  discussed;  in  addition,  such  an 
assumption  makes  determination  of  directional  reflectance  exceedingly 
complex. 

The  diffuse  reflection  of  laser  from  the  target  requires  that 
for  single  scattering  the  Lambert  law  is  limited  to  optical  depths  of 
0.1  or  less  as  described  by  Van  de  Hulst.  For  larger  optical  depths, 
the  full  multiple  scattering  problem  must  be  addressed.  For  extremely 
narrow  receiver  fields  of  view,  there  is  justification  for  use  of  the 
Lambert  law  (Ref  71:202),  but  for  use  in  our  model  it  must  be  acknow¬ 
ledged  that  the  use  of  the  Lambert  law  for  single  scattering  for  the 
atmospheric  attenuation  of  radiation  from  target  to  receiver  is  a  mode! 
limitation.  Because  the  multiple  scattering  problem  is  not  addressed 
for  this  portion  of  the  energy  transport,  then  this  model  may  result 
in  underestimation  of  the  maximum  lock-on  range. 

With  this  caution,  hen  the  assumption  is  made  that  energy 
transport  from  the  target  \  iiver  is  approximated  by  the  Lambert 
Law  and  the  maximum  laser  lock-on  range  equation  becomes 
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where  0^.(2.)  Is  the  atmospheric  extinction  coefficient  per  unit  path 
length  as  a  function  of  position,  di  is  the  incremental  path  length, 
and  Rr  is  the  distance  from  target  to  receiver.  Substituting  the 
design  parameters  constant  (discussed  in  Chapter  II)  into  Eq  (11)  then 
the  laser  target  designator  system  maximum  lock-on  range  equation  be¬ 
comes 
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-  (6.7  x  109m2)e 
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(12) 


which  is  a  transcendental  equation  that  is  dependent  on  the  specific 
atmospheric  extinction  coefficients. 


Atmospheric  Extinction  Coefficient 

The  extinction  coefficient,  as  we  have  discussed,  refers  to 
the  absorption  and  scattering  of  the  laser  radiation  by  atmospheric 
gases  and  particles.  That  is,  the  extinction  coefficient  is  the  sum 
of  all  absorption  and  scattering  coefficients.  It  is  the  sum  of  the 
coefficients  describing  molecular  scattering,  molecular  absorption, 
aerosol  scattering,  aerosol  absorption,  and  the  scattering  and  ab¬ 
sorption  by  clouds,  rain,  ice  crystals,  fog,  snow  and  others.  The 
total  atmospheric  attenuation  is  a  complex  function  of  the  number  of 
molecules,  kinds  of  molecules,  number  of  aerosols,  aerosol  composition, 
aerosol  size  distributions,  number  of  water,  ice,  fog,  snow,  or  cloud 
droplets,  droplet  composition  and  droplet  size  disti  fbution  in  the 
volume  of  air  traversed  by  the  radiation.  Moreover,  these  factors  are 
not  constant  throughout  the  atmosphere,  and  therefore  atmospheric 
attenuation  becomes  a  function  of  the  particular  radiation  path.  In 
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addition,  the  scattering  and  absorption  is  dependent  on  the  particular 
wavelength  of  the  radiation. 

The  problem  of  determination  of  atmospheric  attenuation  of 
1.06y  laser  radiation  may  then  be  reduced  to  the  following  subproblems: 
(l)  determination  of  molecular  scattering  and  absorption,  (2)  deter¬ 
mination  of  the  aerosol  scattering  and  absorption,  and  (3)  determination 
of  the  scattering  and  absorption  of  other  atmospheric  particles  (such  as 
clouds,  rain  and  fog),  and  the  relationship  of  these  to  observables  or 
measurable  quantities.  These  separate  forms  cf  scattering  and  absorption 
are  expressed  in  terms  of  separate  attenuation  coefficients.  In  general, 
these  separate  coefficients  are  additive  yielding  the  total  extinction 
coefficient. 

On  a  clear  day,  for  example,  the  atmospheric  transmission  of 
laser  radiation  is  a  function  only  of  the  aerosols  and  molecules 


where  ag  is  the  aerosol  attenuation  coefficient  (per  unit  length), 
om  is  the  molecular  attenuation  coefficient  (per  unit  length' 
and  L  is  the  path  length.  The  aerosol  attenuation  coefficien'  • 
the  sum  of  the  scattering  and  absorption 


o  -  a  +  8 
a  3  a 


(14) 


where  a  is  the  aerosoi  absorption  coefficient  (per  unit  length),  and 
8_  is  the  aerosol  scattering  coefficient  (per  unit  length).  Likewise, 
the  molecular  attenuation  coefficient  is  the  sum  of  the  molecular 
scattering  and  absorption 
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where  is  the  molecular  absorption  coefficient  (per  unit  length), 
and  is  the  molecular  scattering  coefficient  (per  unit  length). 

The  determination  of  these  individual  coefficients  is  the  subject 
matter  of  the  next  chapters.  We  shall  successively  investigate  each 
of  these  coefficients  by  determining  in  each  case  the  factors  that 
contribute  to  the  radiation  attenuation.  Next,  since  the  atmosphere 
Is  inhomogeneous,  the  factors  must  be  related  to  position  within  the 
atmosphere.  Last,  these  factors  must  be  related  to  atmospheric  obser¬ 
vables  or  measurables,  so  that  the  laser's  radiant  power  at  a 
particular  point  in  the  atmosphere  is  predictable. 
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IV.  Molecular  Absorption  and  Scattering 


Radiation  Is  scattered  by  all  of  the  atmospheric  gaseous  mole¬ 
cules,  but  It  is  selectively  absorbed.  There  are  many  atmospheric 
constituents.  All  are  important  for  the  determination  of  molecular 
scattering;  but  only  some  of  the  molecules  are  effective  absorbers  of 
radiation  and  these  are  effective  only  at  certain  wavelengths  which 
depend  upon  the  structure  of  the  molecule. 


Molecular  (Rayleigh)  Scattering 

Ail  atmospheric  molecules  scatter  radiation.  Their  size  is 
on  the  order  of  10“6p  to  1 0~8p  and  therefore  the  scattering  cross- 
section  is  quite  small.  Particles  which  are  small  relative  to  the 
wavelength  of  the  incident  radiation  experience  the  same  electric 
field  throughout  their  entire  dimensions,  which  are  assumed  to  be 
sphevical ,  This  external  electric  field  establishes  a  dipole  in  the 
small  particle  which  is  a  function  of  its  polarizability.  This  dipole 
then  emits  radiation  in  the  characteristic  dipole  pattern  which  results 
In  the  removal  of  energy  from  the  incident  beam  (Ref  8:184).  This 
phenomenon  is  called  Rayleigh  scattering,  and  the  scattering  coefficient 
can  be  expressed  as 


8ir3  (Mg2-1)2  6+3  f 

8  RAY  ~  6-7f 


(16) 


where  X  is  the  wavelength  of  incident  radiation  (cm),  Mg  is  the  re¬ 
fractive  index  of  the  gas  relative  to  a  vacuum,  N_  is  the  number  of 
molecules  per  cm3,  and  f  is  the  depolarization  factor,  which  for  air 
is  0.042  and  is  often  omitted  (Ref  8:189). 
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It  must  be  emphasized  that  this  equation  is  valid  only  for  spherical 
molecules  and  a  correction  term  must  be  added  for  non-spherical 
molecules.  This  equation,  however,  will  prove  useful  in  yielding  an 
estimate  of  the  molecular  scattering  in  the  atmosphere  provided  the 
number  of  molecules  can  be  determined.  According  to  Avogadro's 
hypothesis,  the  number  of  molecules  of  any  gas  per  unit  volume  depends 
only  on  temperature  and  pressure.  For  the  U.  S.  Standard  Atmosphere 
(15°C  and  1013  mb  mercury)  the  number  of  molecules  in  one  cubic  centi¬ 
meter  is  approximately  2.5^7  x  1019.  Then  for  the  U.  S.  Standard 
temperature  and  pressure,  the  Rayleigh  scattering  coefficient  for 
1.064y  laser  radiation  is 


„  *  frl  c.(| .opo;.)2-!]? _  i  8.16  X  10-S/an 

m  3(1.064  x  I  <r*e ■)*  2-547  *  ’013 

cm3  (17) 

where  M  =  1.000273S  (Ref  27:2). 

Expressed  in  inverse  kilometers,  the  molecular  scattering  coefficient 


(*m  -  x  (105  cm/km)  (18) 

The  Rayleigh  scattering  coefficient  for  1.064y  laser  radiation  at  15°C 
and  1013  mb  mercury  is  8.16  x  10"’4  (Km"1) .  The  molecular  scattering 
coefficient  for  the  visible  (A  =  .55p)  at  15°C  and  1013  mb  is  1.14  x 
10~2 (Km"1) . 

In  addition  at  higher  altitudes  (lower  temperatures  and  pres¬ 
sures)  this  coefficient  becomes  smaller.  For  slant  path  transmissions 
from  altitude,  then,  the  average  molecular  scattering  coefficient  would 


- a 
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be  even  less  significant.  Increasing  the  temperature  to  35°C  and 
pressure  to  1025  mb,  conditions  not  normally  encountered,  does  not 
significantly  affect  the  scattering.  From  these  calculations,  it 
becomes  clear  that  molecular  scattering  is  not  a  serious  limiting 
factor  for  reasonable  path  lengths  of  the  laser  radiation  and  may  be 
safely  ignored  for  1.06%  radiation.  Conversely,  the  Rayleigh  scat¬ 
tering  in  the  visible  may  not  always  be  ignored,  but  for  our  purposes 
now  it  is  rather  insignificant. 

Molecular  Absorption 

The  molecules  are  frequently  excellent  absorbers  of  radiation. 
Molecular  constituents  of  the  atmosphere  that  are  significant  absorbers 
are:  (1)  water  vapor;  (2)  carbon  dioxides;  (3)  ozone;  (4)  nitrous 
oxide;  (5)  carbon  monoxide;  (6)  methane;  and  (7)  oxygen.  Other  trace 
gases  such  as  sulfur  dioxide  may  also  be  good  absorbers  but  they  do 
not  naturally  occur  in  the  atmosphere  and  significant  concentrations 
are  limited  to  areas  with  large  industrial  activity. 

The  absorption  by  molecules  is  due  to  the  changes  in  molecular 
states  caused  by  the  addition  of  energy.  The  molecules  absorb  the 
energy  by  changing  vibrational,  rotational,  or  electronic  energy  levels 
or  a  combination  of  the  three.  Each  single  allowed  transition  gives  a 
particular  absorption  line.  For  each  allowed  vibrational  energy  level 
transition  there  are  many  rotational  energy  levels  possible  thereby 
creating  a  band  structure  or  rather  a  series  of  closely  spaced  lines 
which  overlap. 

A  great  amount  of  work  has  been  done  in  the  field  of  molecular 
absorption  and  molecular  spectroscopy.  For  a  detailed  account  of  the 
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theory  of  molecular  absorption,  the  reader  is  referred  to  References 
33,  35  or  36. 

From  the  exhaustive  studies  on  absorption,  there  have  emerged 
various  band  models  for  prediction  of  average  transmission  for  dif¬ 
ferent  spectral  regions.  Some  of  the  better  known  broadband  models 
are  the  Altshuler,  Statistical  or  Mayer-Goody,  Elssaser,  Random 
Elssaser,  and  the  Air  Force  Cambridge  Research  Laboratories  Lowtran 
li.  There  are  also  two  approximation  schemes  used  in  the  band  models: 
the  weak-line  and  the  strong-line  approximations. 

In  general,  these  band  models  are  based  on  experimental 
measurements  of  the  solar  spectral  transmittance  or  laboratory  data. 
The  data  is  then  smoothed  to  provide  continuous  curves  and  average 
transmittance  in  the  various  spectral  regions.  In  addition,  these 
models  are  statistically  extrapolated  into  spectral  regions  where  data 
is  not  available.  For  purposes  of  predicting  low  resolution  spectral 
transmittance,  a  band  model  is  adequate  if  it  is  based  on  measurements 
in  the  spectral  region.  However,  for  laser  transmittance  where  high 
spectral  resolution  is  required,  then  the  band  models  are  woefully 
Inaccurate.  Laser  wavelengths  are  assumed  to  be  monochromatic  so  an 
average  transmittance  for  a  spectral  region  does  not  provide  laser 
transmittance.  To  determine  laser  transmittance,  a  line  by  line 
calculation  is  necessary. 

Such  a  compilation  of  molecular  absorption  lines  has  been  the 


ongoing  project  of  the  Air  Force  Cambridge  Research  Laboratories  since 
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vibration-rotation  lines  of  seven  naturally  occurring  molecules  of 
significance  (including  their  isotopes)  in  the  terrestrial  atmosphere. 
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These  seven  are  the  same  atmospheric  constituents  listed  earlier  in 
this  section.  With  the  exception  of  oxygen,  all  of  these  molecules 
are  minor  atmospheric  constituents,  but  they  represent  most  of  the 
absorption  lines  in  the  visible  and  infrared.  This  work  assumes 
(and  investigation  seems  to  support  the  fact)  that  all  of  these  gases, 
with  the  exception  of  water  vapor  and  ozone,  are  uniformly  mixed  in 
the  atmosphere.  Each  of  the  lines  in  this  study  (Ref  47)  represents 
either  a  theoretical  calculation  based  on  atomic  and  molecular  theory 
or  actual  solar  or  laboratory  data  or  both. 

In  order  to  compute  the  transmittance  of  a  given  spectral  line 
in  the  atmosphere  it  is  necessary  to  describe  the  absorption  coef¬ 
ficient  as  a  function  of  frequency  for  each  line.  The  four  essential 
parameters  for  each  line  are  the  resonant  frequency,  v0(cm-1),  the 
intensity  per  absorbing  molecule,  S(ctrf  Vmolecule  cm-2),  the  Lorentz 
line  width  parameter,  cx0(cm"1/atm) ,  and  the  energy  of  the  lower  state, 
E'^crn-1).  The  frequency,  vQ,  is  to  a  close  approximation  independent 
of  temperature  and  pressure.  The  intensity,  S,  is  pressure-independent 
and  its  temperature  dependence  can  be  calculated  from  E”  and  v.  The 
line  half-width  at  half  maximum,  a,  is  proportional  to  the  pressure 
and  is  temperature  dependent.  The  precise  line  shape  is  uncertain 
but  theoretical  computations  usually  begin  with  the  Lorentz  line  shape, 
which  has  two  limitations:  (l)  Lorentz  lineshape  for  infrared 


frequencies  requires  modification  for  pressures  lower  than  100  mb 
(which  is  not  a  factor  in  our  case), and  (2)  The  Lorentz  lineshape  may 


be  inaccurate  in  the  distant  wings  of  a  line  \iv— or  when  col¬ 


lision  broadening  forces  are  dipole-quadrupole  (Ref  47:2,3). 
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Reference  47  gives  the  intensity  as  a  function  of  a  temperature 


as 


S(T)  = 


S(Ts)  Qv(Ts)Qr{Ts) 
Qv(T)  Qr(T) 


1.439E"(T-Ts) 


TT, 


(19) 


Where  E"  is  the  energy  of  the  lower  state  of  transition  (cm-1),  Qv 
is  the  vibrational  partition  function,  Qr  is  the  rotational  partition 
function,  T  is  the  temperature  (°K),  and  Ts  is  296°K.  Assuming  a 
Lorentzian  line  shape,  the  line  half  width  at  half  maximum  at  any 
temperature,  T,  and  pressure,  P,  is  given  by 


a  =  a 


o  P. 


(20) 


where  aQ  is  the  measured  or  computed  line  half  width  at  temperature 
Tq  and  pressure  PQ;  PQ  is  the  pressure  at  which  experimental 
measurement  or  theoretical  determination  of  aQ  is  made,  and  TQ  is  the 
temperature  of  the  experimental  measurement.  The  molecular  absorption 
coefficient  due  to  the  absorption  by  one  molecule  of  one  molecular 
species  for  one  absorber  line  is 


e 


S  a 


lm 


tt[(v-v0)2  +  a2  3 


(21) 


where  S  is  the  line  intensity, 
a  is  the  line  half-width, 
vQ  is  the  central  line  frequency, 

and  v  is  the  laser  frequency,  which  is  assumed  to  be  monochromatic. 
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Generally,  however,  a  large  number  of  absorption  lines  belonging  to 
different  molecular  species  contribute  to  the  attenuation  at  a  spe~ 
cific  laser  frequency;  that  is,  the  absorption  coefficient  per 
molecule  for  a  laser  of  frequency,  v,  Is 


Sij  “ij 


aWm  j  5  nE(v-Vjj)2"+  a.j^] 


(22) 


where  subscript  i  refers  to  the  different  spectral  lines  of  each 
molecular  gas  and  subscript  j  refers  to  the  j  different  molecular 
gases.  Multiplying  the  absorption  coefficient  by  the  total  number 
of  molecules  along  with  the  path  length  yields  the  optical  depth. 
Another  more  conventional  method  for  our  purposes  involves  defining 
the  total  absorber  density  along  the  path.  Multiplying  am/n)  by  this 
absorber  density  gives  the  molecular  absorption  coefficient  per  unit 
path  length,  which  is  the  way  we  have  defined  all  absorption  and 
scattering  coefficients.  Then 


S* .  ct . .  m.  # ^ ^ . 

a  .  ZZ  'J  ij  j _  (23) 

m  ttC  (v-v.  j)2  +  ajj23 

where  mj  is  the  number  density  of  each  molecular  gas,  j,  per  unit 
path  length. 

Distribution  of  Molecular  Absorbers.  If  the  density  of  each 
molecular  absorber  were  known  for  all  points  P  in  the  atmosphere  then 
Eq  (23)  would  yield  the  molecular  absorption  coefficient.  The  mole¬ 
cular  concentrations  of  the  atmospheric  gases,  however,  are  not 
uniform.  They  vary  for  the  different  temperatures,  pressures, 
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locations,  and  air  masses.  Several  approximations  will  be  made. 

The  first  assumption  is  that  the  atmosphere  is  horizontally 
stratified.  All  horizontal  paths  are  defined  as  being  through  uni¬ 
form  air  masses.  By  this  stratification,  two  categories  of 
atmospheric  gases  emerge:  (1)  the  uniformly  mixed  gases,  and  (2) 
the  unmixed  gases.  Of  the  seven  molecular  absorbers  under  consideration 
only  two  fail  to  qualify  for  the  first  category.  These  two,  ozone  and 
water  vapor,  are  not  only  a  function  of  atmospheric  conditions  but  they 
are  also  a  function  of  geographical  location  and  their  evolution 
properties. 

For  the  uniformly  mixed  gases,  the  concentration  is  generally 
written  as  the  equivalent  length  of  pure  absorber  in  the  path  as  a 
function  of  the  pressure  and  temperature.  For  the  horizontally 
stratified  atmosphere,  this  concentration,  AL,  is  the  product  of  the 
path  length,  R(cm),  the  fractional  part  of  the  total  atmospheric 
volume  of  absorber,  f,  and  the  pressure  ?n  atmospheres,  P,  This 
concentration  is  most  conveniently  expressed  in  units  of  cm-atm 
reduced  to  standard  temperature  and  pressure  (STP)  (P.ef  49:87). 

The  AFCRL  report,  Optical  Properties  of  the  Atmosphere ,  gives 
this  molecular  concentration  for  a  horizontal  path  at  any  altitude, 
z,  as 


Bfe)  (cm~atm)gTp  m 

where  Ts  **  273°K  and  Ps  =  1  atm.  This  concept  of  equivalent  absorber 
path  at  standard  conditions  may  be  extended  to  vertical  paths  provided 
the  vertical  molecular  distribution  may  be  described.  For  uniformly 
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mixed  gases  of  fractional  concentration  by  volume,  f,  the  vertical 
distribution  is  modeled  using  the  hydrostatic  equation  which  states 
that  the  change  In  pressure  along  a  smaP  vertical  path  is  the  pro- 
duct  of  the  molecular  density,  gravity,  and  path  length.  If  the 
hydrostatic  equation  is  integrated  from  an  altitude,  z,  out  to 
infinity,  then  the  hydrostatic  equation  may  be  expressed  In  terms  of 
the  scale  height  of  the  atmosphere  above  altitude,  z,  H(z) 


P(z)  -  gp(z)H(z)  (25) 

where  P(z)  =  pressure  at  altitude  z, 
p(z)  =  molecular  density  at  z, 
g  «  acceleration  due  to  gravity 

For  a  vertical  column  of  air  at  STP,  H(o)  ®  7-93  Km  (Ref  49:88). 

The  scale  height  in  Km  for  any  altitude  z  may  be  determined  by 

H(z)  =  10.2  (26) 

{*(*>/ 


After  some  manipulation  then  the  amount  of  absorber  for  a  slant  path 

between  altitudes  z  and  z  becomes 

1  2 


f  x  10. .2 


P(z,)-P(z_) 


Sec  3 


where  8  is  the  angle  from  the  vertical  (zenith).  Values  for  the  con¬ 
centration  of  the  uniformly  mixed  gases  along  with  values  in  (cm-otm) 

STP 

for*  Vertical  paths  f rOiii  sea  level  are  giyen  in  Table  1.  These  Values 


for  vertical  paths  when  multiplied  by  the  fractional  pressure  change 
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and  the  secant  of  the  zenith  angle  yield  the  (cm-atm)  of  absorber 

STP 

between  two  pressure  levels  or  rather  two  altitudes. 


Table  i 

Concentrations  of  Uniformly  Mixed  Gases 


Constituent 

Molecular 

Weight 

PPM  by 

Volume 

(cm-atmL.p  in 

Vertical  Path 
From  Sea  Level 

Air 

28.97 

106 

8  x  105 

CO 

2 

44 

330 

264 

NO 

2 

44 

0.28 

0.22 

CO 

28 

0.075 

0.06 

CH 

4 

16 

1.6 

1.28 

0 

2 

32 

2.095  x  105 

1.68  x  105 

This  result  is  useful  in  that  it  allows  us  to  determine  the 


number  of  molecules  per  unit  area  in  the  path.  Since  there  are 


always  Avogadro's  number  of  molecules  in  22 .4  liters  at  STP  then  the 
following  relation  for  the  uniformly  mixed  gases  may  be  derived 


1  (cm-atm)  =  2.69  x  IQ19  molecules/cm2  (28) 

STP 


Unfortunately,  the  gases  ozone  and  water  vapor  do  not  behave 
so  simply.  These  gases  are  more  dependent  on  their  particular 
production  processes  and  their  concentrations  in  the  atmosphere  3re 


utr/rn/  /‘*-h 


subject  to  wide  variation.  Ozone  production  and  distribution  is 
important  for  some  absorption  problems,  but  for  purposes  of  this 
report,  its  contribution  for  molecular  absorption  is  negligible 
for  two  reasons:  (l)  ozone  concentrations  at  or  below  20,000  feet 
are  quite  small,  and  (2)  absorption  lines  of  ozone  of  sufficient 
intensity  are  not  within  the  spectral  region  of  interest.  For  these 
reasons,  this  discussion  will  be  restricted  to  the  atmospheric  con¬ 
centrations  of  water  vapor. 

The  amount  of  water  vapor  in  air  depends  on  such  diverse 
atmospheric  conditions  as  the  wind  speed,  atmospheric  pressure  and 
temperature  and  on  the  history  of  the  air  mass  and  on  geographical 
location.  And  the  concentration  of  water  vapor  does  not  necessarily 
decrease  with  altitude  or  the  formation  of  clouds  at  an  inversion 
would  be  a  rare  phenomena. 

The  amount  of  water  vapor  in  a  short  path  is  determinable 

provided  the  temperature  and  relative  humidity  are  known.  This 

quantity  of  water  vapor  is  given  in  total  length  of  water  (liquid) 

that  may-  be  precipitated  out  of  the  path  per  unit  area.  The 

prectpi table  water,  W,  may  be  defined  as  the  amount  of  water  contained 

in  a  vertical  column  of  air  of  unit  cross  section  extending  between 

two  pressure  levels  P  and  P  .  St  is  given  by 

12 

,  P 

W  =  —  /  2  wdP  (29) 

9  Pi 

where  g  ia  rhe  acceleration  of  gravity  and  w  is  the  mixing  ratio 
(Ref  6l:3“32).  If  w  is  assumed  constant  for  the  small  incremental 
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change  In  P  and  the  hydrostatic  equation  is  used  then  Eq  (29)  becomes 


> 


a 

"A 


W  «  wpwdz'«  •  (30) 

where  w  Is  the  mean  mixing  ratio  in  the  path  dz  and  where  the  density 
of  water,  Is  constant. 

The  mixing  ratio  w  is  defined  as  the  weight  of  water  vapor 

contained  in  mixture  with  a  unit  weight  of  dry  air  and  Is  related  to 

the  relative  humidity  (RH)  by  RH  =  100  x  w/w  where  w  is  the 

s  s 

saturation  mixing  ratio  (Ref  9:1 55) •  For  water  vapor  the  saturation 
mixing  ratio  is  given  approximately  by 


ws  %  ^  (31) 

where  e  =  .622,  a  ratio  of  the  molecular  weights  and  eg  is  the 
saturation  vapor  pressure  at  a  particular  temperature  (Ref  37:59)* 
The  saturation  vapor  pressure  at  a  temperature  t  in  °C  is  given 
approximately  by  the  empirical  relationship  of  Teten  (1930): 


[7  -  5t/ (237 . 3  +  0] 

e  %  6.11  x  10 


(32) 


where  eg  is  in  millibars  and  t  is  the  air  temperature  in  °C  (Ref  57:9). 
Combining  Equation  (30),  (31),  and  (32)  and  using  the  ideal  gas  law, 
the  following  relationship  is  derived  for  the  precipitable  water 
vapor: 


1 1  r  a.  /  /nm 

1/  .?w  \*0/  O 


A.W 

WJ 


W  %  (R.H.) (.458  x  10"6cm"3)  10 


dz 


(33) 
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Dividing  Eq  (33)  by  the  path  length  in  cm  and  converting  to  km  then 
the  precipitable  water  vapor  in  grams  per  cm2  per  km  is 

[7.5t/(237.3  +  t)] 

W  .0458  (10  )  (RH)  (34) 

If  this  incremental  path  is  taken  as  the  vertical,  then  to  determine 
W  for  slant  paths,  this  value  in  Eq  (34)  must  be  increased  by  a 
factor  of  secant  6  where  0  is  the  angle  from  the  vertical  (where 
6  <  80°) . 

The  water  vapor  concentration  per  km  path  length  as  a  function 
of  relative  humidity  and  temperature  is  shown  in  Figure  5.  Further, 
the  number  of  gm  cm~2/km  of  water  vapor  is  directly  convertible  into 
molecules/cm2  per  km  for  use  in  Eq  (23)  by  the  relation  derived  in 
Reference  49: 

1  gm/cm2-km  =  3*34  x  1022  molecules/cm2  (35) 

Distribution  Data 

At  this  point,  a  note  of  caution  is  in  order.  Neither  of  the 
expressions  derived  is  exact.  They  are  valid  only  for  small  path 
lengths  and  are  valid  only  when  the  pressure,  temperature  and  relative 
humidity  profiles  for  the  paths  are  known.  These  profiles  fluctuate 
greatly  in  just  a  few  hours  so  an  exact  determination  of  molecular 
absorption  is  essentially  impossible  even  within  the  bounds  of  these 
simplifications.  Very  few  observations  are  taken  of  the  vertical 
profiles  and  these  are  usually  rawinsonde  balloons  which  are  not 
absolutely  accurate. 
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The  rawinsonde  provides  little  data  on  the  first  200  feet 
above  ground  level  because  of  the  lag  in  lock-on  time.  Further,  the 
temperature  is  accurate  only  within  1°C  and  the  relative  humidity 
within  5  percent  (Ref  45). 

These  balloon  observations  are  relatively  sparce.  The  most 
that  are  taken  by  an  individual  station  is  four  in  a  24  hour  period. 
Most  stations  take  none  and  a  few  take  one  or  two.  For  example,  in 
an  area  the  size  of  Germany,  there  was  but  one  station  that  took  4 
observations  per  day  in  July,  1970  (Ref  51).  Nevertheless,  a  data 
Base  for  vertical  profiles  has  been  developed  and  utilized  by  weather 
services  worldwide  in  support  of  aviation.  Global  Weather  Central 
has  the  additional  capability  of  prediction  of  upper  level  winds, 
temperatures  and  pressures  for  all  hours. 

Part  of  this  data  base  has  been  used  by  AFCRL  in  its  compi¬ 
lation  of  model  atmospheres.  These  model  atmospheres — Tropical , 
Midlatitude  Summer,  Midlatitude  Winter,  Subarctic  Summer,  Subarctic 
Winter,  and  the  U.  S.  Standard— have  been  compiled  from  a  large 
number  of  observations  for  the  specific  regions  described.  These 
are  average  conditions  and  therefore  give  only  mean  conditions. 

These  atmospheres  reported  in  Ref  61  and  Ref  49  were  used  for  the 
theoretical  molecular  absorption  calculations  described  earlier  in 
this  paper. 

Comparison  of  Distributions.  The  vertical  profiles  described 
in  the  Midlatitude  Summer,  Midlatitude  Winter,  and  U.  S.  Standard 
models  are  shown  in  Table  II.  The  profiles  are  shown  from  the 
surface  to  the  500  mb  pressure  level  which  corresponds  to  an  altitude 
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of  approximately  1 7 » 000  feet  to  19,700  feet  depending  on  atmospheric 
conditions. 

Table  l! 

Vertical  Profiles  of  the  Midlatitude  Summer,  Midlatitude 
Winter,  and  1962  U.  S.  Standard  Atmospheres:  Sur¬ 
face  to  500  millibars  (mb) 


Relative  Humidity  (%)/Temperature  (°C) 


Season 

Model 

Surface 

850  mb 

700  mb 

500  mb 

■Winter 

Midlatitude 

80/-1 

75/-5 

60/-11 

50/-2S 

Winter 

(.35)a 

(.22) 

(.127) 

(.0315) 

Summer 

Midlatitude 

79/21 

60/14 

50/6 

32/— 1 1 

Summer 

(M) 

(.70) 

(.33) 

(.066) 

Fall/ 

U.  S. 

60/15 

67/6 

70/-6 

74/-24 

Spring 

Standard 

(.772) 

(.43) 

(.205) 

(.053) 

aValues  in  parentheses  are  the  Water  Vapor  Concentrations  in  gm-cm”2/ 
km  determined  for  the  Temperature  and  Relative  Humidity  using  Figure 
5. 

The  vertical  profiles  for  an  area  of  Germany  of  longitude 
10°E  and  latitude  from  48°  to  49°N  were  provided  by  USAF  ETAC 
(Environmental  Technical  Applications  Center)  that  represent  means 
of  observations  over  a  six  year  period.  This  data,  shown  in  Table 
iii,  was  provided  b"  month  and  then  averages  for  the  seasons  were 
computed.  The  water  vapor  concentrations  for  both  sets  of  data  were 
computed  using  Figure  5.  These  concentrations  are  shown  in 
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parentheses  in  the  tables. 


Table  II! 

Average  Monthly/Seasonal  Profiles  of  Relative  Humidity 
and  Temperature  for  Germany  (^9°  N  10  E) 
Surface  to  500mb* 


Relative 

Humidity(^)  / 

Temperature  (°C) 

Season 

Month 

Surface 

850mb 

700mb 

500mb 

Dec 

93/0 

69/-1 

51/-8 

39/-24 

Winter 

Jan 

80/2 

54/-2 

44/-9 

A2/-26 

Feb 

80/2 

63/-2 

52/-10 

47/-26 

Avg 

84/1 

621-2 

49/-9 

43/-25 

(.42)a 

(.28) 

(.12) 

(.020) 

Mar 

65/6 

69/0 

53/-8 

46/-25 

Spring 

Apr 

61/10 

69/2 

57/-7 

44/-23 

May 

67/14 

65/6 

49/-4 

37/-19 

Avg 

64/(0 

68/3 

53/-6 

42/-22 

(.58) 

(.40) 

(.17) 

(.038) 

Jun 

77/16 

66/10 

59/-1 

50/- 17 

Summer 

Jul 

77/19 

66/12 

59/1 

35/- 14 

Aug 

77/20 

57/12 

65/2 

42/-14 

Avg 

77/18 

63/11 

61/1 

42/- 15 

(1.2) 

(.63) 

■  (.31) 

(.065) 

Sep 

77/15 

61/10 

51/1 

46/- 15 

Fall 

Qct 

87/9 

61/6 

46/-2 

34/- 18 

Nov 

75/4 

69/1 

53/-6 

48/-22 

A.vg 

80/9 

64/6 

50/-2 

43/- 18 

(.66) 

(.45) 

(.20)’ 

(.053) 

*Data  provided  by  USAF  Environmental  Technical  Applications  Center 

^Values  in  parentheses  are  the  Water  Vapor  Concentrations  in 
gffl-crrT2/km  determined  for  the  Relative  Humidity  and  Temperature 
using  Figure  5. 
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Comparisons  of  the  concentrations  show  that  the  Midlatitude  Summer 
reasonably  approximates  the  mean  summer  concentrations  in  Germany. 
Likewise,  the  Midlatitude  Winter  model  approximates  the  winter  pro¬ 
files  for  Germany.  A  correlation  v/as  also  found  between  the  U.  S. 
Standard  model  and  Fall  and  Spring  conditions  reported  for  Germany, 
although  the  correlation  was  not  nearly  as  good  for  the  latter.  For 
the  Spring  Season,  the  model  tends  to  be  overly  pessimistic  for  mean 
atmosphere  conditions  for  this  area  of  Germany.  Nevertheless,  the 
correlation  of  the  three  models  with  the  four  seasons  in  this  area 
was:  sufficiently  good  that  the  models  should  provide  an  approximation 
of  the  atmospheric  molecular  absorption  for  mean  conditions. 

Wavelength  Dependence  of  Absorption.  As  we  discussed 
previously,  for  monochromatic  transmission,  which  is  the  approximate 
case  for  lasers,  a  line  by  line  calculation  of  absorption  is  necessary. 
For  each  laser  frequency  the  transmission  is  dependent  on  the  parti¬ 
cular  absorption  lines  around  the  laser  frequency,  their  strengths, 
and  their  number.  Because  the  molecular  absorption  is  so  sensitive 
to  the  exact  laser  frequency,  then  the  exact  laser  frequency  must  be 
specified. 

Unless  the  temperature  of  a  laser  can  be  continuously 
monitored,  little  can  be  said  about  the  exact  laser  frequency.  The 
frequency  range  can  and  must  be  determined  for  the  temperature  range 
expected  for  laser  operation.  Since  the  frequency  is  so  important 
In  determination  of  the  molecular  absorption,  a  frequency  must  be 
specified.  Calculation  based  on  the  center  frequency  of  the  expected 


operating  range  is  not  correct  and  may  be  entirely  misleading.  On  the 
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other  hand,  a  line  by  line  calculation  01'  absorption  over  a  specific 
range  and  then  taking  an  average  is  also  incorrect  and  can  lead  to 
errors  much  larger  than  those  induced  by  deviations  from  the  mean 
atmospheric  conditions  as  shown  in  the  tables  in  Appendix  A.  As¬ 
suming  that  the  laser  is  not  equipped  with  a  temperature  or  frequency 
monitoring  device,  then  the  exact  laser  frequency  is  unknown.  In 
order  to  calculate  the  molecular  absorption,  a  worst  case  frequency 
must  be  assumed. 

For  the  purpose  of  this  thesis  a  range  of  wavelengths  for 
the  laser  was  selected  extending  from  1.063y  (9^07 • 3377cm  l)  to 
l.Q65u  (9389. 671A  cm  1)  in  an  attempt  to  make  this  thesis  applicable 
to  as  many  users  as  possible.  Within  this  specified  range,  there  are 
two  molecular  absorbers  with  line  intensities  of  1,0  x  10  27  (cm  V 
molecule  -  cm  2)  or  more  (Ref  50).  A  graph  of  these  intensities 
may-  be  found  in  Appendix  A.  Some  mention,  at  this  point,  must  also 
be  made  about  the  determination  of  these  intensities. 

In  general,  line  intensities  of  the  individual  molecules  are 
very  difficult  to  determine  especially  ?p.  the  atmospheric  environment 
where  the  bands  overlap.  Where  the  data  has  been  supplied  by 
laboratory  measurement  under  vacuum  conditions,  then  the  Intensity 
determination  is  much  more  reliable.  For  a  spectral  region  of  few 
lines  (the  lines  overlap  only  slightly),  then  the  Intensities  may  be 
expected  to  be  quite  reliable  even  for  atmospheric  studies.  There 
is  every  indication  that  the  intensities  within  this  region  of  our 
study  are  accurate;  however,  for  other  regions  the  intensities  may 
be  less  accurate. 


GEP/PH/74-4 


Absorption  at  l,063u.  Calculations  using  the  atmospheric 
models  were  made  for  frequencies  throughout  the  interval  1.063u  to 
1.065P.  Within  this  region,  the  molecular  absorption  at  l.0636p 
(9401.760  cm  *)  was  the  greatest;  that  is,  transmission  at  this 
wavelength  was  the  worst  case.  Calculations  for  the  absorption 
coefficient  for  five  wavelengths  in  the  interval  and  two  just  outside 
the  given  spectral  region  are  shown  in  Appendix  A. 

Genera!  results  of  the  molecular  absorption  calculations 
Include  the  following  points.  The  absorption  increases  as  the 
temperature  and  relative  humidity  increase.  Laser  frequencies  that 
fall  precisely  on  absorber  lines  will  be  affected  most  by  absorption 
in  this  region;  and  those  that  do  not  fall  on  lines  are  less  sensitive 
to  atmospheric  conditions.  Frequencies  that  fall  precisely  on  water 
lines  are  especially  sensitive  to  temperature  and  humidity  variations. 
The  absorption  coefficients  within  the  region  vary  over  three  orders 
of  magnitude  for  the  different  frequencies,  but  for  any  given 
frequency  the  variation  due  to  the  mean  temperature,  pressure,  and 
relative  humidity  models  is  on  the  order  of  a  factor  of  eight  or  less 
for  each,  al titude  considered.  The  extreme  sensitivity  of  absorption 
to  the  laser  frequency  is  shown  most  graphically  by  comparison  of 
the  two  frequencies,  corresponding  to  the  wavenumbers  9388.98  cm  1 
and  9407.07  cm  l,  with  the  worst  of  the  frequencies  within  the 
region,  corresponding  to  the  wavenumber  9401.760  cm  1 . 

Calculation  of  the  Absorption  Coefficient 

As  an  approximation  for  laser  moiecuiar  absorption  coefficients 
for  Germany,  the  model  atmospheres  were  used  as  a  basis  for  the 
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calculations.  The  coefficients  were  determined  using  the  AFCRL  pro¬ 
gram  for  the  worst  case  wavelength  (l.0636p)  for  the  mean  atmospheric 
conditions  (model)  at  each  altitude  in  km.  The  absorption  coefficients 
are  then  averaged  for  the  mean  absorption  coefficient  per  layer  (1  km 
thick).  Next,  the  coefficients  (in  km  *)  are  summed  up  to  the 
altitude  of  interest,  multiplied  by  the  altitude  (in  km)  of  the 
designator  or  aircraft,  and  then  multiplied  by  the  secant  0  to  yield 
the  entire  optical  depth. 

These  calculations  are  designed  to  yield  the  approximate 
molecular  absorption  for  atmospheric  conditions  normally  encountered 
in  Germany.  In  order  to  use  the  approximation  scheme,  no  particular 
meteorological  parameters  are  required  for  table  entry.  For  this 
atmospheric  approximation,  the  molecular  absorption  is  found  to  vary 
only  with  the  season.  In  reality,  however,  large  deviations  from 
mean  conditions  may  be  encountered.  Caution  should  be  exercised  when 
using  this  scheme  for  concurrent  high  temperatures  and  high  humidities. 
Calculations  show  that  for  extreme  conditions  in  summer,  transmission 
may-  be  reduced  by  more  than  fifty  percent. 

in  addition,  this  approximation  scheme  should  not  be 
arbitrarily  extrapolated  to  fit  other  similar  mean  meteorological 
profiles  unless  an  analysis  similar  to  the  one  here  has  been  conducted. 
Large  errors  would  be  induced  by  extrapolating  this  absorption  data 
to  more  tropical  climates.  Also,  in  no  rircumstance  should  this  data 
be  used  for  laser  frequencies  outside  of  this  spectral  region. 
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V.  Atmospheric  Aerosol 


In  addition  to  the  gaseous  molecules,  the  atmosphere  contains 
many  suspended  aerosols.  Aerosols  are  defined  as  dispersed  solid  cr 
liquid  matter  in  a  gaseous  medium,  in  this  case,  air.  If  we  neglect 
cloud,  rain,  and  fog  droplets,  and  consider  only  the  aerosols  found 
in  dry  air,  the  aerosols  range  in  size  from  clusters  of  several 
molecules  to  particles  of  radius  of  20p  as  shown  in  Figure  6  (Ref 
till).  In  spite  of  the  fact  that  aerosol  concentrations  are  orders 
of  magnitude  less  than  molecular  concentrations,  the  aerosols  are 
very  important  in  atmospheric  transmittance  of  radiation.  These 
particles  scatter  and  absorb  radiation,  and  in  addition  they  affect 
the  processes  of  atmospheric  condensation  resulting  in  the  formation 
of  clouds,  snow,  rain,  and  fog. 

This  chapter  will  discuss  atmospheric  aerosol  processes  in 
general  and  how  these  affect  the  scattering  and  absorption  of 
radiation.  This  discussion  is  necessarily  limited  to  those  processes 
that  limit  atmospheric  transmission.  Its  purpose  is  to  provide  a 
minimum  background  from  which  later  conclusions  may  be  drawn  about 
atmospheric  aerosols.  For  a  detailed  aerosol  discussion,  the 
reader  is  referred  to  References  *il  and  hk. 


Aerosol  Attenuation.  The  use  of  Hie  Theory  for  calculations 
of  aerosol  absorption  and  scattering  requires  that  the  shape  of  the 
aerosol  be  a  sphere  and  that  certain  other  aerosol  properties  be 
known. 
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Hie  theory  is  the  exact  determination  of  the  intensity  of  radiation 
in  the  far  field  based  on  the.  aerosol  properties  and  incident 
radiation  (see  Appendix  B) .  Parameters  required  for  this  computation 
are  the  size  and  complex  index  of  refraction  of  the  particle,  and  the 
wavelength  of  the  incident  radiation.  When  extending  this  computation 
to  a  population  of  polydisperse  particles,  the  size  distribution  is 
required  as  well  as  the  particle  number  density,  if  all  particles 
are  not  of  the  same  composition,  then  the  composition  of  particles 
must  be  associated  with  several  particle  distributions  which  are 
additive. 

There  are,  however,  three  important  limitations  to  the  use 
of  this  theory.  The  first  is  that  the  assumption  is  made  that  all 
scattered  radiation  has  the  same  wavelength  as  the  incident  radiation. 
All  quantum  transitions  are  prohibited.  The  second  limitation  is  that 
only  independent  scattering  is  considered.  Scattering  by  a  diffuse 
medium  is-  not  included,  in  order  to  qualify  as  an  independent 
scattering  medium,  the  particles  must  be  arranged  in  such  a  random 
or  inhomogeneous  fashion  that  the  scattering  of  the  medium  is  not 
the  cooperative  effect  of  all  the  particles.  In  addition,  the 
particles  must  be  separated  by  some  minimum  distance.  In  practical 
meteorological  problems,  this  limitation  is  not  restrictive  a.  ; 
even  in  the  most  dense  fogs  the  distance  between  particles  is  twenty 
times  the  radius;  estimates  have  shown  that  separation  distances 
of  three  times  the  particle  radius  are  sufficient  (Ref  62:5).  The 


third  restriction  is  that  the  effects  cf  multiple  scattering  are 


neglected.  This  requires  that  the  energy  removed  from  the  original 
beam  be  equal  to  sums  of  that  removed  by  each  individual  particle; 
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or  In  other  words,  each  individual  particle  is  subject  to  the 
radiation  of  the  original  beam. 

Formation  and  Remova 1  of  Particles  in  the  Atmosphere.  Parti¬ 
cles  found  in  the  atmosphere  are  formed  by  a  variety  of  processes. 

Some  of  these  processes  are  man's  activities,  but  others  are  strictly 
natural  phenomena.  Among  the  important  formation  processes  are:  (l) 
grinding  or  impaction  such  as  produced  by  industrial  activities; 

(2)  breakup  of  liquids  such  as  the  mist  caused  by  sea  spray  or  a  tire 
on  wet  pavement;  and  (3)  condensation  of  particles  in  flames.  Other 
processes  serve  as  apparent  aerosol  sources  but  in  fact  are  part  of 
the  aerosol  evolutionary  cycle:  (1)  coagulation  of  smaller  particles; 
(2)  nucleation,  which  Is  important  in  cloud  formation;  (3)  reactions 
on  the  surfaces  of  particles  resulting  in  particle  growth;  (4)  breakup 
of  large  particles  that  occur  mechanically  or  due  to  phase  change; 
and  C 5)  gaseous  interaction  with  water  droplets  (Ref  8:160).  Such  a 
variety  of  formation  processes  produces  a  v/ide  range  of  particle 
sizes  and  physical  composition. 

Other  processes  remove  particles  from  the  atmosphere.  The 
larger  particles  are  affected  by  gravity  and  settle  to  the  earth's 
surface  if  they  are  not  transported  by  advective  and  convective 
currents.  Particles  of  all  sizes  are  removed  by  impaction  on  plants 
and  other  objects.  In  addition,  particles  are  removed  by  washout 
and  rainout.  Rainout  refers  to  the  loss  of  particulate  matter  by 
condensation  and  formation  of  precipitation.  Washout  refers  to  the 
interaction  between  precipitation  droplet  and  aerosol  particle. 
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Particle  Size  and  Size  Distribution  (General).  Such  a  variety 
of  production  and  removal  mechanisms  leads  to  a  range  of  possible  sizes 
and  number  densities.  The  size  distributions  and  densities  are  con¬ 
trolled  by  the  production  and  removal  mechanisms.  Early  measurements 
of  aerosol  led  to  the  conclusion  that  aerosol  distributions  could  be 
approximated  by  an  equilibrium  distribution  (Ref  59:22).  The  distri¬ 
bution  is  the  result  of  continual  production  and  evolution  of 
particles  and  removal  of  particles  by  the  mechanisms  previously 
discussed. 

Unfortunately  the  measurement  of  aerosol  sizes  and  number 
densities  ?s  extremely  difficult  due  to  the  range  of  particle  size. 
Until  recently, ■ measurement  of  the  entire  distribution  required  that 
several  different  measurement  techniques  be  utilized  simultaneously. 
This  led  to  incomplete  measurements  and  the  mistaken  impression  that 
atmospheric  aerosols  were  composed  of  several  distinct  and  separate 
sizes  of  particles:  Aitken  particles  (r  <  O.lp  ),  large  particles 
(O.lji  <  r  <  l.Op)  and  giant  particles  (r  >  l.Op).  In  addition,  many 
of  these  measurement  techniques  employed  impactor  plates  (often 
mounted  on  aircraft)  and  then  subsequent  analysis  in  the  laboratory. 
These  techniques  often  led  to  inaccuracies  because  they  required  that 
the  aerosol  particle  te  removed  from  its  natural  habitat,  the  air. 

The  method  of  collection  often  caused  physical  alteration  of  particle 
shape  and  the  method  of  analysis,  the  electron  microscope  in  the 
laboratory,  often  caused  chemical  alteration.  Recently,  however, 
improved  instrumentation  described  in  Reference  69  has  given  more 
accurate  measurements  of  aerosol  size  distributions. 
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From  the  various  aerosol  measurements  came  attempts  to  repre¬ 
sent  the  aerosol  distribution  as  an  analytical  function.  There  are 
three  such  general  distribution  forms  in  common  use:  the  power  law, 
the  log-normal,  and  the  modified  gamma  distributions.  Junge  (Ref  4l), 
Bull  rich  (Ref  7),  Whitby  (Ref  68),  and  others  have  demonstrated  that 
the  power  law  distribution  of  the  form 


infer  ■  cr'6  (36) 

is  a  good  approximation  of  aerosol  distribution  in  the  atmosphere 
(Ref  8:17*0.  The  constant  C  is  a  function  of  particle  concentration 
and  the  power  3  describes  the  relative  amounts  of  large  and  small 
particles;  N  is  the  total  concentration  of  aerosol  particles  of  radius 
smaller  than  r.  Because  of  the  wide  range  of  particle  sizes  and  con¬ 
centrations,  it  is  more  convenient  to  use  logarithmic  scales.  Hence, 
the  log  radius-number  distribution  is  defined  by  (Ref  *11:115): 


n(r)  « 


d(logr) 


Using  this  notation,  then  the  power  law  is  most  commonly  expressed 


dnjr)  =  c'r-v 
dr 


(38) 


where  v  -  3  +  1  and  3  has  values  ranging  from  2.5  to  4.0  to  fit 
experimental  data.  On  the  average,  3=3  usually  gives  the  best 
results  for  a  range  between  0.04y  and  lOp  (Ref  44:108).  When  this 
distribution  is  used  to  determine  scattering  and  the  upper  aerosol 
radius  limit  is  ignored  (r2  “)  then  this  compares  quite  favorably 

with  Angstroms  well-known  empirical  formula  for  the  wavelength 
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dependence  of  total  scattering  (Ref  41:142).  In  addition,  this  power 
law  is  most  useful  from  the  standpoint  that  it  describes  the  size 
distribution  by  a  simple  mathematical  expression.  One  disadvantage 
to  the  power  law  distribution  is  that  the  boundaries  must  be  specified 
and  have  a  strong  influence  on  the  integrations  and  scattering 
calculation . 

A  more  complex  distribution  is  the  log-normal  distribution; 
it  requires  a  third  parameter,  standard  deviation,  in  addition  to  the 
two  required  by  the  power  law  distribution.  This  distribution,  how¬ 
ever,  has  the  advantages  that  it  can  handle  large  ranges  of  particle 
sizes  and  it  represents  the  extreme  radii  values  more  accurately 
In  the  integration  process.  This  distribution  also  represents  the 
particle  sizes  produced  by  several  comminution  processes  in  nature. 

The  normalized  log-normal  function  is  usually  written 


ink)  =( - J - \ 

d(logr)  \(2ti)*  log  ag  /e 


(logr-logrQ)2 
2  log2  ag 


(39) 


where  ag  and  rg  are  the  geometric  standard  deviation  and  mean  radius 
respectively  and  dn(r)  is  the  normalized  number  of  particles  between 
radii  r  and  r  +  dr  (Ref  8:1 69). 

The  modified  gamma  distribution  has  also  been  used  to  describe 
wet  and  dry  aerosol  size  distributions  (Ref  17  and  18).  The  wet 
aerosols  include  rain,  hail  and  clouds  and  the  dry  aerosols  are  haze. 
In  general,  the  modified  gamma  {s  mathematical ly  as  advantageous  as 
the  log-normal  distribution,  but  the  modified  gamma  has  not  been 
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linked  to  any  specific  natural  processes.  In  addition,  it  has  more 
parameters  which  are  not  relatable  to  any  aerosol  property.  The 
modified  gamma  distribution,  however,  offers  the  advantage  that  the 
one  generalized  distribution  describes  the  various  meteorological 
phenomena.  This  family  of  distributions  is  described  by 


n(r)  «  ara  exp(-bra) 
or 

dn(r)  (a“l).  .  a.  ,  .  o, 

-jp-  =  ar  (ct-obr  )exp(-br  )  (40) 

where  a,  b,  and  a  are  positive  real  constants  and  a  is  a  positive 
integer  (Ref  17:15,76). 

Although  there  are  other  distributions  used  to  describe  aerosol 
measurements,  these  three  are  the  most  common  and  most  widely  accepted 
in  describing  aerosol  distributions.  The  distributions  will  now  be 
discussed  in  terms  of  applicability  to  measured  results  and  physical 
considerations. 


Particle  Size  and  Size  Distribution  (Specific).  The  particle 
size  distribution  and  number  density  is  controlled  by  the  particular 
production  and  removal  mechanisms.  This  statement  necessarily 
implies  that  if  there  are  different  production  mechanisms  associated 
with  various  geographical  locations,  then  it  can  be  expected  that 
there  would  not  be  one  worldwide  aerosol  distribution  function.  In 
fact,  the  distribution  function  at  one  geographical  location  would  be 
to  some  degree  a  function  of  the  local  aerosol  production  mechanisms. 

If  it  were  not  for  the  constant  mixing  of  the  atmosphere,  then 
one  or  more  generalized  mathematical  functions  would  not  be  descriptive 
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of  the  atmospheric  aerosol.  In  fact,  description  of  atmospheric 
aerosol  would  be  a  point  by  point  computation.  Fortunately,  there 
is  continual  smoothing  of  the  aerosol  distribution  function  due  to 
the  transport  of  aerosols  and  continual  atmospheric  modification 
processes.  Nevertheless,  local  aerosol  sources  are  important  but 
not  for  large  scale  distributions. 

There  is,  for  example,  a  difference  in  the  general  aerosol 
distribution  and  number  density  found  over  a  large  land  mass  and  that 
found  over  the  oceans.  Because  of  advective  and  convective  currents, 
there  is  no  sharp  division  between  the  maritime  distribution  and 
continental  distributions,  in  fact,  aerosols  of  maritime  origin  have 
been  detected  thousand  of  miles  inland;  likewise,  continental  particles 
are  measured  in  the  center  of  the  oceans  (Ref  41:176).  In  like  manner, 
an  area  with  a  large  amount  of  industrial  activity  would  have  a  dif¬ 
ferent  aerosol  size  distribution  and  number  density  than  a  region  that 
was  free  of  such  activity.  On  the  smaller  scale,  the  area  around  a 
large  industrial  city  may  be  expected  to  have  a  slightly  different 
distribution  than  a  rural  area  located  some  distance  away,  though 
mixing  will  provide  continual  smoothing  of  the  different  distributions 
and  particle  number  densities. 

There  are  two  distributions  that  have  been  used  extensively 
to  describe  global  dry  aerosol  distributions.  They  are  the 
Deirmer.djian  Continental  and  Maritime  Haze  Distributions  and  variations 
of  these  original  models.  The  Continental  Haze  model  uses  the  power 
law  with  exponent  of  four.  The  number  of  large  particles  compared 
to  the  small  particles  is  described  by  the  exponential  power.  The 
Maritime  model  is  a  modified  gamma  distribution.  It  describes  the 
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aerosol  distribution  of  extremely  maritime  air  and  its  use  over  land 
areas  is  limited  to  coastal  regions.  In  general,  the  differences  in 
the  air  masses  are  as  follows:  (1)  the  number  of  particles  in  the 
martime  air  is  less  (usually  substantially  less)  than  the  continental 
atr;  (2)continental  air  has  more  Aitken  particles  by  orders  of  magni¬ 
tude;  and  (3)  a  greater  portion  of  the  maritime  particles  than 
continental  particles  are  in  the  large  particle  range.  The  existence 
of  two  generalized  distributions,  maritime  and  continental,  as  well 
as  these  general  distribution  characteristics, is  well  accepted.  Such 
model  distributions  are  useful  quantitatively,  but  it  is  important 
to  realize  that  these  are  characteristic  distributions  and  that  in 
specific  instances,  conditions  may  deviate  considerably  from  these 
mean  distributions. 

There  are  two  important  conclusions  that  may  be  drawn  from  this 
discussion.  One  is  that  aerosol  distributions  ir.  the  atmosphere  are 
affected  by  local  aerosol  sources.  The  other  is  that  aerosol  distri¬ 
butions  tend  to  smooth  out  yielding  generalized  aero'cl  distributions 
for  larger  geographical  regions  provided  the  air  mass  circulates. 

Point  sources  of  aerosols  are  still  important,  but  if  there  is  suf¬ 
ficient  circulation  a  smoothing  or  quasi-steady  state  distribution  is 
approached.  Regardless,  it  becomes  apparent  that  geographical  location 
and  air  mass  circulation  are  extremely  important  in  determining  aerosol 
size  distribution  and  number  of  particles. 

Particle  Composition.  Aerosol  chemical  composition  is  important 
in  determining  the  complex  refractive  index  for  Mie  Theory  calculations. 
The  chemical  composition  of  the  aerosoi  is  determined  by  the  source  and 
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evolutionary  process  of  the  aerosol.  Several  extensive  investigations 
have  been  made  of  the  effect  of  both  of  these  processes  on  the  aerosol 
composition.  These  investigations,  like  the  size  measurements,  have 
been  hampered  by  the  difficult  nature  of  the  task  of  measuring 
compos i t i on . 

One  of  the  difficulties  of  determining  composition  is  the  in¬ 
ability  to  measure  composition  of  the  smaller  (Aitken)  particles.  In 
general,  the  composition  of  large  and  giant  particles  is  measured  and 
the  composition  of  the  smaller  particles  is  assumed  to  be  similar. 

Several  workers  have  measured  the  distillate  from  rainwater  or  other 
forms  of  precipitation  and  assumed  this  measurement  was  representative 
of  the  dry  aerosol.  This  method  also  fails  to  measure  the  small 
particle  composition,  because  condensation  nuclei  are  limited  in  large 
to  the  giant  and  large  particles  although  there  may  be  some  small 
particle  contribution  due  to  washout.  Despite  these  limitations,  a 
remarkable  feature  of  these  investigations  is  that  their  findings  are 
surprisingly  similar.  In  addition,  it  is  the  large  and  giant  particles 
that  affect  the  transmission  of  radiation  in  the  visible  and  infrared 
regions. 

Aerosol  is  composed  of  particles  from  a  variety  of  sources. 

Average  global  aerosol  is  composed  of  the  following  general  classes 
of  compounds:  Sulfur,  Nitrogen,  Chlorine,  soils,  and  combustion 
products  (Ref  59:12).  Table  IV  shows  estimates  of  the  magnitude  of 
particle  formation  (r  <  2up)  for  the  surface  layer  of  the  atmosphere. 

The  percentage  composition  of  these  compounds  in  various  aerosol 
measurements  is  found  to  differ  from  location  to  location.  Aerosol 
composition  over  ocean  areas  is  found  to  have  a  large  chloride  component. 
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Table  IV 

Estimates  of  Particle  Formation  In  Millions  of  Metric 
Tons  Per  Year  for  Particles  <  20\i  Radius 
(Ref  140:183) 


Natural 

Soil  and  rock  debris3 

100-500 

Forest  fires  and  slash-burning  debris3 

3-150 

Sea  salt 

(300) 

Volcanic  debris 

25-150 

Particles  formed  from  gaseous  emissions: 

Sulfate  from  H?$ 

130-200 

Ammonium  salts  from  NMj 

80-270 

Nitrate  from  NO 

X 

60-A3C 

Hydrocarbons  from  plant  exudations 

75-200 

Subtotal 

773-2200 

Man-Made 

Particles  (direct  emissions) 

10-90 

Particles  formed  from  gaseous  emissions: 

Sulfate  from  S02 

130-200 

Nitrate  from  N0X 

30-35 

Hydrocarbons 

15-90 

Subtotal 

185-415 
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The  existence  of  chlorides  is  not  limited  to  maritime  areas,  however, 
nor  does  it  suggest  a  particle  of  maritime  origin  since  industrial 
activities  have  been  linked  to  chloride  aerosol  production  (Ref  41: 
170). 

The  physical  constitution  of  aerosol  particles  in  the  atmos¬ 
phere  varies  between  the  extremes  of  a  dry  insoluble  dust  particle 
and  a  droplet  of  completely  soluble  material.  In  general,  though, 
particles  consist  of  mixtures  of  many  materials,  both  soluble  and 
insoluble.  Throughout  the  particle  lifetime,  it  is  subject  to  the 
influence  of  other  particles  through  the  processes  of  coalescence  and 
coagulation.  It  is  also  subject  to  many  atmospheric  gases,  some  of 
them  trace  gases,  which  are  known  to  play  an  important  role  in  aerosol 
evolution.  In  fact,  it  is  now  believed  that  perhaps  the  majority  of 
particles  in  suspension  in  the  atmosphere  are  secondary  products  formed 
from  material  which  entered  the  atmosphere  as  a  gas  (Ref  40:192).  Some 
idea  of  the  complexity  of  the  problem  may  be  illustrated  by  Table  V. 

Since  the  particle  is  changing  in  time  and  can  be  affected  by 
any  of  a  number  of  processes,  and  generally  its  definite  origin  is 
unknown,  the  concept  of  a  mixed  particle  evolved.  This  is  the  most 
general  concept  for  an  atmospheric  aerosol  particle.  This  concept 
attaches  an  average  composition  or  average  index  of  refraction  to  the 
entire  distribution  of  particles.  The  value  of  this  concept  is  the 
limitation  it  places  on  determining  a  distinct  aerosol  distribution 
for  each  chemical  composition. 
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Table  V 

Possible  Atmospheric  Gases  From  Which  Particles  May  Be 
Formed  Through  Secondary  Reactions 
(from  Ref  8:4,5) 


Symbol 

Name 

n2 

Ni trogen 

°2 

Oxygen 

**°A  r 

Argon 

co2 

Carbon  dioxide 

He 

Neon 

He 

Hel ium 

CH4 

Methane 

Kr 

Krypton 

H2 

Hyd rogen 

Xe 

Xenon 

CO 

Carbon  monoxide 

n2o 

Nitrous  oxide 

so2 

Sulfur  dioxide 

NO 

Nitric  oxide 

no2 

Ni trogen  dioxide 

HCHO 

Formaldehyde 

nh3 

Ammonia 

o3 

Ozone 
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Complex  Index  of  Refraction.  Since  particles  are  being  formed 
by  g.  s  conversion  and  are  evolving  through  gas  absorption,  conden¬ 
sation  and  coagulation,  only  a  relatively  few  measurements  have  been 
made  from  which  the  index  of  refraction  can  be  inferred  (Ref  59:15). 
Bull  rich  (Ref  41 : 1 1 4)  calculated  a  real  index  of  refraction  of  1 .5*» 
based  on  assumed  atmospheric  constituents.  Volz  (Ref  65:822)  computed 
a  real  refractive  index  of  1.53  based  on  distillation  residues  from 
rain  water.  Recent  measurements  shown  in  Table  VI  seem  to  indicate 
that  these  earlier  determinations  of  the  real  refractive  index  were 
in  fact  accurate.  The  wavelength  dependence  of  the  refractive  index 
(real  and  imaginary  part)  is  illustrated  in  Table  VI. 

The  determination  of  the  imaginary  (absorptive)  part  of  the 
refractive  index  is  less  well  established.  Values  for  the  absorptive 
part  at  l.Oy  range  from  .072  (Fisher  1971)  to  approximately  .015 
(Ref  66)  as  shown  by  comparing  the  figures  in  Table  VI.  The  deter¬ 
mination  of  the  complex  index  is  extremely  difficult  and  errors  of 
from  10%  to  an  order  of  magnitude  can  be  expected  in  these  values 
(Ref  59:15). 


Particle  Growth.  Some  of  the  difficulty  in  determining  aerosol 
index  of  refraction  and  size  distribution  is  attributable  to  the  growth 
of  atmospheric  aerosol  with  increase  in  relative  humidity.  The  water 
vapor  affects  the  dry  particle  in  two  ways.  One,  it  condenses  on 
the  particles  thus  causing  the  individual  particles  to  increase  in 
size.  This  growth  of  the  particles  affects  the  entire  size  distribution 


second ,  the  water  vapor  condensing  on  the  particle  causes  a  change 


in  the  complex  index  of  refraction. 


/  I  II/  /  *1  I 


”1 


in 

*j 

u 

TO 

Cu 


>. 

J_ 

ta 

c 


D) 

id 


E 


~a 

c 

ID 


«D 

Q) 

CC 


in 

D) 

u 


•a 

c 


4) 

> 


+J 

O 

ID 

U 


4) 

CC 


U 

L. 

4) 

4) 

4) 

03 

u 

►*• 

JZ 

*—  -C 

|M 

-C 

L. 

4) 

u 

4}  O 

4) 

o 

N 

Z3 

c 

v> 

c  in 

C 

W 

f— 

A3 

*— 

03  »— » 

03 

0 

to 

X 

u. 

X  u_ 

m 

llm 

> 

vO 

cn 

* 

•* 

•K 

vp 

*— 

fA  CM 

cn 

0mm 

CM 

r-s 

.  o 

VO 

LA  1%, 

-y 

Cv| 

m 

• 

• 

o 

•  o 

• 

O 

• 

o 

• 

• 

• 

<~s 

Ifl 

C 

o 

v-  o 

•— 

CM 

IA  CO 

cn 

O 

CM 

so 

o  o 

VO 

LA  vo 

CM 

m 

«M*  • 

• 

o 

•  O 

• 

o 

• 

o 

r-»  • 

0mm 

• 

0mm 

X 

4-J 

O) 

c 

V 

4) 

r— 

co 

-y  ia 

o 

VO 

cn 

>  o 

vO 

to 

LA  LA 

LA 

LA 

o 

ID  fv. 

• 

o 

•  O 

• 

o 

• 

o 

?  * 

• 

• 

• 

CM 

p— • 

LA  CO 

LA 

m 

vo 

m 

VO 

to 

la  ,<r 

LA 

»—• 

in 

o 

in 

• 

o 

•  o 

• 

O 

• 

o 

• 

r— • 

• 

•y 

-3- 

Is-  <V| 

m 

m 

cn 

in 

O 

vO 

•N* 

LA  -3" 

LA 

LA 

o 

• 

O 

•  O 

• 

© 

• 

o 

•1 

• 

#*• 

• 

• 

• 

• 

• 

• 

C7> 

—  CD 

O) 

o> 

ro 

03 

ID  ID 

ID 

03 

03 

03 

4) 

E 

4>  E 

4) 

E 

03 

E 

as 

a:  — 

DC 

a: 

cn 

L. 

>* 

>* 

4) 

c 

c 

.a 

ID  VO 

03 

c 

E  vO 

E 

a> 

o 

km 

in 

u 

IA 

a> 

4) 

W 

cn 

C3 

43  O 

—m 

•<c 

*— 

a 

U"\ 

VO  4) 

>» 

•u 

03 

•V 

l* 

"  CL 

O. 

c  u 

M 

6\ 

M  PW 

— .  e~ 

m  *% 

m 

U  r> 

0) 

c 

B 

«  c 

4)  4) 

e  i 

j5 

0)  3 

*0 

«— 

E 

•o  —  c 

•o  X 

u  E 

o 

4J  — 

o 

A3 

D 

O  ID  CU 

o  o 

4J  3 

03  n 

X  X  to 

X  X  “> 

X  X 

O  U) 

u  3  in 

•a 

«J 

•a 

c 

4) 

•M 

X 

4) 

>• 


ID 

V 

C 


4) 

4) 

-Q 

4) 

> 

ID 

X 

VO 

O 


L. 

O 

4- 

01 
4) 
3 
* — ■ 
ID 
> 


3 

O 

o 


>* 

c 

o 

•o 

c 

4) 

+-> 

X 

o 

u 

4) 

X 

O 

in 

IT 

T3 

C 

ID 


4> 

C 

ID 

X 


o 

ID 


m 

O 


4; 


c 

ID 

■w 

-Q 

o 

o 


-M 

>» 


u 

ID 

in 

in 


4) 

O 

4) 

C 

in 

ID 


5 


c 

o 


ID 

o 

Q. 

l_ 


4) 

■u 

c 


o 

o 


o 

r>» 


LO 

LA 


o 

-y 


4-> 

ID 

•M 

o 

c 

4) 

u 

ID 

N 

O 

> 

>» 

-O 

•o 

4) 

X 

in  • 

—  in 

—  4) 
X  3 
3  — 
Q.  ID 

> 

in 

4)  4) 
3  in 

#—  (U 
m  _r^ 

>  *•* 

•t: 

•* 


X 

) 


! 


66 


GEP/PH/7A-A 


Present  studies  indicate  that  particle  growth  is  determined 
by  a  number  of  factors.  One  of  the  most  important  factors  is  the 
original  particle  composition — soluble  or  insoluble  material  or  mixture. 
Another  is  the  particle  shape,  although  the  shape  of  the  atmospheric 
aerosol  is  usually  assumed  to  be  spherical.  Another  factor  affecting 
particle  growth  is  the  relative  humidity  change  and  and  starting  point 
(size) . 

Particle  growth  raises  a  number  of  questions  in  terms  of  Hie 
theory  calculations.  Three  questions  that  must  be  considered  are: 

(1)  How  ooes  the  index  of  refraction  change  v/ith  relative  humidity? 

(2)  How  is  the  particle  size  distribution  affected  by  particle  growth? 
and  (3)  How  is  original  composition  important  to  particle  growth? 

The  first  and  third  questions  may  be  answered  together  because 
the  particle  composicion  is  directly  related  to  the  index  of  refract¬ 
ion.  There  is  a  range  of  particle  compositions  between  totally  water 
soluble  material  and  insoluble  material.  Consider  the  case  of  the 
water  soluble  material  first.  As  water  condenses  on  the  outer  surface 
of  the  particle,  a  shell  of  water  is  formed  on  the  sphere.  As  the 
relative  humidity  increases,  a  point  is  reached  at  which  the  vapor 
pressure  of  the  aqueous  solution  at  the  surface  is  less  than  the  vapor 
pressure  of  water  in  the  ambient  air.  At  this  point,  water  vapor  from 
the  air  is  collected  until  the  substance  is  dissolved  and  in  equili¬ 
brium  with  the  air  (Rer  8:222).  During  this  period  the  aerosol 
experiences  a  tremendous  growth.  This  description  ;s  valid  for  a  pure 
soluble  material  such  as  sodium  chloride.  This  process  is  called 
deliquescence.  If  aerosols  v/erc  all  pur?,  soluble  materials,  this 
particle  growth  would  be  a  very  complex  phenomena.  In  fact,  the 
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measurement  of  growth  rates  that  are  less  rapid  than  pure  soluble 
materials  supports  the  concept  of  the  mixed  particle. 

On  the  other  hand,  the  dry  aerosol  particle  may  be  insoluble 
and  grow  very  little.  In  any  case,  the  question  of  how  to  deal  with 
the  particle  as  water  is  condensed  about  it  remains.  There  are  two 
distinct  possibilities.  One  involves  working  the  problem  of  a  sphere 
covered  with  a  spherical  water  shell.  In  other  words,  a  two  step  pro¬ 
blem  with  separate  indices  of  refraction.  The  second  involves  working 
the  one  sphere  problem  but  changing  the  overall  refractive  index  as 
the  sphere  grows  in  size. 

This  second  method  is  more  acceptable  to  the  mixed  particle 
concept  (not  to  mention  the  increased  difficulty  in  problem  solution 
if  the  first  method  is  utilized),  it  is,  however,  important  to 
realize  the  two  important  conditions  that  must  be  fulfilled  in  order 
to  ignore  shell  structure  for  haze  calculations:  (l)  the  mixing 
ratios  of  the  components  must  differ  only  slightly;  and  (2)  the  radius 
of  the  particle  must  approximate  the  wavelength  of  the  incident 
radiation  (Ref  44:115).  That  is,  if  these  two  conditions  are  met,  We 
assume  the  refractive  index  can  be  determined  from  the  mixing  ratio 
of  its  components. 

The  second  major  question  is  just  as  complex.  Since  the 
particle  sizes  are  different  they  grow  at  different  rates  because  of 
the  difference  in  curvature  at  their  surfaces.  If  their  compositions 
were  different  then  the  problem  would  be  compounded.  Assuming  that 
the  particles  are  mixed  in  composition  allows  an  approximation  to  the 
increase  in  the  size  distribution.  If  no  particles  are  being  added 
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or  removed  from  the  distribution,  then  for  relative  humidities  up 
to95%the  growth  rate  by  particle  size  is  rather  uniform.  This  is 
borne  out  by  experimental  measurement  as  shown  in  Table  VII.  Beyond 
35%  relative  humidity  (R.H.)  the  larger  parti clesgrow  more  rapidly, 
and  hence  the  large  and  giant  particles  become  the  fog  and  cloud 
droplets. 


Table  VII 

Particle  Growth  With  Relative  Humidity  (RH) : 

Ratio  r/ro  of  Radii  of  Wet  (r)  and  Dry  (r0)  Aerosol  Particles 
Radii  (r0)  Given  in  Microns  (Ref  7-5^7) 


rQ(ym) 

0.04 

0.1 

l 

10 

RH 

r/ro 

0 

1 

1 

1 

1 

0.3 

1.038 

1  .o4o 

1  .042 

1.042 

0.6 

1.153 

1.159 

1.165 

1.165 

0.8 

1.400 

1.414 

1.426 

1 .426 

0.9 

1.656 

1.695 

1.720 

1.720 

0.95 

1.93 

2.01 

2.05 

2.06 

0.975 

2.29 

2.48 

2.58 

2.59 

0.99 

2.79 

3.19 

3.45 

3.48 

0.995 

3.13 

3.83 

4.30 

4.36 

0.998 

3.48 

4.69 

5.81 

5.91 

Yet  another  particle  growth  mechanism  is  coagulation.  This 


process  effectively  forms  the  lower  limit  of  particle  sizes.  In  a 
stagnant  or  equilibrium  case  where  the  production  of  small  (Aitken) 
particles  has  ceased,  the  particle  size  distribution  is  shifted  to  the 
larger  particle  sizes  as  shown  in  Figure  7.  This  growth  mechanism  was 
convincingly  demonstrated  in  the  comprehensive  1969  Los  Angeles  Smog 
Experiment  (Ref  68:386-196).  This  growth  mechanism  is  not  as  important 
as  it  might  seem  because  it  affects  largely  the  inefficient  scatterers 
for  our  purposes. 
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Fig.  7*  Particle  Size  Distribution  Shift  Due  to  Coagulation 
Superscript  h  refers  to  hours,  d  to  days  (after  Ref 
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Vertical  Distribution.  The  particle  size  distributions  tend 
to  remain  constant  with  increase  in  altitude  except  for  particles 
with  radii  larger  than  1y..  A  decrease  in  percentage  concentration 
of  these  larger  particles  is  noticeable  in  the  absence  of  convective 
mixing  currents  that  would  force  them  aloft  (Ref  41:127).  In  general, 
there  is  a  decrease  in  particle  concentration  with  increase  in 
altitude,  and  this  is  rather  uniform  for  all  particle  sizes  (Ref  41:128). 

The  vertical  aerosol  distribution,  however,  is  affected  by  a 
number  of  factors.  One  of  the  most  important  is  the  aerosol  size  dis¬ 
tribution  at  the  surface  and  the  particular  production  mechanisms. 
Another  is  the  origin  and  history  of  the  air  mass.  Others  are  the 
vertical  temperature  and  pressure  profiles  and  the  strength  of  the 
vertical  mixing. 

Since  most  particle  sources  are  on  land  masses  then  the  lower 
layers  of  the  continents  are  the  most  polluted  (naturally  and  arti- 
fically)  in  terms  of  aerosol  concentrations.  The  vertical  extent  of 
this  layer  over  the  continental  areas  varies  by  season  and  geographical 
region.  The  height  or  extent  of  the  mixing  layer  depends  on  existing 
and  past  meteorological  conditions.  It  has  been  well  documented  how 
sensitive  particle  distributions  reflect  the  vertical  temperature 
distribution.  Haze  layers  beneath  temperature  inversions  are  a 
frequent,  and  quite  striking,  phenomena  (Ref  4i:l83). 

Over  ocean  areas  and  above  tho  continental  mixing  layer  the 
aerosol  concentrations  are  noticeably  smaller  than  in  the  mixing  layer, 
in  fact,  the  distribution  above  the  mixing  layer  is  fairly  constant 
and  differs  only  from  the  maritime  distribution  in  numbers  of  large 
particles  (Ref  40.202).  This  background  distribution  is  representative 
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of  the  troposphere  above  the  mixing  layer  (over  ocean  and  continent) 
and  is  considered  a  typical  very  clean  air  distribution.  The  vertical 
distribution  of  aerosols,  however,  is  the  subject  of  a  more  lengthy 
discussion  later  in  this  paper. 
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V ! .  Aerosol  Attenuation  for  Horizontal  Path 
at  Ground  Level 

The  previous  chapter  discussed  the  physical  aerosol  mechanisms 
and  their  relationship  to  aerosol  attenuation.  The  purpose  of  this 
chapter  is  to  discuss  aerosol  attenuation  models  for  application  to  a 
specific  geographical  region.  Existing  models  are  evaluated  on  the 
basis  of  their  suitability  for  use  in  describing  aerosol  attenuation 
over  continental  Germany.  Aerosol  attenuation,  as  predicted  by  the 
model,  is  then  linked  to  meteorological  observables. 

Particle  Size  Distribution 

Several  different  particle  size  distributions  have  been  used 
successfully  by  various  authors  to  explain  measured  attenuation  values. 
In  general,  all  distributions  are  the  result  of  aerosol  attenuation 
measurements.  Particle  size  distributions  are  chosen  on  the  basis  of 
the  best  analytical  fit  to  data.  Other  quantities  such  as  the  number 
of  particles  and  the  refractive  index  of  the  particles  have  been  dif¬ 
ficult  to  obtain  for  reasons  discussed  in  Chapter  V.  In  addition, 
these  quantities  are  not  always  measured  concurrently  with  the 
attenuation,  and  therefore  determination  of  a  wxique  size  distribution 
has  been  quite  complicated. 

An  example  of  this  difficulty  is  the  set  of  attenuation 
coefficients  described  in  Reference  1 4 .  This  Naval  Research  Laboratory 
(NRL)  project  involved  attenuation  measurements  taken  over  two  over¬ 
water  paths  in  the  Chesapeake  Bay  area  from  April,  1959  >  to  January, 
I960,  for  various  meteorological  conditions.  Figure  8  shows  typical 
reported  experimental  scattering  curves  from  this  report. 
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The  difficulty  in  describing  analytically  this  set  of  measurements  by 
one  aerosol  size  distribution  is  apparent  from  the  figure.  Reference 
72  analyzes  the  thirty-six  reported  observations  from  this  project 
and  classifies  the  resultant  size  distributions  in  no  less  than  five 
categories. 

To  some  degree,  this  difficulty  has  been  exaggerated,  for  the 
NRL  project  was  accomplished  over  a  large  body  of  water  and  hence  in 
an  area  where  a  mixing  of  two  general  distributions,  maritime  and 
continental,  is  common,  in  fact,  one  of  the  findings  of  the  report 
was  that  the  aerosol  size  distribution  above  or  near  water  can  generally 
be  considered  as  having  a  Junge-type  distribution  plus  a  small  compo¬ 
nent  of  larger  particles  (Ref  15:551). 

For  most  of  Germany,  however,  the  assumptions  of  a  generalized 
continental  particle  size  distribution  is  realistic  (Ref  30).  Areas 
of  Germany  that  border  large  bodies  of  water  such  as  the  Baltic 
Coastline  are  likely  to  have  distributions  that  are  similar  to  the 
composite  distributions  measured  at  Chesapeake  Bay  For  the  remainder 
of  Germany,  the  aerosol  size  distribution  may  be  described  approxi¬ 
mately  using  a  form  of  power  law  described  in  Chapter  V.  This 
generalized  distribution,  of  course,  is  modified  where  local  aerosol 
sources  abound  (large  cities  or  industrial  complexes). 

A  study  of  general  air  mass  movement  and  source  regions  for 
Europe  supports  the  theory  of  a  continental  aerosol  d ? s -  -  "'vtion.  A 
study  of  European  air  masses  and  their  source  regions  with  particular 
emphasis  on  Germany  was  conducted  by  the  German  meteorologist  Schinze 
and  described  in  Reference  9.  Even  the  maritime-polar  air  (mp)  which 
is  a  predominant  air  mass  in  central  and  western  Europe  is  extensively 
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modified  by  transport  across  land  areas  (Ref  9:282). 

Additional  evidence  which  supports  the  use  of  a  continental 
distribution  is  the  set  of  attenuation  models  and  measurements 
reported  by  Reference  7  and  Reference  3?i  for  locations  in  Germany. 

These  measurements  were  modeled  using  a  straight  power  law  and  a 
Deirmendjian  type  continental  haze  model.  The  use  of  both  will  be 
discussed  in  the  next  section. 

Analysis  of  Size  Distributions 

If  a  relationship  Is  desired  between  a  meteorological 
observation  recorded  in  the  visible  and  attenuation  at  another 
wavelength,  then  using  Mie  theory  and  an  assumed  aerosol  size  distri¬ 
bution  the  relationship  is  approximately  determined.  The  use  of  the 
word  approximately  is  important  here,  for  the  index  of  refraction 
as  we  shall  later  see  does  have  some  effect  on  this  relationship. 

Additionally,  as  long  as  the  number  o*  particles  is  unknown, 
then  a  theoretical  determination  of  aerosol  attenuation  using  Mie 
theory  yields  only  a  relative  scattering  coefficient.  Since  the  total 
particle  density  is  essentially  unobtainable  even  with  today's  improved 
equipment,  then  each  size  distribution  is  normalized  such  that  the 
attenuation  is  the  attenuation  given  per  particle  per  unit  volume. 

This  concept  is  used  throughout,  this  analysis  of  the  different  particle 
size  distributions.  For  the  present,  the  total  particle  density  is 
unimportant.  The  important  factors  to  be  analyzed  are  the  size  dis¬ 
tribution  (with  respect  tc  the  wavelengths  of  interest)  and  the 
particle  index  of  refraction. 

Using  exact  Mie  theory  calculations  as  described  in  Appendix 
B,  the  aerosol  attenuation  coefficients  per  particle  per  cm3  per  km 
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path  length  were  computed  for  wavelengths  from  .40y  to  1.0636y. 
Assumed  indices  of  refraction  for  the  different  wavelengths  are  given 
in  Table  Vi!i. 


Table  VI 1 1 

Assumed  Indices  of  Refraction  Values  Obtained 
by  Taking  Mathematical  Average  of  Models  5  and  6  in  Table  VI 


Wavelengths  in  Microns 


.hO 

.55 

.70 

1.00 

1.0636 

1.55 

1.53 

1.52 

1.51 

1.51 

.022 

.031 

.035 

.0  44 

.OhS 

R  *=  Real  Part 
!  Imaginary  Part 


Using  these  indices,  aerosol  attenuation  coefficients  were  compared 
for  the  different  distributions.  The  different  size  distributions 
investigated  were:  (l)  Deirmendj ian  Maritime  Haze;  (2)  Deirmendjian 
Haze  C;  (3)  A  combination  distribution  using  varying  contributions  of 
the  two  components --Mari time  Haze  and  Haze  C;  (*0  The  modified 
Diermencjian  Haze  C  used  by  AFCRL;  (5)  Deirmendjian  Haze  L;  and  (6)  a 
straight  power  law.  The  exact  distributions  used  are  given  in  Table 


GEP/PH/7^ 


re- ^rjT'TT.'S  TT^y.’KWSYITE' "  f^f-^  ^  e*^ 


'**"  v^p 


GEP/PH/7^"1* 


The  upper  and  lower  bounds  of  integration  were  chosen  as  20y  and 
.003n  respectively  except  in  the  case  of  the  power  law  where  .03  was 
used  as  the  lower  limit.  The  published  Dei rmendj ian  Haze  C  distri¬ 
bution  was  used  except  the  upper  bound  was  extended  to  20y. 

Comparison  of  the  distributions  yields  the  following  general 
characteristics  as  shown  in  Figure  9.  A  100%  Maritime  distribution 
gives  the  highest  per  particle  attenuation;  this  attenuation  increases 
as  the  wavelength  is  increased  from  the  visible  to  near  IR  then  begins 
to  drop  off  around  ly  such  that  the  aerosol  attenuation  at  1.0636p 
is  approximately  equal  to  the  attenuation  at  ,55y.  As  the  percentage 
of  maritime  contribution  is  reduced  and  the  continental  increased  then 
the  attenuation  at  1.0636y  becomes  less  than  that  at  .55y>  although 
this  decrease  is  not  significant  until  the  percentage  maritime  is 
less  than  10%.  This  suggests  that  for  regions  such  as  coastal  regions 
and  ocean  areas  the  aerosol  attenuation  at  1 .0636y  may  not  be  signifi¬ 
cantly  less  than  that  in  the  visible  spectrum. 

The  attenuation  per  particle  increases  by  approximately  one 
order  of  magnitude  from  a  distribution  containing  10%  maritime  and 
90%  continental  to  a  distribution  of  100%  maritime.  Haze  L  v/as  found 
to  approximate  a  distribution  of  about  ^0%  maritime  and  60%  conti¬ 
nental  with  the  exception  that  there  was  a  noticeable  slope  to  the 
curve.  That  is,  attenuation  was  less  at  1.0636y  than  at  .55y» 

A  straight  power  law  was  used  from  .03y  to  20y  to  compare  the 
attenuation  coefficients  against  those  computed  using  Deirmendjian 
Haze  C.  The  resulting  attenuation  coefficients  were  approximately  one 
order  of  magnitude  less  than  those  computed  for  Haze  C  but  the  slope 
of  the  resulting  curve  v/as  essentially  identical  to  the  Haze  C  slope. 
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Other  powers  were  not  Investigated  by  this  author,  but  such  an  inves¬ 
tigation  has  been  made  in  Reference  14  and  results  are  shown  in  Figure 
10.  The  power  r  =  -3  in  this  figure  corresponds  to  the  integrated 
power  law  that  is.  the  distributions  ate  identical  with  the 

exception  of  the  tower  limit  of  .Olp  used  in  the  URL  investigation. 
Tiere  is  the  additional  difference  of  the  index  of  refraction  which 
for  the  curves  in  Figure  10  is  assumed  to  be  i. 33  with  no  imaginary 
part*  Although  the  curves  from  Figure  10  do  not  yield  attenuation 
coefficients  that  are  exactly  cmparoble  with  those  computed  here, 
their  importance  is  that  they  display  the  relative  change  with  respect 
to  wavelength  of  the  attenuation  coefficient  as  the  power  is  changed 
on  the  size  distribution. 

An  additional  analysis  involved  the  change  of  the  lower  limit 
of  the  Haze  C  distribution  to  .02p.  The  result  was  a  subsequent 
reduction  in  the  attenuation  coefficients  which  Is  expected  since  the 
distribution  is  normalized.  There  was,  however,  no  significant  change 
in  the  slope  of  the  attenuation  coefficient*  curve  as  shown  in  Figure 
9. 

Two  distributions  were  used  to  check  the  effect  of  changing 
the  upper  limit  of  the  Haze  C  type  distribution:  A  Haze  C  model  with 
a  lower  limit  of  ,02y  (such  as  the  one  used  by  AFCRL;  and  a  Haze  x. 
with  lower  limit  ,03u,  the  lower  limit  used  by  Deirmendjian  (Ref  16: 
*«QS) .  Tiie  assumed  Indices  of  refraction  are  those  of  water  and  are 
given  in  Table  X  with  the  computed  attenuation  coefficients.  This 


computation  allowed  e  check  with  published  values  for  the  Deirmendjian 


Haze  C  {*03p-5p}  (R«f  18:183), 
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RELATIVE  ATTENUATION  COEFFICIENT  (o-km*1) 


ig.  10.  Effect  of  Changing  Power  Lav 
Slopes  of  Attenuation  Curves 

{Rof  i  (i  •  1 9  ^ 

v  —  •  •  •  4  •“/ 
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Table  X 

Sensitivity  Analyses  on  Effect 
of  Changing  Limits  of  Integration* 


Wavelength 

in 

Index  of  Refraction 

Attenuation  Coefficient  of 
Distribution 

Microns 

Real 

Imag 

.02  5u 

.03  +  5m 

1.34 

0 

.478E-04 

.52*»E-0*» 

.55 

1.34 

0 

.396E-04 

MkE-0k 

•  70 

1.33 

0 

.30^-0^ 

•333E-04 

1.0636 

1.33 

0 

.197E-04 

.217E-0A 

1.61 

1.315 

0 

.121E-04 

.132E-04 

.02  IOjj 

.03  10u 

1.34 

0 

MkE-0k 

.531 E-0^} 

.55 

1.34 

0 

.A02E-0*i 

ME~Qk 

.70 

1.33 

0 

.31 0E— 0^1 

,3k0E-0k 

1 .0636 

1.33 

0 

.203E-04 

.223l-0k 

1.61 

1.315 

0 

•  127E-0*» 

.139 E~04 

.02  20y 

.03  +  20y 

.*»5 

1.34 

0 

.ksn-ob 

.534E-04 

.55 

1.3*1 

0 

MkE-Ok 

.WiE-0*» 

.70 

1.33 

0 

.312E-0*» 

.343E-04 

1.0636 

1.33 

0 

,206E-0*i 

.226E-0I} 

1.61 

1.315 

0 

.  130E— 04 

J42E-04 

*Values  determined  for  Continental  Haze  Model 


As  the  upper  limit  is  increased  the  ratio  of  the  attenuation  coefficient 
at  .55m  to  the  one  at  1.Q636m  decreases.  The  addition  of  larger 
particles  to  the  distribution  increases  the  attenuation  more  at  the 
longer  wavelength. 
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Analysis  of  Refractive  Index 

The  two  distributions,  .03y  to  20y  and  .02y  to  20y,  computed 
with  the  refractive  index  of  water,  yield  a  comparison  with  the  same 
distributions  computed  with  a  different  refractive  index  assumed  in 
the  previous  section.  A  comparison  of  the  resulting  attenuation 
coefficients  demonstrates  the  fact  that  for  an  assumed  size  distri¬ 
bution  an  increase  in  the  refractive  index  results  in  a  decrease  in 
the  ratio  of  the  coefficient  at  .55 y  to  the  one  at  1.0636y,  a  fact 
explainable  by  Mie  Theory. 

Additional  sensitivity  studies  were  conducted  on  the  imaginary 
part  of  the  refractive  :ndex  holding  the  real  part  constant.  Results 
are  shown  in  Table  XI. 

Table  XI 

Effect  of  Changing  Imaginary  Part  of  the 
Refractive  Index  Holding  Real  Part  Constant 


Index  of 

Refraction 

Total  Attenuation 

Absorption 

Wavelength 

Real 

imag 

Coefficient 

Coefficient 

.55 

1.54 

.005 

.657E-04 

.034E-04 

.55 

1.5*1 

.011 

.658E-04 

.063E-04 

.55 

1.54 

.035 

.664E-04 

.141E-04 

.55 

1.5*1 

.048 

.667E-04 

.173E-04 

I.O636 

1-5** 

.017 

.356E-04 

.049E-04 

I.O636 

1.5** 

.022 

.358E-04 

.059E-04 

I.O636 

1.54 

.033 

.362E-04 

.079E-04 

1,0636 

1  -5*» 

.045 

.367E-04 

•098E-04 

1.0636 

1.5*i 

,068 

.377E-04 

.128E-04 

The  values 

selected  for 

the  imaginary 

part  were  those  reported  in 

references  66  and  34,  the  later  of  which  were  taken  from  measurements 
from  two  locations  in  Germany.  In  addition,  calculations  were  made 
for  mean  imaginary  indices  of  refraction.  The  results  show  that  the 
various  reported  values  for  the  Imaginary  part  do  not  significantly 
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affect  the  overall  computed  attenuation,  it  must  also  be  noted  that 
the  reported  values  for  this  wavelength  region  are  not  wildly  varying 
so  a  comparison  of  lowest  reported  values  at  one  wavelength  (.55 p) 
to  the  highest  reported  values  at  a-other  (1.0636y)  is  not  valid. 

A  more  important  analysis  was  that  conducted  using  the  refrac¬ 
tive  indices  reported  by  Wane]  in  Reference  3*1  for  measured  German 
aerosol  and  values  provided  by  Volz  for  water  soluble  aerosol.  The 
modified  Deirmendjian  Haze  C  size  distribution  was  used  (.03 ji  to  20ji) . 
Results  of  the  computation  show  that  later  reported  values  (Model  5, 
Model  6,  V0I2)  are  rather  consistent.  The  computed  ratios  of  the 
coefficient  at  .55p  to  the  coefficient  at  1.0636u  showed  a  trend.  The 
values  computed  from  the  indices  reported  at  Mainz,  Germany  Showed 
close  correlation  as  shown  in  Table  Xj  1.  The  computed  values  for 
indices  provided  by  Volz  and  those  from  Hanel 's  Model  6  showed  some 
similarity.  These  later  values  are  probably  from  a  more  natural  aerosol 
since  the  Model  6  site  was  elevated  and  away  from  urban  influence,  and 
the  values  p-ovlded  by  Volz  are  from  a  variety  of  sources  (Ref  6/}). 

The  difference  in  refractive  indices  reported  by  Volz,  Hanel, 
and  Fischer  (Ref  3iiS5)  may  well  be  due  to  th°  fact  that  European 
aerosols  are  modified  by  carbons  (soot)  and  fossi>]  Tv'S1  wastes  (Ref  60) . 
An  investigation  of  the  aerosol  absorption  coefficient  in  New  York  City 
reported  values  similar  to  those  reported  in  Germany  (Ref  12), 

Based  on  the  analysis  of  the  refractive  indices  previously 
described  and  comparison  of  various  reported  values,  this  author 
elected  to  use  an  average  of  the  indices  reported  for  Models  5  and  6 
(Ref  3*1:373)  as  the  aerosol  index  for  continental  Germany. 
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Table  Xf  I 

Per  Particle  Aerosol  Attenuation  Coefficients 
for  Continental  Haze  Model  (.Q3p  to  20p) 
for  Reported  Indices  of  Refraction 


Wave length 

Index  o^  Refraction* 

in 

Attenuation 

Model 

Microns 

Real 

Imag. 

Coefficient 

,4o 

1.64 

* 

.044 

.952E-04 

.55 

1 .62 

.051 

•737E-04 

Model  1 

•  70 

1 .61 

.058 

.595E-04 

(Mainz) 

1.00 

1.61 

.072 

.436E-04 

1.0636 

1.61 

.077 

•414e-o4 

.40 

1.57 

.042 

.883E-04 

.55 

1.55 

.048 

.676E-04 

Model  5 

.70 

1.54 

.055 

.544E-04 

(Mainz) 

1.00 

1.53 

.068 

.394E-04 

1 .0636 

1.55 

.072 

•374E-04 

.40 

1.53 

.013 

.842E-04 

Model  6 

.55 

1.51 

.015 

.630E-04 

(1000m  above  .70 

1.50 

.016 

.499E-04 

sea  level 

1  .00 

1.43 

.020 

.35  IE-04 

Germany) 

1.063& 

1,43 

.021 

•331E-04 

.40 

1.53 

.005 

.843E-04 

.55 

1.53 

.006 

.64/E-04 

Vol  2 

•70 

1*53 

.007 

•519E-04 

(Water 

i.oo 

1-52 

.016 

•366E-04 

Soluble) 

1.0636 

1,52 

.01? 

.345E-04 

,40 

1.55 

,022 

,  .864E-04 

Average  of 

.55 

1.53 

.031 

,654F.~y4 

Models  5 

.70 

1 .52 

.035 

•522E-04 

pnd  6 

1.00 

1.51 

.044 

.3732-04 

1.0636 

1.51 

,046 

,3536-04 

^Values  of 

refractive 

indices  are  those 

used  in 

Table  VI. 
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These  reported  values  were  extrapolated  to  1.0636p  by  this  author. 

When  these  average  values  are  used  in  the  modified  Haze  C  distribution, 
the  resulting  attenuation  coefficients  closely  approximate  those  com¬ 
puted  using  the  Volz  indices: 

Aerosol  Model 

This  computation  model  assumes  a  size  distribution  of  the  form 
of  the  Deirmendjian  h’aze  C  with  the  exception  of  the  upper  limit  which 
Is  extended  to  20p.  The  index  of  refraction  varies  with  wavelength 
but  is  an  average  of  the  values  reported  by  Hanel .  Resulting  Mie  Theory 
computations  shows  the  aerosol  attenuation  coefficient  at  1.0636p  to  be 
related  to  the  aerosol  attenuation  coefficient  at  ,55y  by 

/ 

cr  (1 .0636p)  =  .5394  a  (.55*0  (41) 

O  3 

where  the  constant  .5394  is  determined  uniquely  by  Mie  Theory  for  an 
assumed  distribution  and  the  average  indices  of  refraction  shown  in 
Table  XII.  Constants  for  other  distributions  are  given  in  Table  XIII. 

Effect  o_f  Relative  Humidity  on  Model 

The  effect  of  increasing  relative  humidity  on  aerosol  size 
distributions  and  index  of  refraction  was  discussed  in  Chapter  5*  This 
effect  of  particle  growth  (and  shift  in  size  distribution)  and  change 
in  index  of  refraction  has  been  investigated  by  several  workers.  These 
investigations  may  be  classified  as  either  computer  simulation  programs 
or  as  experimental  determinations. 

Notable  among  the  former  investigations  are  the  works  of 
Ba.rnhardt  and  Streete  (Ref  2)  and  John  A.  Hodges  (Ref  38), 
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Table  XII J 

Ratios  of  1.0636m  Attenuation  Coefficient  to  .5 5m  Attenuation 
Coefficient  for  Several  Aerosol  Distributions  with 
Same  Assumed  Refractive  Indices* 


Aerosol 

Distribution 

a(l .0636y) 

o'(  .55u) 

Haze  L 

.796 

Power  Law 

5T4(.1m  ->  20y) 

.538 

Maritime 
(.003u  ■*  20m) 

1.00*» 

25%  Maritime 
+  75%  Continental 

.932 

50%  Maritime 
+  50%  Continental 

.977 

^Refractive  Indices  used  are  the  Average  Model  Indices  from 
Table  XII. 
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Both  provide  a  growth  factor  computation  based  on  the  fractional  growth 
of  a  particle  developed  by  Wright  in  1939.  Using  the  growth  factor,  a 
new  distribution  is  computed  based  on  the  relative  humidity  assum¬ 
ptions  of  the  physical  rate  of  growth.  Computations  of  one  growth 
factor  for  the  entire  distribution  also  assume'-  that  all  particles  grow 
at  the  same  rate.  In  addition,  these  programs  compute  a  real  index  of 
refraction,  which  is  an  analytical  fit  to  the-  data  of  Bull  rich  (Ref  2: 
1339). 

This  method  was  rejected  for  use  in  this  report  for  three 
reasons:  (l)  the  technique  involved  assumptions  that  were  not  neces¬ 
sarily  applicable  to  Germany  (indices  of  refraction);  (2)  the  results  of 
the  technique  have  not  yet  been  compared  to  measured  data;  and  (3)  the 
technique  involved  lengthy  computer  programming.  Instead,  the  results 
of  aerosol  measurements  were  investigated. 

Only  one  of  the  reports  investigated  offered  the  relationship 
of  particle  growth  at  various  wavelengths.  This  information  is 
imperative  for  use  in  the  aerosol  model  since  the  model  relates 
attenuation  at  1.0636ii  to  that  at  .55u.  This  investigation  (Ref  7) 
reported  several  significant  results.  First,  it  was  found  that  the 
growth  rate  of  all  particle  sizes  was  approximately  the  same  up  to 
a  Relative  Humidity  of  95%.  Beyond  $5%  the  large  particles  grow  at  a 
faster  rate  skewing  the  distribution.  In  addition,  the  aerosol 
attenuation  coefficient  was  measured  for  various  relative  humidities 
as  shown  in  Figure  11.  The  wavelength  dependence  of  extinction  is  shown 
at  the  left  in  Figure  11  for  different  relative  humidities.  The  findings 
of  this  report  are  largely  substantiated  by  similar  measurements  by 
Winkler  and  Junge  (Ref  70). 
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The  increase  in  relative  humidity  does  have  some  effect,  however, 
on  attenuation  as  shown  in  the  figure.  The  ratio  of  these  coefficients 
is  given  as  a  function  of  the  relative  humidity.  Then  as  the  relative 
humidity  increases  there  will  be  an  increase  in  the  attenuation  coef¬ 
ficient  which  is  equal  for  wavelengths  .55u  and  }. 0636*1.  The  important 
concept  at  this  point  is  that  the  ratio  between  coefficients  at  the  two 
wavelengths  stays  the  same. 

Aerosol  Attenuation  from  Meteorological  Observables 

Sophisticated  equipment  for  measuring  the  attenuation  of 
radiation  due  to  aerosols  is  not  widely  available  nor  nractical.  Be¬ 
cause  of  this  fact,  it  becomes  necessary  to  relate  some  atmc  iheric 
quantity  to  the  aerosol  attenuation.  Such  a  relationship  for  daylight 
was  first  expressed  by  Koschmieder  (Ref  *»3)  and  later  repeated  more 
elaborately  by  Middleton  (Ref  53)-  Koschmieder 's  lav/  is  expressed  by 

m  d  _  3.912  (42) 

fl  •  K.  —  * 

a 

a 

where  M.R.  =  Meteorological  Range  (Km)  and  a  =  aerosol  attenuation 

d 

coefficient  (Km-1). 

Meteorological  range  in  white  light  is  the  maximum  distance 
at  which  an  observer  can  barely  detect  a  large  dark  object  against  a 
white  background.  It  is  the  distance  at  which  the  object  to  background 
contrast  (inherent  contrast)  is  reduced  to  .02  of  its  original  value 
at  the  eye  of  the  observer  (apparent  contrast).  In  addition,  the  large 
dark  object  must  subtend  so  large  an  angle  at  tlu  eye  of  the  observer 
that  if  the  subtended  angle  were  greater,  the  reported  value  of 


5) 


GF.P/PH/74-4 


meteorological  range  would  not  be  changed  (Ref  23:186). 

It  must  be  emphasized  that  meteorological  range  is  not  visi¬ 
bility  as  reported  by  the  meteorologist.  The.  meteorologist  reports 
prevailing  visibility,  which  is  the  greatest  horizontal  visibility 
prevailing  throughout  at  least  half  of  the  horizon  circle  wich  may 
not  necessarily  be  continuous  (Ref  29:A-6~3).  1°  addition,  distant 

dark  objects  of  sufficient  angular  size  are  rarely  available  to  the 
meteorologist  (Ref  23:186).  These  objects  when  available  are  usually 
limited  in  number  and  not  randomly  distributed  so  as  to  reduce  the 
statistical  probability  of  selecting  the  exact  meteorological  range 
(Ref  52:40).  Because  of  these  factors  the  reported  prevailing 
visibility  (v)  is  usually  less  than  the  actual  meteorological  range 
(MR). 

A  comparison  of  reported  prevailing  visibilities  'with  measured 
meteorological  ranges  by  the  Tiffany  Foundation  and  the  National 
Bureau  of  Standards  concludes  that,  on  the  average,  the  reported  pre¬ 
vailing  visibility  was  75%  of  the  true  meteorological  range  (Ref  23:186). 

Equation  (42)  now  becomes 


V 


2.934 


(43) 


where  V  -  Reported  visibility,  prevailing  visibility,  or  "yisjbi 1 ity." 

For  determination  of  the  aerosol  attenuation  this  relationship  is 
assumed  to  be  valid  at  the  discrete  wavelength  «55u  which  is  the  approxi¬ 
mate  center  of  eye  sensitivity  for  photopic  vision  (ordinary  color 
vision  in  daylight). 
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Visibility  at  night  is  a  more  complex  phenomena.  In  the  ab¬ 
sence  of  a  lighted  object,  the  eye  detects  contrast  of  the  object  but 
the  wavelength  sensitive  region  of  the  eye  has  shifted  to  approximately 
•53p  (scotopic  or  rod  visition).  The  threshold  detectable  contrast  at 
the  eye  for  night  vision  is  complicated  by  a  chemical  process  commonly 
called  dark  adaptation  which  is  the  process  by  which  the  eye  adjusts 
to  lower  levels  of  illumination.  Dark  adaptation  for  the  normal  eye 
(scotopic  region)  takes  approximately  30  minutes,  although  dark 
adaptation  time  for  the  photopic  region  generally  occurs  in  approximately 
8  minutes  (Ref  5^:39). 

Because  of  dark  adaptation,  meteorologists  generally  use  light 
sources  as  meteorological  range  markers.  The  meteorological  range  at 
night  as  determined  by  observing  a  point  source  of  light  is  given  by 


NMR 


-  2  £n  (NMR) 


where  NMR  =  Night  Meteorological  Range, 

Et  =  Threshold  flux  density  of  the  eye, 
and  i  =  Intensity  of  light  source 


m 


This  transcendental  equation  known  as  Allard's  Law  is  best  solved 
graphically.  The  proper  value  of  £t  depends  on  the  luminance  of  the 
background  and  on  the  probability  of  detection  (Ref  53:138).  A 
graphical  solution  to  this  equation  is  presented  in  Figure  12.  This 
solution  Is  based  on  a  95%  detection  probability  and  assumes  no  other 
light  within  the  field  of  view. 


93 


GEP/PH/7M 


As  is  evident  from  Figure  12,  the  changes  in  threshold  with 
variations  in  surrounding  brightness  are  much  more  pronounced  at  night 
than  daytime  (using  1-103  footlamberts  for  daytime).  As  we  have  also 
seen,  the  state  of  dark  adaptation  of  the  observer's  eye  is  an  import¬ 
ant  factor.  For  this  reason,  the  adoption  of  s  standard  threshold  of 
flux  density  at  night  has  little  practical  value  (Ref  5H:56).  In 
addition,  meteorological  observers  are  rarely,  if  ever,  dark  adapted 
for  night  visibility  observations  so  that  these  observations  must  be 
considered  quite  unreliable  and  useless  for  comparative  studies  of 
visibility  (Ref  54:40). 

Because  night  visibilities  as  reported  are  considered  less 
reliable  than  day  visibilities  and  because  there  are  additional  compli¬ 
cating  factors  in  measuring  visibility  at  night,  night  observations 
were  not  used  in  determining  the  aerosol  attenuation  coefficient.  The 
computations  within  this  report  are  limited,  therefore,  to  daylight 
hours. 

Use  of  the  Model 

As  we  have  discussed  previously,  the  total  particle  density 
is  essentially  unmeasureable.  For  this  reason,  we  have  limited  our 
discussion  to  the  per  particle  attenuation.  The  per  particle 
attenuation  is  useless  unless  there  is  a  method  of  measuring  the  number 
of  particles.  This  number  is  determined  by  the  visibility. 

in  the  last  section  we  discussed  a  method  (based  on  visual 
observations)  to  determine  the  attenuation  coefficient  at  . 55v .  Given 
the  reported  visibility,  Eq  (43)  determines  the  coefficient  at  .55v. 

If  one  divides  this  attenuation  coefficient  by  the  one  determined  by 
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Hie  calculations  at  .55d  then  the  number  of  particles  may  be  deter¬ 
mined.  This  exercise  is  not  necessary,  however,  since  the  relationship 
between  attenuation  coefficients  at  the  two  wavelengths  (-S5y  and 
1.0636y)  has  been  established  by  Eq  (41*1.  Given  the  reported  visibility 
then  one  may  determine  directly  t,.  uation  coefficient  at  1 .0636P 

by  combining  Eqs  (41)  and  (43) 

a  (1.0636m)  =  (45) 

3  V 

where  oa(1.0636y)  =  aerosol  attenuation  coefficient  in  Km-1  for  a 
horizontal  path. 

Other  Model;; 

During  the  course  of  this  study  several  models  for  determination 
of  aerosol  attenuation  coefficients  were  investigated.  These  included: 

(1)  the  AFCRL  Report  Optical  Properties  of  the  Atmosphere  (Ref  49); 

(2)  1964  AFCRL  Report  by  Elterman  (Ref  26);  (3)  Handbook  of  Geophysics 
and  Space  Environments  (Ref  61);  (4)  Electro  Optics  Handbook  (Ref  24); 
(5)  1968  AFCRL  Report  by  Elterman  (Ref  27);  and  (6)  1970  AFCRL  Report 
by  Elterman  (Ref  28).  Briefly,  this  investigation  focused  on  differ¬ 
ences  in  reported  attenuation  coefficients,  in  basic  model  assumptions, 
and  on  differences  with  the  model  derived  in  this  report. 

The  1964  and  1 968  AFCRL  reports  by  Elterman,  the  Electro-Optics 
Handbook,  and  the  Handbook  of  Geophysics  and  Space  Environments,  are 
all  based  on  measurements  made  by  Curcio  (Ref  14)  in  the  NRL  project 
discussed  earlier.  The  attenuation  coefficients  given  in  these  reports 
are  based  on  the  measurements  reported  by  Curcio  for  visibility  of  25 
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km  (although  the  Electro-Optics  Handbook  reports  this  visibility  as 
23.5  km).  The  Electro-Optics  Handbook  then  apparently  extrapolates 
these  values  to  other  meteorological  ranges  using  Eq  (42)  as  shewn  in 
Figure  13.  Curves  of  the  original  values  are  shown  in  Figure  8  and 
again  grouped  in  Figure  14  (Ref  14:12).  The  curvature  of  the  lines  of 
equal  visual  range  in  the  NRL  reports  is  not  the  same;  in  fact,  an 
analysis  of  the  data  shows  only  Meteorological  Ranges  between  6  and  12 
miles  and  approximately  2 5  miles  give  identical  slopes.  In  addition, 
it  is  recalled  that  these  data  were  taken  in  an  area  where  the  aerosol 
attenuation  is  influenced  to  some  degree  by  a  maritime  distribution 
component.  In  fact,  the  data  indicates  a  general  correlation  between 
meteorological  range  and  slope  of  the  particle  size  distribution  curve 
such  that  a  change  from  high  to  lew  meteorological  range  generally 
indicates  a  steeper  slope  to  the  particle  size  distribution  curve 
(Ref  14:12).  Because  of  these  factors,  these  reports  may  not  offer  a 
method  for  determination  of  the  aerosol  attenuation  coefficient  for 
continental  geographical  locations. 

The  1970  AFCRL  Report  by  Elterman  provides  a  method  for 
determining  aerosol  attenuation  coefficients  for  eight  surface  meteor¬ 
ological  ranges  from  2  to  13  km.  This  report  has  utilized  the 
measurements  taken  by  Curcio  for  4  km  and  extrapolated  from  2  to  13  km 
by  using  this  relation  from  Reference  28: 


(46) 


where  o  (V  ,,\)  =  aerosol  attenuation  coefficient  the  visual 

51  n 

range,  q,  and  the  wavelength,  A 
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ca(V4 ,X)  =  the  measured  aerosol  attenuation  coefficient  at  visual 

range  of  4  km  and  wavelength,  A,  and 

3  (A  rr)  =  Rayleigh  (molecular)  scattering  coefficient.  The  use  of 
r  .55 

this  relation  assumes  a  constant  particle  size  distribution  for  this 
meteorological  range  which  apparently  is  not  supported  by  Curcio's 
measurements  (Ref  14:11,12).  Elterman's  report  also  has  the  same 
disadvantage  for  use  at  continental  geographical  locations  as  discussed 
previously. 

The  1972  AFCRL  Report  titled  Optical  Properties  of  the  Atmos¬ 
phere  uses  a  Continental  Haze  model  after  the  one  originally  used  by 
Ds i rmendj ian  (Ref  49:8).  The  distribution  used  in  this  CRL  report 
differs  from  the  one  developed  here  in  upper  and  lower  limits  of 
Integration;  the  CRL  distribution  is  defined  from  .02y  to  lOp  (Ref  49:8), 
In  addition,  the  index  of  refraction  used  is  1,50  for  X  =  ,55u  and 
1.50  -  ,033  for  A  -  1.0636p.  These  differences  cause  a  reduction 
in  the  ratio  of  the  attenuation  coefficients  at  1 .0636p  to  ,55p.  In 
other  words,  this  AFCRL  report  tends  to  be  slightly  more  optimistic 
than  the  model  developed  here-  This  AFCRL  model  was  discussed  previously 
in  the  size  distribution  and  index  of  refraction  analysis  section. 

Model  Lim?  tations/Assumptions 

The  model  developed  in  this  report  for  determination  of  the 
aerosol  attenuation  coefficient  for  a  horizontal  path  at  ground 
(observation)  level  is  restricted  for  use  in  areas  where  the  particle 
size  distribution  is  described  approximately  by  the  assumed  continental 
haze  model.  It  is  also  restricted  for  use  in  areas  where  the  Index 
of  refraction  is  approximately  the  one  used  here,  although  Index  of 
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refraction  is  not  as  sensitive  a  factor  in  establishing  the  relation¬ 
ship  between  wavelengths  as  the  size  distribution. 

This  model  also  assumes  that  the  reported  visibilities  are  75% 
of  the  true  meteorological  range.  The  concept  of  meteorological  range 
further  assumes  that  Rayleigh  scattering  and  molecular  absorption  are 
negligible.  This  assumption  is  valid  for  meteorological  rano.es  on  the 
order  of  those  used  in  this  analysis  for  tactical  application.  The 
model  is  further  limited  to  relative  humidities  at  or  below  95%. 
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VII.  Ae rosol  Attenuation  for  a_  Slant  Path 

Chapter  VI  determined  the  aerosol  attenuation  coefficient  at 
ground  level  based  on  surface  observations.  This  chapter  discusses 
the  distribution  of  aerosols  with  respect  to  altitude  and  relates 
this  distribution  to  atmospheric  phenomena.  It  then  discusses  the 
feasibility  of  determining  aerosol  attenuation  coefficients  from 
ground-based  observations. 


Particle-Size  Pi stribution  (Function  of  A1 t i tude) 

We  have  discussed  briefly  the  particle  size  distribution  as  a 
function  of  altitude  in  Chapter  V.  The  particle  size  distribution  is 
found  to  remain  relatively  constant  with  altitude,  with  the  exception 
of  the  giant  particles  which  are  more  affected  by  sedimentation  (Ref 
41:182).  There  is  a  decrease  ?n  the  percentage  concentration  of  giant 
particles  with  increasing  altitude  in  the  absence  of  convective  cur¬ 
rents  that  force  them  aloft.  These  larger  particles  settle  into  more 
dense  air  where  the  increased  Brownian  motion  will  keep  the  particles 
ai rborne. 


There  is  evidence  to  support  the  existence  of  a  world-wide 
background  aerosol  distribution.  This  aerosol  distribution  is  the 
aerosol  distribution  for  the  troposphere  outside  the  contamination 
layers.  These  contamination  layers  are  the  layers  closest  to  the 
earth's  surface,  and  therefore  they  are  affected  by  the  local  production 


mechanisms.  This  background  aerosol  distribution  is  represented  as  a 


modified  continental  aerosol  distribution  (Ref  4i:i8y).  That  is,  the 


distribution  differs  significantly  from  the  continental  distribution 


only  In  the  relative  concentration  of  larger  particles. 
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There  is  some  disagreement;  also,  as  to  the  exact  lowei 
particle  limit  just  as  in  the  case  of  the  continental  distribution. 

Some  contend  that  the  lower  limit  should  be  somewhat  higher,  and  the 
concentration  of  the  smaller  particles  (radius  <  .ly)  be  less  than 
the  continental  distribution  (Ref  40:201).  Recent  laboratory  studies, 
however,  have  established  that  small  particles  may  be  produced  in  the 
atmosphere  by  sunlight  ionization  of  trace  gases  (Ref  40:202).  The 
lower  limit  is  not  important,  however,  as  long  as  this  limit  is  less 
than  .ly  since  particles  smaller  than  this  limit  are  ineffective 
scatterers  of  visible  and  infrared  radiation. 

The  background  aerosol  distribution,  then,  differs  from  the 
continental  aerosol  distribution  in  the  concentration  of  giant  particles, 
particles  whose  radii  are  larger  than  ly  (Ref  41:127).  Because  the 
measurement  of  aerosols  at  altitude  presents  formidable  problems  and 
because  of  the  expense  of  such  measurements,  the  change  in  the  relative 
concentration  of  these  larger  particles  with  altitude  is  not  v/el 1 
documented. 


Total  Particle  Concentration  (Function  of  Alti tude) 

In  contrast,  the  decrease  in  the  concentration  of  measureable 
aerosols  is  well  documented.  The  results  of  early  investigations  show 
that  the  average  decrease  in  particle  number  concentration  is  approxi- 
materly  exponential  with  altitude  (Ref  41:184).  Some  of  these 
investigations  show  an  exponential  decrease  to  a  certain  altitude; 
above  this  altitude,  the  number  concentration  decreases  less  rapidly 
or  is  constant  as  shown  in  Figure  15. 
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extensive  vertical  aerosol  study  conducted  by  AFCRL  in  1 966  (Ref  25). 
This  I966  study  included  a  total  of  119  profiles  of  the  troposphere 
(and  stratosphere).  These  profiles  taken  over  New  Mexico  were  com¬ 
prised  of  absolute  values  of  aerosol  attenuation  coefficients.  A 
comparison  of  the  two  resulting  models  (Elterman,  1964,  and  Elterman, 
1968)  is  shown  in  Figure  36. 

Mixing  Depth  Concept 

Individual  attenuation  measurements  and  average  profiles  have 
led  to  the  mixing  layer  concept  in  an  attempt  to  physically  explain 
the  vertical  aerosol  distribution.  As  we  discussed  earlier  in  this 
chapter,  the  average  aerosol  concentration  above  a  certain  altitude 
is  either  constant  or  gradually  decreasing.  Below  this  altitude, 
the  time-average  aerosol  concentration  profile  shows  a  rather  steep 
exponential  decrease.  This  lower  layer  described  by  the  steep  ex¬ 
ponential  profile  is  called  the  mixing  layer  (although  this  is  not  the 
strict  meaning  of  the  term). 

Basically,  the  mixing  layer  is  the  vertical  expanse  above 
the  surface  that  is  affected  by  the  aerosols  being  produced  at  the 
surface.  This  boundary  layer  as  it  is  also  called  is  the  layer  in 
which  vigorous  vertical  and  horizontal  m  ixing  of  the  air  takes  place 
(Ref  55:1).  Thus,  the  boundary  or  mixing  layer  is  largely  affected 
by  the  aerosol  production  mechanisms.  Above  this  layer,  the  aerosol 
concentration  is  related  to  the  air  mass  itself. 

The  depth  of  the  mixing  layer  (or  the  maximum  altitude  of  the 
boundary  layer)  has  been  the  subject  of  several  investigations. 


ALTITUDE  (km) 
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The  studies  of  Siedentopf  (1944) ,  Penndorf  (1954),  Faraponova  (1965), 
Rosen  (1967),  and  others  reported  in  Reference  58,  show  that  the 
average  boundary  layer  depth  is  approximately  5  km,  although  the 
reported  average  height  varies  from  3  km  to  6.5  km  (Ref  28:7,5).  Such 
large  values  of  the  mixing  depth  do  not  agree  with  meteorological 
computations  of  the  mixing  depth  based  on  convective  and  advective 
currents  some  of  which  are  discussed  briefly  later. 

AFCRL  Verti cal  Aerosol  Model 

Based  on  the  evidence  available  at  the  time,  Dr.  Elterman 
devised  an  extensive  vertical  aerosol  model.  This  1970  AFCRL  model 
(Ref  28)  assumes  the  mixing  depth  to  be  5  km.  The  aerosol  attenuation 
coefficients  above  5  km  are  those  reported  in  the  1968  vertical  aerosol 
model  revision  (Ref  27).  Below  5  km  the  aerosol  attenuation  coef¬ 
ficient  exponential ly  decreases  with  altitude  up  to  5  km.  Eight 
surface  meteorological  ranges  are  used  and  these  surface  meteorological 
ranges  are  based  on  the  aerosol  attenuation  coefficients  and  ranges 
reported  by  Curcio  (Ref  14)  as  we  discussed  earlier  in  Chapter  V!. 

Each  of  these  eight  surface  attenuation  coefficients  is  exponentially 
decreased  to  agree  with  the  same  reported  5  km  value.  Hence,  different 
exponential  scaling  factors  are  used  as  shown  in  Figure  17. 

The  most  recent  AFCRL  model  of  the  aerosol  attenuation  coef¬ 
ficient  (Ref  49)  also  uses  the  exponential  scaling  factor  concept. 

This  model,  however,  uses  different  attenuation  coefficients  at  ground 
level — these  based  on  computations  using  a  Deirmendjian  Haze  C  model 
as  we  discussed  earlier.  The  aerosol  attenuation  coefficients  above 
5  km  are  the  same  as  those  in  the  1968  AFCRL  Elterman  Report  (Ref  27). 
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Fig.  17.  Family  of  Scaling  Factors  Used  by 
Elterman  to  Relate  Different  Ground  Level 
Attenuation  Coefficients  to  5  km 
Attenuation  Coefficient.  The  Scale  Heights, 
Hp,  are  Computed  at  Reference  Frequency  .55 y 
(Ref  28:9). 


Latest  Attenuation  Measurements 

Recent  aircraft  aerosol  attenuation  measurements  made  for 

AFCRL  (Ref  20  and  21)  do  not  support  the  model  developed  by  Elterman 

(Ref  30  and  Ref  67).  These  extensive  total  volume  scattering 

measurements  were  made  in  Mew  Mexico  under  project  name  Atom  and  in 

Southern  Germany  under  project  name  Haven  View.  These  measurements 

show  that  the  mixing  layer  is  somewhat  less  than  the  5  km  computed 
/  . 

from  averages  of  many  measurements.  In  fact,  the  measurements  show 
that  the  mixing  layer  extends  from  ground  level  to  altitudes  from 
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approximately  five  hundred  meters  up  to  around  2.5  km  with  levels  of 
1.5  km  predominating  (Ref  30). 

These  measurements  also  show  that  in  the  mixing  layer  the  at¬ 
tenuation  coefficient  shows  a  very  slow  exponential  decrease  (on  the 
order  of  7.99  scale  height  used  with  the  standard  atmosphere  for 
molecular  computations)  or  it  is  approximately  constant.  At  the  mixing 
layer  boundary,  the  aerosol  attenuation  coefficient  shows  a  sharp 
decrease  to  a  clean  air  mass  attenuation  value  or  a  gradual  transition 
to  such  a  clean  air  mass  value  through  a  200  meter  layer.  Above  the 
boundary  layer,  the  attenuation  coefficient  again  shows  a  gradual 
exponential  decrease  but  this  decrease  is  frequently  less  than  the 
density  lapse  rate  (i.e.,  7"8  km  scale  height)  (Ref  22 : 7~1 7) * 

These  whole-volume  attenuation  measurements  described  in 
References  20,  21  and  22  are  well  documented,  and  extensive  sets  of 
meteorological  data  were  recorded  in  the  course  of  these  studies. 
Consequently,  various  attenuation  coefficients  are  related  to  specific 
meteorological  processes.  Measurements  of  high  values  for  aerosol 
attenuation  in  the  mixing  layer  are  associated  with  temperature 
inversions  or  layers  of  high  stability.  The  aerosol  attenuation  coef¬ 
ficient  is  relatively  constant  throughout  the  mixing  layer,  then  at 
the  mixing  height  there  is  a  sharp  (step  function)  decrease  to  the  air 
mass  attenuation  coefficient  above  the  layer  as  shown  in  Figure  18. 

When  the  mixing  layer  attenuation  coefficient  is  low  (that  is, 
the  visibility  is  good  indicating  very  clean  air)  then  measurements 
show  that  there  is  a  slow  exponential  decrease  in  the  attenuation 
coefficient  with  increasing  altitude. 
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Fig.  18.  Vertical  Temp  "-9*ure  and  Whole  Volume  Scattering  Coefficient 
Distributions  Measured  by  AFCRL  Over  Germany.  The  top  set  is 
from  Might  C-137  and  the  lov/er  set  is  from  Flight  C-139.  Filter 
Number  5  approximates  the  photopic  (.557v). 
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This  slow  exponential  decrease  Is  evident  even  in  cases  where 
a  strong  inversion  exists.  The  inversion  then  has  relatively  little 
effect  on  clean  ai r. 

In  addition,  there  is  some  variance  in  the  value  of  the  aerosol 
attenuation  coefficient  reported  for  the  air  mass  above  the  mixing 
layer  as  shown  in  Figure  18.  This  indicates  that  the  air  mass  above 
the  mixing  layer  may  be  affected  to  some  degree  by  the  aerosol  concen¬ 
tration  of  the  boundary  layer.  Some  mechanical  or  convective  transport 
must  occur  (Ref  30).  However,  generally  the  air  mass  above  the  mixing 
layer  approximates  a  clear  day  or  better. 

Average  Aerosol  Prof i les 

Individual  vertical  aerosol  profiles  as  discussed  in  the  previous 
section  show  certain  characteristics  which  are  lost  when  average  values 
are  taken.  The  fact  that  these  average  profiles  tend  to  show  some 
exponential  decrease  is  probably  due  to  a  correlation  between  a  changing 
height  of  the  mixing  layer  and  the  concentration  of  the  particies 
within  the  layer  (Ref  41:183,184).  in  addition,  the  fact  that  earlier 
measurements  of  the  vertical  aerosol  distribution  yielded  higher  mixing 
height  values  is  possibly  due  to  the  fact  that  aircraft  flights,  and 
particularly  those  of  several  years  ago,  usually  occur  in  good  weather 
flying  conditions.  This  higher  average  mixing  layer  height  may  also 
be  due  to  the  changing  individual  mixing  layer  height. 

The  difficulty  in  using  such  average  vertical  profiles  may 
best  be  Illustrated  by  the  following  reducio  ad  absurdum.  Suppose 
that  in  the  course  of  one  year,  a  station  recorded  6  months  of  fairly 
good  weather.  That  is,  the  visibility  is  excellent  and  the  aerosol 
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attenuation  coefficients  are  low  vaiues.  Further,  suppose  that  during 
the  other  6  months  the  visibilities  are  extremely  low  (on  the  order 
of  2  km)  and  aerosoi  attenuation  coefficients  are  high,,  if  the  missile 
system  being  evaluated  is  effective  only  in  clear  air,  then  an  average 
profile  might  unjustifiably  lead  to  the  conclusion  that  this  missile 
would  be  ineffective  at  this  location  when  in  fact  it  might  be  quite 
effective  50%  of  the  time. 

Vertical  Aerosoi  Model 

For  reasons  we  have  discussed  in  the  previous  section,  it  was 
decided  that,  a  more  representative  vertical  aerosol  model  would  be  a 
model  that  accurately  represented  individual  characteristic  aerosol 
profiles.  These  profiles  are  based  on  the  published  results  of 
Project  Haven  View  and  Project  Atom  (Ref  20,  21,  and  22)  and  discussions 
with  the  project  manager  Dr.  Robert  Fenn  (Ref  30). 

These  measurements  demonstrated  that  there  is  probably  a  mixing 
or  boundary  layer  in  which  the  particles  are  fairly  uniformly  mixed 
and  that  above  the  layer  the  air  is  clear.  In  fact,  measurements  of 
the  aerosol  attenuation  coefficient  in  the  air  mass  above  the  boundary 
layer  showed  the  coefficient  to  be  characteristic  of  air  in  which  the 
meteorological  range  is  b5  km  or  more.  The  individual  measurements 
gave  ranges  from  approximately  k5  km  to  approximately  100  km  (Ref  20 
and  21),  although  there  may  be  some  error  in  the  exact  value  due  to 
difficulty  in  calibrating  equipment  in  very  clear  air.  That  is,  the 
Rayleigh  scattering  component  is  no  longer  negligible  in  the  visible 
region  at  such  meteorological  ranges.  In  fact,  the  Rayleigh  scattering 
component  at  aV nwteoro logical  range  of  bO  km  is  12%  and  it  increases 
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as  the  meteorological  range  increases.  These  measurements  were  taken 
•with  a  filter  centered  at  .557u  (an  approximation  of  photopic  vision), 

A  clear  air  model  is  thus  representative  of  the  air  mass  above 
the  mixing  or  boundary  layer.  An  air  mass  whose  meteorological  range 
Is  approximately  equal  to  40  km  (prevailing  visibility  of  30  km)  was 
chosen  for  the  air  '•  ss  model  because  this  value  represented  a  worst 
case  among  the  reported  measurements. 

The  mixing  or  boundary  layer  is  assumed  to  be  uniformly  mixed 
vertically  and  horizontally.  Its  depth  is  dependent  on  a  number  of 
factors.  These  factors  are  the  thermal  and  mechanical  forces  of  the 
atmosphere.  The  measurements  conducted  by  AFCRL  in  Southern  Germany 
indicate  the  mixing  depth  extends  to  1.5  km  during  late  May  and  early 
June  when  the  measurements  were  taken.  A  recent  study  by  USAF  ETAC 
(Ref  55)  based  on  theoretical  meteorological  processes  (thermal  and 
mechanical)  for  one  location  in  Germany  indicates  that  during  this 
period  (May  to  June)  the  mixing  depth  varies  from  150-1175  km.  The 
wide  range  of  reported  values  is  due  to  the  fact  that  the  ETAC  study 
covered  the  2k  hour  period  while  aircraft  measurements  were  made  during 
daylight  hours.  However,  the  maximum  theoretical  computations  do  not 
predict  the  mixing  height  measured  by  AFCRL. 

This  same  USAF  ETAC  study  for  Stuttgart,  Germany,  indicates 
that  the  mixing  depth  extends  to  lower  altitudes  in  winter  than  during 
other  times  of  the  year  (see  Figure  19).  This  corresponds  to  what 
others  have  observed  to  be  true:  Winter  is  the  low-level  haze  season 
in  Germany  (Ref  30) .  There  have  been  too  few  aircraft  measurements; 
however,  to  substantiate  this  seasonal  variation. 
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Moreover,  the  theoretical  model  does  not  show  significant  variation  in 
the  mean  seasonal  mixing  depths  (Fig  19). 

Boundary  (Mixing)  Layer  Depth 

There  are  many  methods  for  calculation  of  the  mixing  depth. 

Some  are  based  on  aircraft/balloon  measurements;  others  are  theoretical 
In  nature.  Still  others  combine  both.  The  measurements  of  Elterman 
and  AFCRL  have  been  discussed;  the  measurements  of  others  (Siedentopf, 
etc.)  have  also  been  mentioned.  This  section  discusses  two  theoretical 
models  and  discusses  the  meteorological  processes  involved. 

One  of  the  theoretical  methods  of  determining  mixing  depth 
height  is  the  method  of  Miller  and  Holzworth  (Ref  39)*  This  method 
is  currently  being  used  by  the  Environmental  Protection  Agency  (EPA) 
and  the  National  Weather  Service  (NWS)  together  with  radiosonde  data 
to  predict  the  thickness  of  the  mixing  layer.  This  method  is  based 
on  the  maximum  and  minimum  surface  temperatures  and  the  radiosonde 
data.  The  maximum  mixing  depth  height  is  the  height  to  v/hich  a  parcel 
of  ground  air  (at  maximum  temperature)  will  ascend  if  lifted  dry 
adiabatical ly  until  its  temperature  is  in  equi 1 ibrium  wi th  the  environ¬ 
ment  (Ref  5 1 3)  •  This  method,  however,  fails  to  account  for  moisture 
In  the  air  and  it  is  restricted  to  occasions  when  there  is  no 
precipitation  and  no  air  mass  change.  It  is  also  restricted  for  use  to 
locations  where  radiosonde  balloon  observations  are  made. 

The  method  proposed  by  USAF  ETAC  (Ref  55)  uses  weighted 
contributions  of  thermal  and  mechanical  parts.  The  weighting  factors 
are  Pasquill  stability  factors  which  are  based  on  meteorological 
conditions  as  shown  in  Figure  20. 
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categories  are  used  to  weight  the  thermal  and  convective  parts  of  theoretical 
mixing  height.  Category  A  contributes  all  its  weight  to  the  thermal  part  and 
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The  mixing  height  is  then  calculated  based  on  ground  observations 
and  estimates  of  upper  level  winds. 

Observations  using  these  and  other  methods  show  that  there  is 
great  variability  of  the  mixing  depth  during  the  day  and  also  from  day 
to  day.  On  the  other  hand,  aircraft  measurements  do  not  show  such 
wide  variations  in  the  mixing  height  on  an  hour  to  hour  basis.  This 
may  best  be  explained  by  considering  the  history  of  the  air  mass. 

If  the  mixing  occurs  in  the  lower  levels  of  a  clean  air  mass,  then  the 
aerosols  are  distributed  to  the  maximum  mixing  layer  height.  Later 
in  the  day  (evening  and  night)  when  the  convective  mixing  processes 
diminish,  the  aerosols  remain  aloft  if  there  are  no  convective  cur¬ 
rents  forcing  them  downward.  The  particles  are  influenced  by  gravity 
but  the  settling  rate  is  slow.  Thus,  on  subsequent  days,  the  aerosol 
remain  mixed  to  nearly  the  height  of  the  previous  day.  Subsequent 
mixing  height  measurements  yield  a  vertical  particle  distribution  that 
is  an  indicator  of  the  previous  mixing  height  maximum  and  not  neces¬ 
sarily  the  mixing  height  at  the  time  of  measurement  (Ref  46:1183). 

Other  processes  are  often  involved,  however.  During  the  night 
a  surface  inversion  frequently  develops.  These  inversions  are  formed 
largely  by  turbulent  transfer  of  heat  down  toward  the  ground  (Ref  55:13). 
Thus,  in  many  cases  the  aerosol  layer  would  follow  the  theoretical 
pattern  of  a  decreasing  mixing  layer  depth.  Because  of  the  paucity 
of  aircraft  measurements  in  varying  meteorological  conditions  and 
because  of  the  complex  atmospheric  processes  involved,  this  author  feels 
that  a  real  time  determination  of  mixing  layer  height  is  impossible 
unless  there  is  detailed  knowledge  of  the  history  of  the  air  mass  and 
probably  upper  air  soundings. 
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Mixing  layer  depth  may  be  determined  for  such  a  study  as  this, 
however,  based  on  a  knowledge  of  the  meteorological  processes  for 
this  geographical  location,  Germany.  One  would  certainly  expect  that 
this  mixing  layer  would  not  apply  to  other  locations  such  as  Southeast 
Asia  where  the  thermal  and  mechanical  atmospheric  processes  differ  in 
magnitude.  This  model  also  fails  to  predict  the  hour  by  hour  and  day 
to  day  variations  of  the  mixing  layer  depths  but  it  is  representative 
of  the  real  vertical  aerosol  profile  based  on  mean  meteorological 
conditions . 

Aerosol  Mc/lei  and  Use 

Tie  model  consists  of  two  parts:  the  mixing  layer  aerosol 
attenuation  coefficient  and  the  aerosol  attenuation  coefficient  of  the 
clear  air  mass  above.  Key  to  the  use  of  the  model  is  ground-level 
(horizontal)  prevailing  visibility.  If  the  prevailing  visibility  is 
30  km  or  more  f  ;oteoro logical  range  hO  km  or  more)  then  the  vertical 
attenuation  coefficient  decreases  as  the  density.  The  aerosol 
attenuation  coefficient  at  any  altitude,  H  (in  km),  is  given  by 


where  oa(o)  is  the  ground  level  attenuation  coefficient.  There  is  a 
slight  difference  in  the  determination  of  the  ground  level  attenuation 
coefficient  for  use  in  this  vertical  model.  Since  there  is  an 


indication  the  upper  particle  size  limit  decreases  with  altitude  then 


the  continental  distribute 
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the  upper  particle  limit  of  lOp  (Ref  41:118).  The  relationship  between 
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the  aerosol  attenuation  coefficient  at  1 .0636p  and  .55y  is  modified 
only  si ightly 


0a(l.o636p)  =  .535  aa(.55y)  (48) 

where  the  particle  composition  (index  of  refraction)  is  assumed  to  be 
the  same.  Combining  Equations  (43)  and  (48),  the  ground  level  at¬ 
tenuation  coefficient  at  1 .0636p  is 

0a(0)  =  ijSZ-  (49) 

where  V  is  the  prevailing  visibility.  The  aerosol  attenuation  coef¬ 
ficient  at  any  altitude,  H,  then  becomes 


1.57 

V 


(*^7.99) 


(50) 


Using  Eq  (50)  then  the  attenuation  coefficients  are  determined  for  the 
different  altitudes.  This  model  follows  the  technique  used  by  others 
(Ref  49)  in  averaging  the  two  computed  values  to  determine  an  average 

I* 

attenuation  coefficient  for  a  1  km  layer.  These  average  values  are  then 
added  and  multiplied  by  sec  6  for  use  in  slant  path  problems.  The 
angle  0  is  the  angle  between  the  vertical  to  the  ground  and  the  angle 
to  the  transmitter/receiver.  Use  of  this  model  is  explained  in  Appendix 
C. 

If  the  prevailing  visibility  is  less  than  30  km  then  the  upper 
air  mass  is  the  clear  air  model  for  a  visibility  of  30  km,  but  the 
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ground  layer  attenuation  coefficient  is  determined  by  the  horizontal 
prevailing  visibility.  This  layer  is  assumed  to  be  uniformly  mixed. 
Thus,  the  attenuation  coefficient  for  the  layer  is  determined  by  Eq 
(45) •  The  depth  of  the  layer  varies,  but  useful  computations  may  be 
made  by  using  a  mean  mixing  depth  of  1500  m  (Ref  30).  These  are 
characteristics  daytime  values  which  are  probably  very  realistic  for 
midmorning  to  late  afternoon  cases.  For  early  morning  and  late  evening 
calculations  where  a  low  inversion  is  known  to  exist,  a  mixing  depth 
of  300  m  offers  a  good  approximation,  although  this  may  be  high  for 
some  nocturnal  inversions  (Ref  60) . 
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This  chapter  discusses  briefly  the  effects  of  other  atmospheric 
phenomena  on  laser  radiation  at  1.0636p.  This  discussion  is  limited 
to  the  attenuating  effects  of  clouds,  fog  and  rain.  Other  phenomena 
are  mentioned  but  are  not  covered  in  any  detail.  In  addition,  the 
sensitivity  analyses  accompanying  previous  discussions  in  this  report 
are  not  included  in  this  chapter.  The  purpose  of  this  chapter  is,  then, 
to  provide  a  useful  method  for  estimating  the  attenuation  at  1 ,063bp 
due  to  clouds,  rain,  and  fog.  Further,  it  is  desirable  that  the 
attenuation  be  related  to  some  meteorological  observable. 


Attenuation  by  Rain 

Attenuation,  as  we  have  discussed,  is  dependent  on  droplet 
size  and  droplet  size  distribution.  For  rain  both  size  and  size 
distribution  are  highly  variable.  Jn  fact,  different  rain  droplet 
distributions  occur  within  the  same  rain  shower.  Since  these  droplet 
size  distributions  are  also  virtually  unmeasurable,  exact  Mie  cal¬ 
culations  have  very  little  meaning. 

For  distributions  where  the  droplet  radius  is  large  compared 
to  the  wavelength  of  radiation,  the  geometrical  approximation  to  Hie 
theory  may  be  used.  For  rain  drops,  whose  radius  is  approximately  1 
mm  or  larger  (Ref  13:727),  the  extinction  cross  section  asymptotically 
approaches  two.  In  addition,  the  imaginary  index  of  refraction  of 
water  is  on  the  order  of  10“G  for  wavelengths  in  the  visible  and 
1.0636jj;  hence,  the  absorption  is  neglible. 
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Because  size  distributions  were  not  definable  and  the  geo¬ 
metrical  approximation  is  valid  for  rain,  Chu  and  Hogg  (Ref  13)  found 
instructive  the  calculation  of  attenuation  coefficients  for  rain 
droplet  distributions  having  a  single  radius.  For  this  ideal  case, 
the  attenuation  coefficient  becomes 

a  =3.2.5^extW  (51) 

R  r 

where  c  =  total  attenuation  in  dB/km, 

R 

r  =  drop  size  radius  in  microns, 

Cext  =  scattering  cross  section  efficiency, 
and  W  -  liqsid  v/ater  content  in  mg/m3. 

Using  the  asymptotic  value  of  two  for  Qext  and  the  conversion  factor 
for  dB/km,  then  Eq  (51)  becomes 

oR  =  w  (km-1)  (52) 

which  is  in  terms  agreeable  with  those  for  other  attenuation  coefficients 
in  this  report. 

The  energy  intercepted  by  precipitation  particles  is,  however, 
partially  scattered  forward  (toward  the  receiver)  and  this  reduces 
the  attenuation  coefficient  determined  by  Eq  (52).  When  scattering  is 
a  maximum  then  so  is  the  forward  scatter  component.  The  larger  the 
scattering,  then  the  more  that  Eq  (52)  is  in  error.  This  problem  is 
more  pronounced  in  cases  of  heavy  rain  (Ref  13:739).  A  forward 
scattering  correction  factor  has  been  derived  in  Reference  13  assuming 
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a  gaussfan  beam  profile  with  an  effective  aperture  equal  to  the  geo 
metrical  cross  section  of  the  water  sphere.  This  correction  factor 
is 


$ 


I  Qsca  rl  _ : _  d2 

L  .  Qext  ^  w  2  +  w2  /L\  2 
o  s  (t) 


(53) 


where  L  =  path  length, 

z  =  distance  to  the  scattering  particles  from  the  source, 
to  **  half-width  of  the  transmitted  beam  pattern  (in  cm), 

(i)  «=  radius  of  the  water  droplet  (in  cm), 

3 

and  =  half-width  of  the  beam  pattern  at  the  receiver  (in  cm) 
if  no  attenuation  is  assumed.  For  visible  and  near  infrared,  there 
is  no  absorption  so  that  Qsca/Qext  becomes  one.  Chu  and  Hogg  in 
Reference  13  have  made  calculations  of  the  correction  factor  assuming 
the  minimum  beam  radius  at  the  uniphase  aperture  to  be  .25  cm  (Ref 
13:738) .  An  additional  assumption  of  this  forward  scattered  correction 
factor  is  that  this  scattered  energy  rejoins  the  original  beam  (Ref  13: 
736).  Integrating  over  different  rain  drop  sizes  then  one  obtains  the 
attenuation  curves  depicted  in  Figure  21.  The  correction  factor  at 
1.0636y  will  be  the  same  as  that  calculated  at  .63p  by  Chu  and  Hogg 
since  there  is  no  absorption  in  either  case. 

The  exact  correction  factor  for  the  laser  designator  problem 
■Is  Incalculable  because  of  the  complex  geometry  in  each  case.  If  one 
assumes  a  similar  correction  factor  applies  in  this  problem,  then  a 
crude  multiple  scattering  attenuation  coefficient  may  be  calculated 
for  rain. 
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Fig.  21.  Rain  Attenuation  Coefficient  Curves  for  Dif¬ 
ferent  Rain  Rates,  Curves  Calculated  from  Measure¬ 
ments  Made  over  2.6  Km  Path  (Ref  13:734). 

Note:  [dB/Km  *  (4. 34)  =  <y(km_1)j 


Fig.  22.  Deirmendjian  Cloud  Size  Distribution, 
Cl  (Ref  17:80) 
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Meteorological  Parameters  for  Rain  Attenuation 

To  calculate  the  attenuation  in  rain,  a  knowledge  of  both 
rain  drop  size  and  liquid  water  content  (mg/m3)  is  needed.  As  we 
discussed  previously,  precise  droplet  size  is  unknown,  but  a  median 
droplet  size  may  be  related  to  the  source  mechanism.  The  source" 
mechanisms,  convective  showers  or  stratus  drizzle  or  rain,  can  exist 
separately  or  in  combination. 

It  is  an  accepted  atmospheric  phenomena  that  the  droplet  size 
and  size  distribution  is  modified  as  the  droplet  falls  from  cloud  to 
ground.  The  large  droplets  (2500p  or  more)  break  up  and  become 
smaller  droplets;  the  very  small  droplets  evaporate  shifting  the  drop¬ 
let  size  distribution  measured  at  the  ground  level  to  larger  median 
values.  For  the  purpose  of  making  an  approximation,  however,  the 
droplet  size  wili  be  assumed  to  be  constant  from  cloud  to  ground. 
Typical  median  droplet  sizes  as  a  function  of  source  mechanism  are 
given  in  Table  XIV. 


Tuble  XIV 


Typical  Median 
as  a  Function 

Rain  Droplet  Sizes 
of  Source  Mechanism 

Cloud  Condition 

Drop  Size  (y) 

Low  Stratus 
ft  2000  ft  AGL) 

500* 

Middle  Stratus 
(ffc  12000  ft  AGL) 

800 

Convection 

1500 

*These  droplet  sizes  are  valid  for  vicinity  of  Cloud  Base,  but  for 
this  model  these  sizes  are  assumed  uniform  to  ground  level. 
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These  values  from  Reference  10  are  the  result  of  a  comprehensive 
literature  search  by  Captain  Caratula.  The  median  stratus  values 
were  obtained  from  Reference  1  (Page  14)  and  Reference  11;  the  con¬ 
vective  values  were  obtained  from  Reference  44  (Page  357)*  It  must 
be  emphasized  that  these  values  are  necessarily  somewhat  arbitrary, 
but  nevertheless  representative  of  these  source  mechanisms. 

The  water  content  is  determined  by  the  reported  rain  rate  and 
cloud  type  (Source  Mechanism).  The  water  content,  W(mg/m3),  for 
stratiform  rainfall  is  (Ref  44:369) 

W  =  72(R)0*88 

where  R  =  Rainfall  rate  in  mm/hr.  A  similar  relationship  for 
convective  rain  is  (Ref  44:369) 

/  n°-97 

W  »  52  (R) 

where  R  =  Rainfall  rate  in  mm/hr.  Rainfall  reported  by  the  meteoro¬ 
logist  is  categorized  as  light,  moderate,  or  heavy.  These  categories 
are  given  in  Table  XV  where  the  metric  equivalent  values  are  noted. 

Since  the  rainfall  categories  encompass  a  range  of  rain  rates,  a  mean 
attenuation  coefficient  for  each  category  is  meaningless.  Calculations 
were  made  using  the  best  and  worst  cases  in  the  particular  categories; 
the  results  of  those  calculations  are  shown  in  Table  XVI.  These  values 
are  uncorrected  for  the  forward  scattered  component  since  the  correction 
factor  depends  on  the  exact  rain  rate.  For  the  light  rain  category 
the  correction  factor  varies  from  0  to  approximately  10%  (Figure  21); 


(54) 


(55) 
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for  the  moderate  rain  category  ft  varies  from  10  to  about  20%  (Figure 
21).  In  extremely  heavy  rain,  the  correction  might  reduce  the  at*“ 
tenuation  by  almost  50%. 

Table  XV 

Reported  Rainfall  Categories 

Classification  Rainfal 1  Rate  Metric  Eg 

Light  R  <  .1  in/hr  R  <  2.54  HL 

Moderate  .  1  III  <  R  <  .3  in/hr  2.54  HI  1  R  <  7-62  — 

hr  “  hr  nr 

Heavy  R  >  .3  in/hr  R  >  7*62 

Table  XVI 


Attenuation  Coefficients  for  Rain  (Km-1) 

Ey  Source  Mechanism  and  Reported  Rainfall  Category 


Classif ication 

Low  Stratus 

Middle  Stratus 

Convection 

Light 

0-9.226 

0-5.766 

0-2.415 

Moderate 

9.226-24.26 

5.766-15.162 

2.415-7.011 

Heavy 

>  24.26 

>  15.162 

>  7-011 

Values  calculated  from  values  in  Tables  XIV  and  XV  using  Eqs  (52),  (54), 
and  (55). 


Attenuation  by  Fog 

Fog  is  by  international  agreement  a  water  droplet  aerosol  that 
reduces  visibility  beiov.  i  km  (Ref  8:223).  The  fog-haze  boundary, 
however,  has  been  identified  by  Foitzik  (1938),  and  later  confirmed  by 
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Neiburger  and  Chien  (I960),  Bullrich  (1963)  and  Eldridge  (1969) >  as 
being  at  a  meteorological  range  of  about  1.2  km  (Ref  28:3).  For  this 
report,  a  fog  is  identified  as  a  water  aerosol  reported  by  the  meteor¬ 
ologist  regardless  of  the  visual  range. 

Fogs  in  general  are  quite  inhomogeneous  which  Is  due  in  large 
part  to  their  formation  mechanism.  In  addition,  the  vertical  extent 
of  fog  may  vary  from  a  few  feet  to  hundreds  of  meters.  Also  the 
vertical  extent  of  fog  is  rarely  reported  by  the  meteorologist  except 
as  a  general  classification  of  fog  (i.e.,  shallow  fog).  For  these 
reasons,  fog  as  reported  by  the  meteorologist  is  of  little  practical 
value  to  this  operational  problem  especially  if  the  sky  is  totally 
obscurred. 

Some  calculations  were  accomplished,  however,  for  the  purpose 
of  establishing  representative  attenuation  values  that  could  be 
compared  to  other  attenuation  coefficients  previously  calculated. 

This  model  of  the  attenuation  coefficient  assumes  that  fog  is  homo¬ 
geneous  in  both  the  horizontal  and  vertical  directions,  it  also 
assumes  that  the  visibility  recorded  in  fogs  is  the  true  visibility 
or  meteorological  range. 

In  addition,  the  attenuation  coefficient  in  fog  is  assumed  to 
be  independent  of  wavelength  for  the  visible  and  near  infrared.  This 
is  an  assumption  that  was  found  to  be  reasonably  accurate  by  Reference 
h.  This  author  theoretically  verified  this  fact  using  several  fog 
models  used  by  others  (Ref  4:*»7)  •  For  the  four  models  investigated, 
the  attenuation  at  1.0636p  v/as  larger  than  at  ,55u,  but  by  only  about 
3%.  The  attenuation  coefficient  for  fog,  Op,  may  be  determined  by 
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3.912  (56) 

F  ~  V 

where  V  Is  the  visibility,  it  has  been  determined  experimentally  that 
this  relationship  is  valid  only  for  visibilities  of  more  than  200  m 
(Ref  4 : 65) .  For  visibilities  less  than  this,  multiple  scattering  be¬ 
comes  extremely  important. 

Attenuation  by  Clouds 

Clouds  like  fog  are  the  greatest  attenuators  of  later  radiation. 
Like  rain  and  fog,  the  cloud  droplet  size  distribution  is  quite  dif¬ 
ficult  to  describe.  The  water  content  and  other  properties  of  clouds 
vary  with  cloud  type  and  formation  mechanism.  There  are,  however, 
several  cloud  models  available  and  one,  in  particular,  in  wide  general 
use.  This  model  cloud  is  the  Deirmendjian  Cloud  Cl,  which  represents 
a  cumulus  cloud  of  moderate  thickness. 

This  cloud  model  size  distribution  is  shown  in  Figure  22.  The 
liquid  water  content  of  this  cloud  model  is  0.063  gm~3  (Ref  17:79). 

It  would  have  been  desirable  to  investigate  other  cloud  distributions. 
Time  constraints  on  this  project,  however,  limited  this  investigation 
to  the  one  cloud  size  distribution  described  by 

n(r)  =  2.373  r^exp(-1.5r)]  (57) 

where  r  is  the  droplet  radius  in  microns  and  the  particle  density  is 
assumed  to  be  100  cm-3.  Mie  calculations  using  this  droplet  size 
distribution  and  refractive  indices  of  water  gave  approximately  the 
same  attenuation  coefficient  as  was  calculated  for  visibility  of  250 
m  in  fog  (see  Table  XV! I). 
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Table  XVII 


Cloud  and  Fog  Attenuation 
Coefficients  (Km-1) 


Wavelength (u)  Attenuation  Coefficient 


.45 

16.34 

.55 

16.57 

.70 

16.82 

1.0636 

17.34 

Fog 

Attenuation  Coefficient 

1.956 
2.608 
3.912 

4.347 
4.890 
5.583 
6.520 
7.824 
9.780 
13.040 
19.560 

Cloud  Attenuation  Coefficients  were  calculated  using  Hie  calculations 
and  the  Dei rmendj ian  Cl  distribution.  Fog  attenuation  coefficients 
were  calculated  using  Eq  (56). 

This  cloud  model  indicates  that  clouds  may  in  general  be 
strong  attenuators  of  radiation.  In  any  case,  transmission  paths  of 
several  hundred  meters  through  such  a  cloud  would  greatly  reduce  the 
amount  of  radiation.  For  example,  a  path  length  of  100m  through  this 
cloud  reducesthe  energy  to  about  18%  of  the  original.  Such  attenuation 
Is  an  important  factor  in  determining  the  maximum  acquisition  and 
tracking  range  of  the  laser-guided  weapon  system. 
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Other  Attenuating  Mechanisms 

Several  other  attenuators  are  worthy  of  mention  at  this  point. 
Some  such  as  snow  are  highly  relevant  to  this  discussion  because  of 
their  occurrence,  but  others  such  as  hail  were  not  mentioned  because 
of  their  rarity. 

Snow,  freezing  rain,  and  sleet  are  of  the  most  important  of 
these  other  attenuating  mechanisms.  They  are  important  because  of 
their  occurrence  in  this  geographical  region.  Unfortunately,  the 
assumption  we  have  made  for  our  solution  using  Mie  Theory — that  is, 
that  the  particle  is  spherical — is  not  valid  for  snow.  In  addition, 
there  are  the  other  complicating  factors  (such  as  snow  rate)  like 
those  encountered  for  rain.  Other  phenomena  such  as  blowing  dust, 
hail,  and  others  were  not  discussed  because  of  their  rarity  in  Germany. 

Backscatter 

This  chapter  would  be  quite  incomplete  without  the  mention  of 
backscatter.  This  term  is  generally  used  to  describe  the  problem  of 
that  component  of  the  scattered  radiation  that  is  scattered  directly 
backward.  Backscatter  becomes  a  serious  problem  when  the  laser 
designator  (transmitter)  and  receiver  are  in  close  proximity.  The 
problem  is  perhaps  best  illustrated  by  the  following  example. 

Backscatter  is  a  common  problem  while  motoring  in  heavy  fogs 
at  night.  If  the  driver  actuates  his  high  beams  then  he  is  able  to 
see  considerably  less  than  when  he  was  using  his  low  beams.  The  reason 
for  this  is  backscattering;  that  is,  more  light  is  being  backseat te red 
into  the  driver's  eyes  reducing  his  overall  sensitivity.  Backscattering 
is  improved  by  separating  the  viewing  angle  and  the  illumination  plane 
(lights).  Separation  of  the  receiver  and  transmitter  is  the  precise 
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reason  for  fog  lights  on  motor  vehicles. 

This  same  problem  exists  in  the  event  the  receiver  and  trans¬ 
mitter  for  the  laser  target  designator  system  are  in  close  proximity. 
This  problem  increases  in  magnitude  as  the  amount  of  scattering 
increases;  the  larger  the  scattered  radiation  the  larger  the  back- 
scattering.  Backscattered  radiation  is  detected  by  the  receiver  and 
therefore  becomes  an  undesirable  background  signal.  Unfortunately 
backscattering  is  largest  when  detector  sensitivity  must  be  an  optimum. 
When  the  scattering  is  large  as  in  fogs  and  clouds,  the  magnitude  of 
the  received  radiation  (reflected  from  the  target)  is  quite  small. 

The  receiver  then  is  not  likely  to  sense  the  target  signal  at  all. 

The  magnitude  of  the  backscatter  component  is  dependent  on 
the  angle  subtended  at  the  receiver  by  the  scatterer  (droplet).  Since 
the  geometry  of  each  backscatter  problem  is  extremely  complex,  then 
the  backscatter  is  generally  given  in  terms  of  a  coefficient  expressed 
in  km_1/steradian.  Representative  values  of  the  backscatter 
coefficient  are  given  in  Table  XVII I . 

Use  of  Attenuation  Coefficients  ? ri  Model 

In  this  chapter,  we  have  derived  attenuation  coefficients  that 
estimate  the  attenuation  of  laser  radiation  in  clouds,  rain  and  fog. 

All  of  these  coefficients  are  not  additive  to  yield  the  extinction 
coefficient.  Clouds  and  fog  are  formed  on  condensation  nuclei,  for 
example,  and  therefore  reduce  the  number  of  aerosols.  In  fact,  as  a 
good  approximation,  the  attenuation  of  aerosols  is  neglected  in 
calculating  the  atmospheric  attenuation  in  fogs  and  clouds.  The 
atmospheric  extinction  coefficient  in  clouds  and  fog,  then,  is  the  sum 


131 


GEP/Pli*/7/*“Ii 


Table  XVII! 

Representative  Per  Particle  Backscattering  Coefficients 


for  Various  Models 

at  -55y  and  1  .063u 

Refractive 

Wavelength 

Backscattering 

1  ndex 

in 

Coefficient  in 
(Knrtyster .) 

Model 

R 

1. 

Microns 

Continental 

1.55 

.048 

.55 

.0153 

Haze 

1.53 

.072 

1,0636 

.0146 

{.03p-»20p) 

1.53 

.031 

•55 

.0194 

1.51 

.046 

1.0636 

.0167 

Marctime 

1.53 

.031 

.55 

.0300 

Haze 

1.51 

.046 

1.0636 

.0141 

Haze 

1.53 

.031 

.55 

.0262 

L 

1.51 

.046 

1,0636 

.0116 

Cloud 

1.334 

.0 

.55 

.0505 

Cl 

1.325 

.00001 

1.0636 

.0625 
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of  the  molecular  attenuation  coefficient  and  the  fog  or  cloud  at¬ 
tenuation  coefficient  discussed  in  this  chapter.  Since  the  former 
is  generally  quite  small  compared  to  the  later,  a  reasonable  approxi¬ 
mation  of  the  atmospheric  attenuation  coefficient  in  clouds  and  fogs 
Is  given  by  the  value  from  Table  XVII  and  Eq  (56)  respectively. 

Unlike  the  cloud  and  fog  attenuation  coefficients,  the 
attenuation  coefficient  for  rain  is  additive.  Although  there  is 
evidence  that  rain  does  remove  atmospheric  particulate  matter  by  the 
processes  of  rainout  and  washout  discussed  earlier,  the  efficiency  of 
these  removal  processes  is  relatively  low  (Ref  32:161).  One  of  the 
findings  of  a  six  year  study  in  the  United  States  of  atmospheric 
attenuation  coefficients  is:  precipitation  in  itself  does  not  appear 
to  have  a  noticeable  effect  in  cleaning  the  air.  Such  noticeable 
atmospheric  clearings  are  associated  with  frontal  passage  occurring 
simultaneously  with  the  rain  (Ref  39-961, 962) .  Calculation  of 
atmospheric  attenuation  in  rain,  then  is  the  result  of  the  attenuation 
by  the  molecules,  aerosols,  and  raindrops. 
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3 X .  Summary,  Conclusions,  and  Recommendations 
Summary  and  Model  Lise 

In  summary,  this  report  develops  a  model  for  prediction  of 
maximum  laser  target  designator  system  range  or,  stand-off  distance 
where  atmospheric  attenuation  of  1 .0636)1  laser  radiation  Is  the 
limiting  system  factor.  This  report  develops  an  optimum  lock-on 
range  equation  in  which  atmospheric  conditions  and  system  geometry 
are  the  only  variables.  Major  assumptions  made  in  the  development 
of  this  equation  are:  (l)  the  laser  designator  and  receiver  have 
specific  electrical  and  optical  design  characteristics,  (Z)  the  laser 
beam  is  reflected  only  from  the  target,  (3)  the  target  is  a  diffuse 
Lambertian  surface  with  ten  percent  ref lectivi ty, (4)  the  receiver  is 
designed  to  electrically  discriminate  against  continuous  background 
signals,  and  (5)  the  target  is  within  the  field  of  view  of  the 
receiver.  Jn  addition,  it  is  assumed  that  the  Beer-Lambert  Law  des¬ 
cribes  the  atmospheric  transmittance  of  laser  radiation. 

Maximum  lock-on  range  or  standoff  distance  is  reduced  to  a 
determination  of  atmospheric  attenuation  due  to  each  of  the  atmospheric 
attenuating  mechanisms;  aerosol  scattering  and  absorption,  molecular 
absorption,  molecular  (Rayleigh)  scattering,  and  scattering  of  clouds, 
rain,  and  fog.  Theoretical  calculations  of  the  coefficients  ‘which 
describe  these  attenuating  mechanisms  are  dependent  upon  four 
important  assumptions:  (1)  the  laser  radiation  is  monochromatic, 

(2)  the  atmosphere  is  characteristic  of  Central  European  atmospheric 
conditions,  (3)  the  atmosphere  is  horizontally  stratified  above  a  flat 
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earth,  and  (4)  the  one  kilometer  thick  vertical  layers  are  assumed  to 
be  homogeneous.  In  addition,  calculation  of  attenuation  due  to 
scattering  is  limited  to  single  scattering  theory. 

Theoretical  calculations  show  that  the  greatest  attenuators 
of  1.063611  laser  radiation  are  clouds,  rain  and  fog.  For  dry  atmos¬ 
pheric  conditions,  the  aerosol  is  the  largest  attenuator,  followed  by 
(usually  by  an  order  of  magnitude  or  more)  molecular  absorption  and 
by  molecular  (Rayleigh)  scattering  which  is  negligible.  Each  of  these 
coefficients  is  related  to  the  physically  measureable  or  observable 
ground  level  meteorological  quantities  listed  below: 

a  (Molecular  attenuation  coefficient)  Relative  Humidity 

Temperature 

Pressure 


ag  (Aerosol  attenuation  coefficient)  Prevailing 

visibility 

Op  (Fog  attenuation  coefficient)  Prevailing 

visibility 

oR  (Rain  attenuation  coefficient)  Rain  Rate 

cloud  type 

oc  (Cloud  attenuation  coefficient)  None 

(For  this  model - 
the  attenuation 
coefficient  is 
constant  for  al 1 
clouds.) 


Vertical  distributions  of  the  molecular  and  aerosol  coefficients  are 
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Use  of  this  atmospheric  model  requires  knowledge  of  the  fol 
lowing  parameters  which  are  listed  together  with  the  attenuation 
coefficient  they  determine: 


Attenuation  Coefficient 
Molecular 

Aerosol 

Fog 

Rain 

Cloud 


Parameter 
Season  of  year 

Prevai 1 ing  visibi 1 ity 
(Real  Time  Mixing  Height 
Desired) 

Prevailing  visibility 
Rain  Rate 

(cloud  type  desired) 

None 


The  molecular  attenuation  coefficients  are  determined  using  Table 
XXIV,  Page  152.  The  boundary  (mixed)  layer  aerosol  attenuation 
coefficients  are  determined  using  Equation  (45)  and  the  aerosol  at¬ 
tenuation  coefficients  for  the  air  mass  above  are  determined  using 
Table  XXVIII,  Page  164.  The  cloud  attenuation  coefficient  is  obtained 
from  Table  XVII,  Page  129.  Fog  attenuation  coefficients  are  obtained 
from  either  Table  XVII,  Page  129, or  Equation  56,  Page  128.  Rain 
attenuation  coefficients  are  obtained  from  Table  XVI,  Page  126. 

Each  individual  attenuation  coefficient  is  summed  over  the 
vertical  path.  The  total  individual  sums  (which  are  now  individual 

vertical  optical  paths  and  no  longer  coefficients)  are  new  added  to 

(_ 

yield  a  totai  attenuation  vertical  optical  path,  oily,  as  chown  below: 
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oHt 

=  a H  molecular  +  oH  aerosol 

For  Clouds : 

<jHt 

=  aH  molecular  +  aH  cloud 

For  Fog: 

CjHy 

=  aH  molecular  +  aH  fog 

For  Rain: 

oHt 

=  crH  molecular  +  aH  aerosol  +  aH  Rain 

where  aH  is  the  attenuation  coefficient  summed  vertically  to  the 
required  vertical  height,  H. 

To  calculate  the  slant  path  attenuation,  then  multiply  the 
resulting  aHj  by  secant  0,  where  0  is  the  angle  between  the  vertical 
and  the  line  connecting  the  target  and  designator  or  receiver.  For 
cases  in  which  the  receiver  and  designator  are  not  co-located  then 
each  optical  path  must  be  calculated  separately.  The  total  optical 
depth  for  the  receiver  cfcr  (where  crfcr  =  crHy  x  secant  0)  and  total 
optical  depth  for  the  designator  a £(j  (where  atd  =  oHy  x  secant  0)  are 
then  substituted  into  Equation  (67),  the  laser  lock-on  range  equation. 

Use  of  the  model  is  demonstrated  in  Appendix  C. 

Concl us  ions 

Several  conclusions  may  be  drawn  from  the  results  of  this 
Investigation: 

1.  Current  models  for  the  radiative  transfer  of  laser  radiation 
through  the  atmosphere  are  limited  in  their  predictions  to  single 
scattering  solutions  with  crude  multiple  scattering  correction  factors 
in  some  cases.  The  model  developed  in  this  report  is  limited  to  single 
scattering  theory  and  hence  may  underestimate  the  maximum  laser  lock-on 
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2.  The  study  of  attenuation  of  laser  radiation  in  the 
atmosphere  due  to  atmospheric  molecular  gases  Is  quite  far  advanced 
compared  to  the  study  of  aerosol  attenuation.  Molecular  scattering 
may  be  ignored  for  1.0636p  radiation.  Molecular  absorption  in  this 
small  spectral  region  (1 .O63-I .065p)  is  not  extensively  complicated 
by  overlapping  molecular  absorption  lines,  but  is  a  function  of  water 
vapor  absorption  lines.  Hence,  molecular  absorption  is  dependent 
upon  atmospheric  water  vapor  concentration. 

3.  Aerosol  attenuation  of  1.0636y  radiation  is  by  comparison 
to  molecular  attenuation  a  more  important  factor  in  limiting  maximum 
laser  lock-on  range.  Precise  knowledge,  however,  concerning  aerosol 
particle  size  distributions,  number  densities,  and  vertical  distri¬ 
butions  is  lacking  because  of  the  difficulty  in  measuring  such 
quantities.  Of  particular  consequence  for  this  tactical  problem  is 
the  vertical  aerosol  distribution.  The  model  developed  in  this  report 
employs  the  boundary  layer  concept  which  is  suggested  by  recent  ex¬ 
tensive  aircraft  measurements.  The  boundary  layer  is  an  aerosol 
particle  dense  layer  where  the  particles  are  approximately  uniformly 
mixed.  Above  the  boundary  (mixed)  layer  the  air  mass  is  relatively 
clear  and  the  particle  concentration  markedly  less.  The  height  of  such 
a  layer  is  variable  and  difficult  to  measure  from  ground  based 
observations.  In  fact,  hourly  predictions  of  aerosol  attenuation 
coefficients  are  impossible  without  a  detailed  history  of  air  mass  and 
current  surface  and  upper  level  meteorological  conditions. 
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k.  The  aerosol  attenuation  of  1.0636p  radiation  is  quite 
sensitive  to  the  particle  size  distribution  used, for  particle  size 
distribution  establishes  the  primary  relationship  among  aerosol  at¬ 
tenuation  coefficients  at  various  wavelengths.  Particle  composition 
(index  of  refraction)  and  upper  and  lower  particle  size  limits  are 
less  important  factors. 

5.  The  total  aerosol  attenuation  coefficient  at  1.0636p  is 
related  to  atmospheric  particle  concentration  which  is  determined  by 
measuring  prevailing  visibility.  Prevailing  visibility,  however,  is 
a  rather  inexact  meteorological  observable;  it  depends  on  many 
scientifically  unmeasureable  factors  including  the  local  terrain, 
visual  objects  available,  and  visual  acuity  of  the  observer.  Such 
visual  observations  would  not  be  highly  reliable  for  real-time 
operational  calculation  of  attenuation.  In  fact,  such  observations 
most  probably  would  not  be  available  for  hostile  areas. 

6.  Calculation  of  attenuation  due  to  tin  scattering  by  clouds, 
rain,  and  fog  are  based  ori  assumptions  of  homogeneity  of  the  scattering 
media  and  other  assumptions  concerning  particle  size  distribution. 

Such  calculations  then  are  necessarily  crude,  but  can  be  relied  upon 

to  yield  relative  attenuation  magnitudes  and  estimates  of  maximum 
lcck"on  range. 

Recommendations 

The  following  are  problem  areas  and  should  be  investigated  to 
improve  knowledge  In  th>s  area  and  increase  the  reliability  of 
estimating  laser  lock-on  ranges: 
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1.  The  model  developed  In  this  report  is  limited  because 
of  the  assumptions  made  in  its  development.  These  assumptions  may 
only  be  checked  and  the  model  validated  through  extensive  meteoro¬ 
logical  measurements  made  during  the  course  of  flight  tests.  V/eapon 
systems  may  be  properly  designed  and  evaluated  only  with  a  complete 
knowledge  of  the  environment  in  which  they  operate.  The  lack  of 
proper  and  reliable  meteorological  data  for  flight  weapons  tests  to 
date  is  quite  evident. 

2.  Since  clouds  are  a  common  atmospheric  phenomena  and  they 
are  also  a  chief  attenuation  mechanism  then  atmospheric  data  records 
of  cloud  cover  and  height  should  be  correlated  with  current  cloud-free 
1 ine-of-s ight  methodology  to  contribute  to  determination  of  lock-on 
range. 

3.  The  effectiveness  of  any  v/eapon  system  in  any  geographical 
region  may  be  evaluated  as  a  function  of  atmospheric  conditions  only 
in  light  of  the  relative  occurrence  of  the  different  conditions.  Such 
information  with  cloud  data  included  was  requested  by  this  author, 
but  due  to  technical  delays  was  not  available  for  inclusion  in  this 
report. 

1}.  This  model  like  other  laser  attenuation  models  is  limited 
to  the  assumption  of  single  scattering.  This  assumption  may  limit  the 
effectiveness  of  these  models  to  predict  laser  lock-on  ranges  when  a 
diffuse  Lambertian  target  is  assumed.  If  typical  military  targets  are 
diffuse  reflectors  (as  current  literature  suggests),  then  the 
applicability  of  the  Beer-Lambert  law  to  such  reflected  radiation  must 
be  investigated. 
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APPENDIX  A 

Molecular  Absorption  Calculations  and  Tables 

For  absorption  coefficient  calculations  in  the  1.0636y  case, 
Eq  (23)  reduces  to 


where  t.ie  i  subscript  indicates  the  i  lines  of  oxygen  and  k  the  k 
lines  of  water  vapor.  The  w  and  o  indicate  water  vapor  and  oxygen 
respectively. 

The  rotational  partition  function  is  proportional  to  (T/236) 
in  the  case  of  oxygen  and  (T/296)1’5  in  the  case  of  water  vapor. 

This  has  been  included  in  Eq  (58).  The  various  values  for  the 
vibrational  partition  function  may  be  determined  from  Table  1XX.  The 
intensities,  ground  energy  levels,  and  line  half  widths  at  half 
heights  are  given  in  Table  XX. 
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Table  IXX 

Vibrational  Partition  Functions  (Ref  47:4) 


Molecule 

i 

Temperature  175 

200 

225 

250 

275 

296 

325 

[ 

• 

h2o 

1.5 

1.000 

1.000 

1 .000 

1.000 

1 .000 

1.000 

1  .001 

f 

°o 

1.0 

1.000 

1 .000 

1 .000 

1 .000 

1.000 

1.000 

1  .001 

Table  XX 


Absorption  Line  Parameters 
(Ref  50) 
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Wavenumber (cm-1) 

S(cm_1/molecm“2) 

^(cm-1) 

E'^cnr1) 

9385.570 

.1092-24 

.051 

1255.160 

9385.920 

.233E-24 

.076 

586.480 

9386.220 

•75SE-26 

.048 

81.581 

9386.960 

.106E-24 

.081 

224.830 

9388.236 

.596E-26 

.048 

79.565 

9388.980 

.542E-24 

.080 

885.620 

9390.556 

•702E-26 

.048 

130.437 

9392.596 

.580E-26 

.048 

128.398 

9392.604 

.432E-26 

.048 

16.388 

9393.460 

.660E-26 

.053 

1006.120 

9394.592 

.602E-26 

.048 

190.775 

9396.653 

.514E-26 

.048 

188.713 

9398.220 

.227E-25 

.063 

882.930 

9398.325 

.480E-26 

.048 

262.583 

9400.407 

.420E-26 

,048 

260.501 

9401 .756 

.357E-26 

.048 

345.850 

9401.760 

.345E-25 

.039 

1327.140 

9403.175 

.453E-26 

.048 

42.224 

9403.460 

.104E-25 

.038 

1327.140 

9403.857 

.319E-26 

.048 

343.749 

9404.882 

.249E-26 

.048 

440.565 

9407.002 

.225E-26 

.048 

438.444 

9407.701 

.163E-26 

.048 

546.712 

9409.070 

•54iE-24 

,061 

586.260 

9409.840 

. 14SE-26 

.048 

544.572 

Ckh  t  A  0  1  n 

^rs  •  v  *  iv 

OftOr 

.048 

DD4.Z/6 

9411.540 

.954E-26 

.065 

744.090 

9412.340 

•331E-26 

,041 

1033.510 

9412.368 

•921E-27 

.048 

662.118 

9412.407 

.576E-27 

.048 

793.239 

9412.820 

.417E-24 

.069 

704.220 

.^.  B . gjlbgj _ 
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Absorption  Lines 

All  the  absorption  lines  with  intensities  greater  than 
1.0  x  10“27(cm-1/molecules  cm-2)  are  shown  in  Figure  23.  The  number 
one  on  the  line  indicates  water  vapor;  a  seven  represents  oxygen. 

Absorption  Coefficient  Tables 

Tables  XXI  through  XXVII  show  the  absorption  coefficients 
(in  km”1)  for  the  different  altitudes  and  model  atmospheres  for 
seven  wavenumbers.  Tables  XXi  through  XXV  are  for  wavenumbers  in 
the  spectral  region  of  interest.  Table  XXVI  and  Table  XXV; 1  are 
for  wavenumbers  outside  the  spectral  region. 

The  following  relationship  between  wavenumber,  k,  and  wave¬ 
length,  X,  is  given  for  convenience  of  the  reader: 

X  =  loVk 

where  X  is  expressed  in  microns  (l  x  10"6m) 
and  k  is  expressed  In  cm”1 
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Table  XXI 

Molecular  Absorption  Coefficients 
for  K  =  9389.671*»  Wavenumbers 


Altitude 

(Km) 

Atmosphere  Models 

Midlat.  Midlat. 

Summer  Winter 

U.S. 

Standard 

Spring 

Fall 

0 

1 .480E-03 

4.210E-04 

6.700E-04 

0  -  1 

1.179E-03 

3.518E-04 

5.494E-04 

1  -  2 

6.908E-04 

2.378E-04 

3.503E-04 

2  -  3 

3.834E-04 

1 .61  OE-O^i 

2.183E-04 

3  *  *1 

2.062E-04 

1.057E-04 

1 .341E-04 

*i  -  5 

1.166E-04 

6.966E-05 

8.529E-05 

5  •  6 

7.005E-05 

4.879E-Q5 

5.669E-05 

6  •  7 

4.668E-05 

3.508E-05 

3.956E-05 

7  -  8 

3.233E-03 

2.585E-05 

2.857E-05 

Table  XXII 

Molecular  Absorption  Coefficients 
for  K  =  9398.220  Wavenumbers 

Atmosphere  Models 

U.S. 

Standard 

A1 ti tude 

Midlat. 

Midlat. 

Spring 

(Km) 

Summer 

Wi nter 

Fall 

0 

8.188E-03 

4.177E-03 

5.054E-03 

0  -  1 

7.187E-03 

3.777E-03 

4.542E-03 

1  -  2 

5.364E-03 

3.047E-03 

3.599E-03 

2  -  3 

3 .841 E-03 

2.426E-03 

2.784E-03 

3  “  4 

2.695E-03 

1 . 875E-03 

2.11  IE-03 

4-5 

1.920E-03 

1.426E-03 

1 .595E-03 

5  -  6 

I.393E-03 

1.090E-03 

1 .204E-03 

6-7 

1  nl.ic.ni 

0  ooorpAl. 

A  AA<I  r  _  A  l« 
j  .u;ut" u* 1 

7  -  8 

7.833E-04 

6.177E-04 

6.850E-04 
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Table  XXIII 

Molecular  Absorption  Coefficients 
for  K  =  9398.4962  Wavenumbers 


Altitude 

(Km) 

Midi  at. 
Summer 

Atmosphere  Models 

Midlat. 

Winter 

U.S. 

Standard 

Spring 

Fall 

0 

1.532E-03 

1  .432E-03 

1.393E-03 

0  -  1 

l  .358E-03 

1 .284E-03 

1  .255E-03 

1  -  2 

1.054E-03 

1 .017E-03 

1.005E-03 

2  -  3 

8.212E-04 

8.023E-04 

8.012E-04 

3  -  4 

6.432E-04 

6.312E-04 

6.353E-04 

4  ~  5 

5.078E-04 

4.961E-04 

5.021E-04 

5  -  6 

4.0UE-04 

3.878E-04 

3.949E-04 

6  -  7 

3 . 1 58E-0/I 

3.01  IE-04 

3.087E-04 

7  -  8 

2.475E-04 

2.320E-04 

2.398E-04 

Table  XXIV 

Molecular  Absorption  Coefficients 

for  K  = 

9401.780  V/avenumbers 

Altitude 

(Km) 

Midlat. 

Summer 

Atmosphere  Models 

Midlat. 

Wi  nter 

U.S. 

Standard 

Spring 

Fall 

2.498E-02 

i.406E-02 

1 .700E-02 

0  -  1 

2.298E-02 

1.376E-02 

1 .624E-02 

1  -  2 

1.939E-02 

1 .323E-02 

1.489E-02 

2  -  3 

1  .649E-02 

1.277E-02 

1 .378E-02 

3-4 

1  .444E-02 

1 .228E-02 

1.290E-02 

4-5 

1.319E-02 

1.183E-02 

1.227E-02 

5  -  6 

1 .24  IE -02 

1.148E-02 

1.180E-02 

6-7 

1.192E-02 

1.115E-02 

J.141E-02 

7  -  8 

1.150E-02 

1 .083E-02 

1 .105E-02 
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Table  XXV 

Molecular  Absorption  Coefficients 
for  K  =  9403.857  Wavenumbers 


A1 ti tude 
(Km) 

Atmosphere 

Midlat 

Summer 

Mode  1 s 

Midlat. 

Wi  nter 

U.S. 

Standard 

Spring 

Fall 

0 

1.116E-02 

1.092E“Q2 

1 .108E-02 

0  -  1 

! ,  1 12E-02 

1.088E-D2 

1 . 103E-02 

1  -  2 

< .  104E-02 

1 .080E-02 

1 .094E-02 

2  -  3 

1  .096E-02 

1 .073E-02 

1.084E-02 

3  -  4 

1 .038E-02 

1.063E-02 

1.074E-02 

4  -  5 

1  - 079E-02 

1.051E-02 

1 .062E-02 

5  -  6 

1.069E-02 

1 .038E-02 

1.050E-02 

6  -  7 

1 .059E“02 

1.024E-02 

1.036E-02 

7  -  8 

1  .046E-02 

1.008E-02 

1.020E-02 

Table  XXVI 

Molecular  Absorption  Coefficients 
for  K  =  9388.98  Wavenumbers 
(Outside  Spectral  Region) 


A1 ti tude 
(Km) 

Atmosphere 

Midlat. 

Summer 

Models 

Midlat. 

Wi  nter 

U.S. 

Standard 

Spring 

Fall 

0 

9.892E-02 

1.891E-02 

3.901E-02 

0  -  1 

8.470E-02 

1.676E-02 

3-397E-02 

1  -  2 

5.893E-02 

K299E-02 

2.485E-02 

2  -  3 

3-759E-02 

9.792E-03 

1 .710E-02 

3  -  4 

2.226E-02 

6.464E-03 

1 .096E-02 

4  -  5 

1  .297E-02 

3.779E-03 

6.819E-03 

5  -  6 

7.53&E-03 

2.238E-03 

4.140E-03 

6  -  7 

4.790E-03 

1.170E-03 

2.475E-03 

7  -  8 

2.940E-03 

5.054E-04 

l .430E-03 

1 
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Table  XXVI I 

Molecular  Absorption  Coefficients 
for  K  =  9^07.07  Wavenumbers 
(Outside  Spectral  Region) 


Altitude 

(Km) 

Midlat. 

Summer 

Midlat. 

Winter 

U.S. 

Standard 

Spring 

Fall 

0 

1.306E-01 

2.792E-02 

5.2S8E-02 

0  -  1 

1 . 138E-01 

2.A06E-02 

4.082E-02 

1  -  2 

8.021E-02 

1.975E-02 

3-5^5E-02 

2  -  3 

5.Z58E-02 

1.520E-02 

2.526E-02 

3  -  h 

3.218E-02 

1.029E-02 

1.6815-02 

k  -  5 

I.937E-02 

6.237E-03 

1.087E-02 

5  -  6 

1 .167E-02 

3.830E-03 

6.869E-03 

6  -  7 

7.692E-03 

2.06^E-03 

A.277E-03 

7  -  8 

4.910E-03 

9. 1 13E-04 

2.579E-03 

15  U 
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APPENDIX  B 

Summary  of  Mie  Theory  and  Calculation.** 

Outl ine  of  Hie  Theory 

In  1908,  Gustov  Mie  solved  the  problem  of  a  monochromatic 
plane  wave  incident  on  a  homogeneous,  isotropic  sphere  of  radius  r 
surrounded  by  a  transparent  homogeneous  and  isotropic  medium.  The 
incident  wave  causes  oscillations  of  the  free  and  bound  charges  in 
phase  with  the  applied  radiation.  These  oscillating  charges  (dipole) 
in  turn  create  secondary  electric  and  magnetic  fields.  These  fields, 
both  inside  and  outside  the  sphere,  are  expressed  by  Maxwell's  Equations 
and  a  formal  solution  of  the  problem  involves  solving  these  equations 
v/ith  the  appropriate  boundary  conditions. 

An  abbreviated  derivation  of  the  fields  inside  and  outside  the 
sphere  is  presented  in  Van  de  Hulst,  Chapter  9.  Basically,  the  field 
outside  the  sphere  is  composed  of  the  scattered  v/ave  and  the  incident 
v/avc.  Applying  the  boundary  conditions  and  the  conditions  to  be  satis¬ 
fied  at  infinity,  then  the  outside  scattered  and  incident  waves  are 
solved  in  terms  of  Legendre  Polynomials  and  Bessel  Functions.  Likewise, 
the  v/aves  inside  the  particle  are  represented  by  Legendre  Polynomials 
and  Bessel  Functions,  due  to  the  geometry  of  the  problem.  Matching  the 
boundary  conditions  of  the  tv/o  functions  now  completes  the  solution. 

The  field  at  any  point  inside  or  outside  the  sphere  may  be  expressed  in 
terms  of  known  functions. 

The  intensity  of  the  scattered  radiation  at  any  point  that  is 
a  great  distance,  R,  from  the  center  of  the  sphere  where  9  is  the 
scattering  angle  is  given  by 
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1o 

R2  87 r2 


^(m.x.O)  +  i  (m,x,e) 


(59) 


where  m  «  complex  refractive  index  of  the  sphere  with  respect  to 
the  medium 

x  =  27ir/.\,  the  dimensionless  size  parameter 
r  =  radius  of  the  sphere 

8  =  angle  between  the  incident  and  scattered  radiation 

!  =  Intensity  of  the  incident  radiation  (Ref  68:83). 

o 


The  intensity  functions  of  the  scattered  light,  i  and  i  ,  are 
respectively  the  components  perpendicular  and  parallel  to  the  plane 
formed  by  the  incident  and  scattered  beams.  The  determination  of 
the  intensity  by  Hie  Theory  is  generally  then  the  computation  of  the 
intensity  functions: 


i  -  |s(o)J2 

1  1 

1  -  |s  (e)  j2 

2  2  1 


(60) 


where  S  (o)  and  S  (6)  are  the  amplitude  functions. 

1  2 

The  field  components  are  di  reedy  proportion-!  to  these 
amplitude  functions.  The  amplitudes  are  given  by 


Si (8)  -  Z 
n=l 


2n+l 
n(n+l ) 


a  tt  (cos  e)+  b  T„(co.  8) 
n  n  n  n 


S2{0) 


(61) 

CO 

l  2n+l  _ 
n=l  n(n+l) 


b  tt.  (cos  8)+  a  m  (cos  0) 
n  n  n  n 
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where  a„  and  b  are  Riccati-Bessel  functions,  which  may  be  written 
in  terms  of  spherical  Bessel  functions,  and  where  7rn  and  Tn  are 
scattering  angle  functions.  These  scattering  angle  functions  are 
the  result  of  the  derivation  of  the  tangential  field  components. 
These  functions  are  written 


rn  (cose)  =  -  sin  6  Pn*  (cose) 


Tn  (cose)  = 


~  Pn*  (cose) 
d6  n 


where  Pn  are  associated  Legendre  Polymonials. 

The  efficiency  factor  is  defined  by  Van  de  Hulst  as  the  ratio 
of  the  outgoing  flux  to  the  incident  flux  per  unit  geometrical  cross- 
sectional  area  of  the  sphere.  These  factors  are  given  in  Chapter  9 
of  Reference  62  as: 


5  (2n+l )  C ja  I2  +  |b 

x2  n=l  n 


Qext  =  — ,  I  (2n+l)  Re  (an  +  bn)  (63) 

^  n=l 

Qabs  =  Qext  "QSCa 

where  Qsca  =  Scattering  efficiency  factor 

Qext  =  Extinction  efficiency  factor,  and 
Qabs  =  Absorption  efficiency  factor 
These  are  the  efficiency  factors  for  one  particle  of  unit  size, 
x  **  2ur/A.  The  foregoino  expressions  now  define  the  Mie  parameters 


from  which  all  other  parameters  needed  to  describe  intensity  end 
polarization  produced  by  a  single  sphere  can  be  derived. 

The  scattering  due  to  N  particles  of  radius  r  in  a  volume 

is  then 


&  -  tivr2  Q.sca 


where  the  area  is  mr2. 


Generally,  however,  there  are  many-sized  particles  in  a  volume  and 


the  scattering  coefficient,  f>,  for  a  volume  of  polydisperse  spheres 


|  ■“*  Qsca(r)n(r)  r2dr 

J 

r. 


where  n(r)  =  the  size  distribution  of  spheres, 


Q(r)  «  scattering  efficiency  at  radius  r. 


This  scattering  coefficient  is  generally  written  in  terms 


of  the  size  parameter,  x  =  2ur/A.  Then  Eq  (6j>)  becomes 


-gpr~  j  xn(x)  Qsca(x)dx 


The  size  distribution  function  n(x)  can  be  of  any  form.  For  aerosol, 


rain  and  fog  distributions,  the  function  n(x)  is  a  smoothly  varying 


function  such  as  the  log-normal,  modified  gamma  or  power  law  distri¬ 


bution  or  it  may  be  o  discontinuous  distribution.  These  distributions 


are  discussed  in  Chapter  5»  The  scattering  efficiency  function, 


Qsca,  on  the  other  hand.  Is  an  oscillatory  function  which  is 
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dependent  on  the  complex  index  of  refraction  and  the  size  parameter, 
x.  Figure  2h  shows  this  function  for  various  indices  of  refraction 
as  a  function  of  the  size  parameter. 

Mie  Theory  is  an  exact  computation  of  the  scattering  of  a 
plane  wave  incident  on  a  sphere.  There  are  two  approximations  to 
exact  Mie  Theory:  Rayleigh  scattering  and  Geometric  scattering.  The 
Rayleigh  is  an  approximation  to  Mie  Theory  for  the  case  r«A;  that 
is  the  radius  of  the  particle  is  much  smaller  than  the  wavelength 
of  the  incident  radiation.  The  basis  of  this  approximation  is  that 
the  contribution  of  the  efficiency  factor  Qsca  is  small  in  this 
region  so  that  the  contribution  to  the  total  scattering  is  small 
regardless  of  the  magnitude  of  n(x).  The  other  approximation  is  the 
Geometrical  approximation ,  valid  for  the  region  r>>A.  In  this  region, 
the  scattering  efficiency  factor  Qsca  approaches  a  constant  value  as 
shown  in  Figure  2*j. 

Description  of  Program 

The  computer  program  used  for  calculations  in  this  report 
was  originally  designed  by  Major  Paul  D.  Try  of  Headquarters,  Air 
Weather  Service  for  use  in  the  visible  spectral  region.  The  program 
was  converted  for  use  on  the  CDC  6600  at  Air  Force  Institute  of 
Technology  by  this  author  and  modified  for  use  in  the  near  infrared 
as  well  as  the  visible  spectral  region.  A  program  listing  is  not 
included  in  this  document.  A  complete  program  is  on  file,  however, 
at  the  Department  of  Physics  (ENP),  Air  Force  Institute  of  Technology, 
Wright -Patterson  Air  Force  Base,  Ohio.  A  short  description  of  this 
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program  follows. 

The  program  contains  two  subroutines.  One  of  the  subroutines 
does  exact  Hie  theory  calculations.  This  subroutine  was  originally 
developed  by  J.  V.  Dave  in  IBM  Report  No.  360D-17.2* .002  and  is  based 
on  Hie  theory  as  explained  in  Reference  62.  The  other  subroutine 
Is  an  approximation  method  to  the  exact  Mie  calculations  which  was 
programmed  by  this  author  based  on  articles  by  John  A.  Hodges  (Ref 
38:2310)  and  by  D.  Deirmendjian  (Ref  17*29,30).  This  approximation 
scheme  is  specifically  designed  to  limit  lengthy  machine  calculations. 
The  main  program  uses  the  Romberg  Integration  technique  to  integrate 
over  the  particle  size  distribution. 
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APPENDIX  C 
Use  of  the  Mode) 

In  order  to  determine  the  lock-on  range  of  the  laser  target 
designator  system  receiver,  the  geometry  of  the  problem  must  be  des¬ 
cribed.  Once  the  problem  geometry  is  described  then  solution  of  the 
problem  is  reduced  to  a  two-step  process:  (1)  determination  of  the 
total  atmospheric  attenuation  coefficients  based  on  atmospheric 
conditions  (observables),  and  (2)  substitution  of  these  coefficients 
into  the  laser  lock-on  range  equation 

Rr2  e°tr  =  (6.7  x  10V)e  "atd  (67) 

which  is  a  simplification  of  Eq  (12)  and  where  o  is  the  total 

tr 

optical  depth  (attenuation  coefficient  times  path  length)  from  target 

to  receiver  and  a  is  the  total  optical  depth  from  designator  to 
td 

target. 

The  first  step  is  the  subject  of  this  work,  and  therefore 
once  the  problem  is  reduced  to  Eq  (67)  the  problem  will  be  considered 
solved.  At  that  point,  the  remaining  transcendental  equation  may  be 
solved  with  the  aid  of  the  computer. 

Before  illustrating  the  use  of  the  model  the  following  important 
rules  for  determining  the  total  attenuation  coefficient  are  reiterated: 

(l)  In  rain,  the  total  attenuation  coefficient  is  the  sum 
of  the  molecular  absorption  coefficient,  the  aerosol  attenuation 
coefficient,  and  the  rain  attenuation  coefficient; 


CL!V.‘il/ 

(2)  In  fogs,  the  total  attenuation  coefficient  is  the  sum 
of  the  molecular  absorption  coefficient  and  the  fog  attenuation 
coefficient; 

(3)  In  clouds,  the  total  attenuation  coefficient  is  the  sum 
of  the  molecular  absorption  coefficient  and  the  cloud  attenuation 
coefficient;  and 

(k)  In  air  with  no  precipitation  falling,  the  total 
attenuation  coefficient  is  the  sum  of  the  molecular  absorption 
coefficient  and  the  aerosol  attenuation  coefficient. 

The  following  summary  of  important  tables  and  equations  is 
also  provided: 

(l)  Table  XXIV  provides  the  molecular  absorption 
coeff icientsj 

(2)  Table  XXVI II  provides  trie  aerosol  attenuation  coefficients 
for  the  background  air  mass; 

(3)  Eq  (A5)  provides  the  aerosol  attenuation  coefficient  for 
the  boundary  layer; 

(l{)  Table  XVI  provides  the  attenuation  coefficients  for  rain 
as  a  function  of  source  mechanisms  and  reported  rair.  category. 

(5)  The  attenuation  coefficient  of  fog  is  calculated  from 
Eq  (5&)  or  determined  from  Table  XVI J;  and 

(61  The  atteriuatior-  coefficient  of  clouds  is  found  in  Table 

XVII. 

Consider  now  the  following  example:  The  laser  designator 
is  at  an  altitude  of  2  km  at  a  look  angle  of  30°.  The  designator  is 
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Table  XXVIII 

Aerosol  Attenuation  Coefficients 
for  1.0636  Micron  Radiation 
for  Vertical  Clear  Air 
Profile  (Meteorological  Range 
of  40  KM  at  Ground  Level) 


Altitude 

Attenuation  Coefficient 

(km) 

(km-1) 

0 

.0523 

0  -  1 

.0*193 

1  -  2 

.0*135* 

2  -  3 

.0383 

3  -  *» 

.0338 

4  -  5 

.0299 

5-6 

.0263 

6  -  7 

.0232 

7  -  8 

.0205 

^Average  Attenuation  Coefficient  from  1.5  km  to  2  km  is  .0*121  km"1. 
However,  since  the  value  is  used  for  only  1/2  km  then  the  value  to 
be  used  in  the  vertical  summation  is  .0211  km. 
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illuminating  the  target  through  a  summer  afternoon  ha2e  where  the 
reported  visibility  is  5  km.  The  boundary  layer  extends  from  ground 
level  to  1.5  km.  The  receiver  is  located  at  an  altitude  of  6  km 
and  at  a  40°  look  angle  where  the  term  lock  angle  refers  to  the 
angle  between  the  vertical  (zenith)  and  line  of  vision  connecting 
target  to  receiver  or  designator. 

First,  let  us  calculate  the  atmospheric  transmission  for  the 

designator.  Table  XXIV  yields  the  vertical  molecular  absorption 

optical  path  (for  summer)  of  .04237.  Caution  should  be  exercised 

when  summing  coefficients  from  the  tables  so  that  the  zero  altitude 

value  is  not  included  in  the  sum.  Next,  the  boundary  layer  aerosol 

attenuation  ?s  given  by  Eq  (45)  as  .3165  km-1.  This  value  times  the 

% 

extent  of  the  boundary  layer  of  1.5  km  plus  the  average  clear  air 
value  from  1.5  km  to  2  km  yields  the  vertical  aerosol  optical  path 
of  0.4959.  Since  there  are  no  clouds  or  precipitation  then  the 
optical  depth  of  the  designator  becomes  .5382  times  secant  (30°). 
Then  the  right  hand  size  (R.H.S.)  of  the  laser  range  equation 
becomes 

R.H.S.  =  (6.7  x  109m2)e  '*6215  (68) 


The  receiver  optical  depth  is  computed  in  like  manner.  The  vertical 
molecular  absorption  optical  depth  is  .0985.  The  vertical  aerosol 
attenuation  optical  depth  is  .62415.  The  total  vertical  optical 
depth  is  then  .72305  and  the  total  receiver  optical  depth  becomes 
.72305  times  secant  (40°) .  The  left  hand  side  of  the  laser  range 
equation  becomes 
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Since  the  laser  range  equation  is  a  transcendental  equation  then 
these  values  must  be  checked  to  ascertain  if  this  is  a  solution  or 
if  the  receiver  is  outside  range.  This  can  best  be  done  with  the 
aid  of  a  computer. 


